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I.  Introduction 


All  chemical  reactions  take  place  at  definite  rates, 
arid  the  deteniiination  and  interpretation  of  these  rates 
has  long  been  one  of  the  central  problems  of  chemical 
research.  Fundamental  to  the  complete  treatment  of  this 
problem  is  an  understanding  of  the  mechanism  of  energy 
transfer  between  molecules  on  collision,  for  only  by 
learning  the  details  of  this  process  will  it  be  possible 
to  obtain  an  understanding  of  how  molecules  go  from 
their  normal  states  to  activated  states,  and  then  to  form 
products. 

This  thesis  is  concerned  with  the  study  of  very  rapid 

reactions  wherein  the  slowest  step,  and  hence  the  factor 

>diich  J-irdts  the  rate  of  approach  to  equilibrium,  may 

well  be  the  rate  of  energy  transfer  betv.-een  molecules  on 

collision.  For  this  purpose  the  isenLroplc  flow  experi- 

1 

ment  of  Kantrowitz  has  been  extended  for  the  case  of 
carbon  dioxide,  and  applied  to  the  study  of  the  pressure 
dependent  chemical  reaction  (NjO^;  ^==^2^02). 


A.  A Review  of  the  General  Theory  of  Energy  Transfer 
During  Molecular  Collision 

Early  studies  of  the  rate  of  energy  equilibration  be- 
tween the  Internal  and  external  degrees  of  freedor:i  of 

molecules  were  the  outgrowth  of  attempts  to  explain  the 

2 

absorption  of  sound  by  gases.  In  1870  Boltzmann  suggest- 
ed that  it  was  within  the  realm  of  possibility  that  the 
transfer  of  energy  between  molecular  motion  (translation) 
and  the  motion  of  the  atoms  within  a molecule  with  respect 
to  one  another  (vibration)  might  take  a rather  long  time. 
This  notion  was  expressed  in  a iouch  more  definite  form  by 

3 

H.  A.  Lorentz  in  a paper  published  initially  in  1331, 

/ 

and  entitled,  "Les  Equations  du  mouvement  des  gaz,  et  la 
propagation  du  son  suivant  la  theorie  cine'tique  des  gas". 
This  paper  contains  the  first  detailed  application  of  the 
then  new  kinetic  theory  of  gases  to  the  problem  of  sound 
propagation  through  a gas.  Application  of  the  law  of 
conservation  of  energy  to  the  system  gave  rise  to  a new 
term  vdiich  Lorentz  interpreted  as  being  due  to  a temporary 
displacement  on  rapid  compression  and  rarefaction  of  the 
ecuilibriui'i  between  the  translational  energy  of  the 
molecules  aaid  their  intramolecular  energy.  The  existence 
of  such  a delay  in  escablishing  equilibriiira  between  the 
external  and  internal  degrees  of  freedom  of  a molecule  is 
equivalent  to  a reduction  of  the  effective  heat  capacities 
of  the  gas.  Since  this  increases  the  effective  y,  the  net 
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result  is  an  increase  both  in  the  velocity  and  absorption 
of  sound  in  the  gas. 

A 

The  question  in  its  general  form  was  raised  by  Hertz 

5 

in  1332,  end  in  1903  Jeans  pointed  out  that  the  general 
solution  of  this  problem  is  not  obvious  even  for  a highly 
simplified  model  utilizing  classical  mechanics.  Jeans  in- 
vestigated the  effect  of  having  a portion  of  the  internal 
energy  lagging  behind  the  translational  energy,  and  showed 
that  both  the  velocity  and  absorption  would  be  increased; 
and  further,  that  the  new  absorption  would  be  proportional 
to  the  square  of  the  frequency.  Then,  using  as  a model  a 
sphere  with  its  center  of  mass  displaced  from  its  geometri- 
cal center,  he  attempted  to  calculate  the  magnitude  of  the 
effects  for  the  case  in  which  the  rotational  energy  lagged. 
Jeans  concluded  that  the  effect  would  be  small,  especially 
for  the  velocity,  but  that  the  observed  high  absorption  in 
air  could  be  accounted  for  in  this  fashion.  This  work  was 
undertaken  by  Jeans  in  an  attempt  to  reconcile  the  values 
of  Y found  by  sound  velocity  measurements  (the  Kundt-tube 
method)  with  those  calculated  from  the  principle  of  equi- 
partition  of  energy.  However,  the  quantum  theory  had  al- 
ready removea  the  major  discrepancy  between  the  measured 
heat  capacities  and  those  computed  from  equipartition,  and 
therefore  this  calculation  wr-s  dropped  from  Jeans'  "Dynajii- 
cal  Theory  of  G3l°s"  (Cc.Dibridge)  in  editions  subsequent  to 
the  first,  in  1904. 


4. 


No  new  theoretical  work  was  done  for  some  time,  but 

6 

in  1914  J.  G.  Stewart  published  a paper  entitled  ”The 
Inapplicability  of  Boltzmann's  Equipartition  Hypothesis  to 
Gases  in  a State  of  Change  of  Internal  Energy;  and  its 
bearing  on  the  experimental  determination  of  the  Specific 
Heat  of  Gases”.  He  pointed  out  that  in  experiments  on  the 
expansion  of  gases  through  nozzles  measurement  of  pressure 
or  temperature  after  a specified  expansion  gave  values 
lower  than  those  calculated  from  the  adiabatic  theorem 
using  the  equilibrium  values  of  y>  and  that  the  measure- 
ments of  amounts  of  material  discharged  give  values  great- 
er than  is  possible  with  reversible  expansion.  Also,  he 
noted  that  values  of  y taken  by  the  explosion  method, 
wherein  the  change  of  translational  energy  is  even  more 
rapid,  are  greater  than  those  obtained  by  static  methods. 
To  explain  these  discrepancies  he  invoked  the  hypothesis 
of  slow  energy  exchange  between  extcriial  ajid  internal 
modes.' 

The  next  great  impetus  to  the  discussion  of  this 
problem  was  the  invention  of  the  aco\:stic  interferometer 

7 

by  G.  W.  Pierce  in  1925.  For  the  first  time,  it  became 
possible  to  measure  the  velocity  of  sound  in  gases  with  an 
accuracy  of  0.1^.  For  carbon  dioxide  Pierce  found  a 
velocity  increase  of  several  per  cent  with  increasing  fre- 
quency; a value  far  in  excess  of  the  immeasurably  small 
increase  predicted  by  the  classical  theories  of  Stokes  and 
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Kirchoff.  He,  hiuself,  attributed  the  large  value  to  the 

effect  of  thermal  conductivity,  and  it  remained  for  Herz- 
0 

feld  and  Rice  to  give  the  correct  explanation  based  on 
the  slow  rate  of  exchange  of  energy  between  the  vibrational 

( 

and  translational  degrees  of  freedom. 

From  then  on  a wealth  of  sonic  data  became  available, 

and  coupled  with  this  advance  has  gone  an  ever  growing 

iinderstanding  of  the  process  of  energy  transfer  between  *• 

[i 

molecules  on  collision.  We  shall  not  attempt  here  to  make  ' 

! 

a complete  survey  of  rnis  field.  Rather  we  refer  those 
interested  in  the  utilization  of  sonic  phenomena  to  the  . 

9 

admirable  and  thorough  thesis  of  E.  H.  Freedman  , and 

those  interested  in  the  theoretical  aspects  of  relaxation 

phenomena  to  the  complete  discussion  presented  in  the 

1 o 

thesis  of  B.  Widom.  VJe  shall  site,  however,  some  of 
the  major  papers  which  are  landmarks  along  the  way  to  our 
present  understanding  of  the  process  of  energy  transfer, 
and  outline  briefly  the  conclusions  drawn  therein. 

The  first  lengthy  analysis  of  inelastic  molecular 

collision  processes  in  the  language  of  quantum  mechanics 

1 1 

was  published  by  C.  Zener  in  1935.  As  a first  approxi- 
mation his  model  pictured  each  molecular  system  of  atoms 
as  bound  together  by  elastic  forces,  while  atoms  of  dif- 
ferent molecules  interacted  as  hard  elastic  spheres.  Con- 
sidering the  one  dimensional  case,  he  then  calculated  the 
energy  interchange  for  a single  head  on  impact  from  the 
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principles  of  energy  and  raonentum  conservation.  It  is 
clear  that  the  finite  tine  required  to  establish  contact 
between  the  two  systems  .acts  os  a severe  restraint  upon 
the  exchange  of  internal  energy.  The  interaction  between 
colliding  molecules  is  actually  so  gentle  that  in  general 
many  vibrations  occur  during  the  effective  time  of  a col- 
lision, and  theory  shows  that  the  energy  transfer  dies  off 
exponentially  with  the  number  of  these  vibrations  which 
can  take  place  during  such  an  interval.  Zener  then  im- 
proved the  rigid  sphere  model  by  assuming  that  the  atoms 

of  separate  systems  repel  one  another  according  to  the  . 

-r/d 

simple  law  V = ce  . Classical  mechanics  applied  to 
this  model  shows  that  for  d not  too  small  there  is  but  a 
small  tendency  for  the  energy  of  a rapid  oscillator  to  be 
transferred  to  a less  rapid  oscillator  or  to  translational 
degrees  of  freedom.  This  treatment  points  up  the  fact  that 
a specified  amount  of  internal  energy  will  bo  dissipated 
much  more  rapidly  if  it  can  be  given  up  to  translational 
energy  as  several  small  quanta  rather  than  as  one  large 
quantum.  Hence  for  polyatomic  molecules  we  ’.rould  expect, 
all  other  things  being  equal,  that  t .osc  vibrations  with 
lower  frequencies  will  tend  to  have  shorter  relaxation  times. 
An  interesting  sidelight  is  Zener's  conclusion  that  the  com- 
paratively large  quanta  which  Hz  can  absorb  as  energy  of 
rotation,  which  are  still  small  compared  to  vibrational 
quanta,  makes  it  unique  in  deactivating  efficiency. 


7. 


In  principle  the  calculation  of  the  extent  of  energy 

transfer  on  collision  recuires  the  use  of  a coriolete  cuantura 

1 2 

mechanical  picture.  However,  Zener  and,  independently, 

1 3 

Lsndau  and  Teller  have  proposed  an  approximate  semi- 
classical  theory  of  these  inelastic  molecular  collision 
processes  which  glosses  over  a number  of  the  difficulties 
encountered  in  the  quantum  mechanical  problem.  In  this 
semi-classical  approach  the  vector  which  gives  the  sepa- 
ration of  the  centers  of  mass  of  the  two  particles  is  taken 
to  be  an  explicit  function  of  the  time  as  given  by  classical 
dynamics.  The  interaction  potential  is  then  a known 
function  of  the  time,  and  its  effect  on  the  internal  states 
of  the  system  can  be  found  by  application  of  the  first 
order  time-dependent  perturbation  theory.  Using  this  ap- 

1 3 

proach  Landau  and  Teller  have  published  a general  theory 
for  the  relaxation  time  of  a real  non-dissociating  poly- 
atomic gas  wherein  several  modes  of  vibration  may  contribnite 
to  the  energy  absorbing  process.  They  established  that 
large  velocities  play  a dominating  role  and  predicted  that 

the  logarithm  of  the  relaxation  time  of  a gas  expressed  in 

-1/3  * 

molecular  collisions  should  be  rroportional  to  T . Re- 
cently a thorough  discussion  of  the  relationship  between 
the  completely  quantum  mechanical  and  semi-classical 


See  section  I.  E.;  cf.  ref.  1(a) iii,  figure  6. 
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theories  has  been  presented  by  Uidom. 

A rather  complete  sumi.iary  of  work  on  the  problem  of  i 

molecular  ener^jy  transfer  up  until  1949  is  provided  by  I 

1 4 

H.  S.  \i.  Massey’s  survey  article  entitled  "Collisions 
Between  Atoms  and  Molecules  at  Ordinary  Temperatures”.  i 

He  sums  up  the  present  status  of  the  problem  as  follows:  | 

”A  considerable  volume  of  evidence  is  available  con-  j 

'i 

cerning  the  ease  with  which  energy  exchange  may  occur  bo-  ; 

tween  that  of  relative  trcinslation  and  various  degrees  of 
internal  freedom  of  the  colliding  systems,  nevertheless, 
although  it  is  possible  to  give  general  rules  concerning 
these  energy  transformations,  we  are  still  far  from  possess- 
ing an  interpretation  so  complete  a.s  to  raalce  possible  the  ‘ 

prediction  of  the  rates  of  a particular  reaction  excent  in  | 

a very  few  cases.” 

1 5 

In  a recent  paper  Curtis  and  Adler  have  applied 
grouo  theoretical  arguments  to  the  treatment  of  three- 

! 

dimensional  collisions  following  Zener's  work  for  the  one  j 

’ ’ 

dimensional  case.  The  scattering  coefficients  were  ob-  f 

tained  by  an  analysis  of  the  solutions  of  the  internal  wave 

1 

equations  involving  the  three  coordinates  specifying  the 
internal  configuration.  The  angular  transition  probabili- 
ties were  obtained  for  the  case  in  >/hich  transitions  are 
relatively  infrequent  (weak  couoling)  in  terms  of  integrals 
of  the  interaction  potential  and  a set  of  phase  shifts. 

This  treatment  is  applicable  to  normal  low  energy  thermal 
collisions,  but  not  to  the  high  energy  collisions  which 
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occur  in  hot  atom  chemistry.  Since  they  consider  the 
question  of  energy  transfer,  molecular  dissociation  or 
rearrangements  are  not  directly  part  of  this  problem  and 
hence  were  not  considered.  The  essential  difference  is 
that  in  the  energy  transfer  process  only  weak  coupling  is 
postulated  between  colliding  \inits. 

1 6 

In  a brief  recent  paper  Herzfeld  has  succeeded  in 
calculating  directly  the  rate  of  energy  exchange  between 
internal  and  external  degrees  of  freedom  for  liquid  benzene. 
The  Lennard-Joncs  interaction  potential  and  model  of  the 
liouid  are  used;  the  fourth-order  terms  in  the  interaction 
provide  the  coupling  between  internal  vibrations  and  Debye 
waves.  This  results  in  a value  three  times  too  rapid, 
vrtiich  is  in  astonishingly  good  agreement  with  the  experi- 
ment, in  view  of  the  approximations  made, 

1 7 

Yet  more  recently  Schwartz,  Slawsky,  and  Kerzfeld 
have  attempted  a cuantitative  calculation  of  the  relaxation 
time  in  pure  gases  and  in  chemically  nonreoctive  gas  mix- 
tures. They  followed  the  wholly  quantum  mechanical  :iethod 
1 ^ 

of  Zener  , using  an  exponential  repulsion  in  a one-dimen- 
sional model,  and  concluded  that  in  some  cases  the  great 
effect  which  impurities  have  may  be  accounted  for  either  by 
their  low  masses  and  resultctnt  high  velocity,  or  by  ’’near 
resonance"  transfers  in  which  the  vibrational  cusn-um  of  the 
substratum  is  used  partly  to  excite  the  vibration  of  the 
impurity,  only  the  difference  being  transferred  to  trens- 
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lation.  However,  there  are  a number  of  interesting  cases 
for  which  neither  of  these  explanations  apply.  The  theo- 
retical values  here  obtained  for  the  relaxation  times  are 
10  to  30  times  shorter  than  the  experimental  values,  but 
the  authors  express  the  opinion  that  this  difference  may 

have  arisen  because  of  the  use  of  a one-dimensional  model. 

1 o 

Vidom  has  extended  2ener»  s original  proposal  for 
S3nsmetrizing  the  theory  in  the  case  of  head-on  collisions 
to  cover  the  more  general  case  of  non- zero  impact  parameter 
and  has  provided  proof  for  the  assertions  of  Rosen  and 

IS  13' 

Zener  and  of  Landau  and  Teller,  to  the  effect  that  the 
semi-classical  theory  is  an  adequate  approximation  to  the 
quantum  mechanical  theory  whenever  the  relative  trans- 
lational energies  -v/hich  are  involved  are  much  greater  than 
the  amount  of  energy  exchanged.  He  discussed  the  statis- 
tics of  semi-classical  collisions  in  order  to  derive  ex- 
pressions for  temperature-dependent  probabilities  of  inel- 
astic encounters  and  presented  a detailed  discussion  of 
classical  motion  in  a spherically  symmetric  inter-molecular 
potential.  The  system  CO2-H2O  was  treated  semi-classically, 
assuming  that  the  interaction  potential  was  spherically 
S3anmetric,  and  corresponded  to  that  interaction  between  the 
CO2  and  H2O  molecules  which  would  occur  were  they  to  approach 
in  a certain  most  favorable  orientation.  The  resulting 
probabilities  were  then  corrected  by  multiplying  by  a suit- 
abl<=>  orientation  factor.  The  relevant  potentials  were  based 
on  considerations  of  the  chemistry  of  the  acid- base  inter- 
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action  between  CO 2 and  H2O.  A potential  of  the  type  (a) 
was  used  rather  than  type  (b),  which  is  that  for  the  inter- 


(a)  (b) 


action  of  molecules  between  which  there  can  be  no  chemical 
interaction.  For  a molecule  with  a given  collisional 
energy  (dotted  arrow)  it  is  clear  that  for  (a)  the  minimum 
distance  of  approach  is  of  the  order  of  a bonding  distance 
(r_),  and  that  this  interaction  occurs  nore  sharply  (great- 
er  slope")  than  fc’'  b).  Thus  in  (a)  the  effective  time  of 
collision  is  less,  and  hence  there  is  an  increased  proba- 
bility of  energy  transfer.  On  the  basis  of  the  s’emi- 
classical  treatment,  and  the  statistics  of  generalized 
impact  parameters,  Widom  found  that  the  effective  collision 
diameter  is  a monotonically  decreasing  function  of  the 
temperature.  He  then  showed  that  when  corrected  for  the 
low  energy  failure  of  the  semi-classical  theory,  the  col- 
lision diameter  also  decreases  at  low  temperatures,  and 
therefore,  in  agreement  with  experiment,  it  passes  through 
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a maximum.  This  theory  quantitatively  accounts  for  the  ob- 

1 9 

servation  of  Eucken  and  Becker  that  the  probability  of 
transfer  of  energy  on  collision  is  larger  the  greater  the 
chemical  reactivity  of  the  two  colliding  molecules,  even 
though  no  chemical  change  actually  takes  place. 

20 

The  transition  state  approach  of  Laidler  and  Eyring 
has  been  of  some  use  in  correlating  qualitatively  certain 
general  observations.  The  above  mentioned  observation  of 

1 9 

Eucken  and  Becker  was  qualitatively  explained  as  follows. 

Consider  a potential  energy  surface  which  represents  a 

collision  between  X and  Y ■vrtiich  could  react  chemically. 

However,  in  an  energy  transfer  process  we  are  concerned  not 

with  major  transitions  between  valleys  in  the  potential 

energy  surface,  vdiich  occur  when  a chemical  change  takes 

place,  but  rather  with  transitions  between  the  vibrational 

levels  within  one  valley.  For  proper  presentation  of  the 

potential,  the  surface  should  be  constructed  with  skewed 

2 1,22 

axes  , and  the  degree  of  dissymmetry  of  the  valley 
will  depend  on  the  tendency  for  chemical  reaction.  V/hen 
a chemical  reaction  is  possible  the  valley  is  definitely 
not  symmetrical  but  distorted,  particularly  at  positions 
corresponding  to  close  approach  of  the  particles,  and  when 
the  system  enters  this  distorted  region  the  probability  of 
transfer  to  another  vibrational  lovel  becomes  appreciable, 
and  may  be  high.  Thus  the  transference  of  energy  will  be 
associated  in  an  approximate  manner  with  chemical  reactivity. 


It  is  clear,  however,  that  the  number  of  possible  initial 
and  final  states  may  bo  very  large  and  that  even  between 
two  given  states  the  number  of  likely  reaction  path.^  :.u-y  be 
very  great.  It  is  this  complexity  which  makes  any  detailed 
microscopic  treatment  quite  difficult,  and  predictions 
nearly  impossible. 

As  an  extension  of  the  potential  energy  surface  ap- 
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proach.  Castellan  and  Hulbert  have  calculated  reflection 
coefficients  for  a matter  wave  in  a parabolic  channel  and 
in  a square  channel,  each  channel  closed  at  one  end  by  a 
vertical  plane  tilted  slightly  from  the  normal  to  the  axis 
of  the  channel.  This  enabled  tho-j'  to  find  theoretically 
the  number  of  collisions  necessary  to  dissipate  one  quantum 
of  vibrational  energy  end  comparison  with  the  data  from 
measurements  on  sound  dispersion  shows  good  agreement  for 
N;0  and  COj.  They  also  arrived  at  some  tentative  con- 
clusions regarding  t.,>e  teituperature  dependence  or  the 
average  probability  of  tran.sition  per  collision. 

The  correlation  of  the  probability  of  colli sional 
energy  transfer  with  chemical  reactivity  is  not  strong  in 
the  dispersion  of  souna  experiments,  since,  owiap  to  the 
low  temperatures,  the  vibrational  levels  are  far  down  the 
potential-energy  valleys,  in  which  region  the  distortion 
due  to  reactivity  is  very  small.  Hence  the  number  of  col- 
lisions required  to  effect  a change  of  level  is  very  Icrge. 
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In  the  spectroscopic  work  on  energy  level  trans:, tions. 


on  the  other  hand,  the  vibrational  levels  involved  are 
higher,  their  spacing  is  much  smaller,  and  the  number  of 
collisions  required  is  therefore  less;  the  correlation 
with  chemical  reactivity  is  much  more  pronounced. 

With  respect  to  unimolecular  reaction,  it  seems  clear 
from  the  experimental  observation  that  no  amount  of  t.Jded 
gas  can  increase  the  specific  rate  to  a value  greater  than 
its  limiting  value  at  high  pressures,  that  the  added  mole- 
cules do  not  interfere  with  the  over-all  chemical  reaction, 
end  that  their  function  is  to  maintain  the  equilibrium  con- 
centration of  activated  molecules.  If  the  property  possess- 
ed by  foreign  gases  of  maintaining  the  equilibrium  concen- 
tration of  activated  complexes,  and  thus  preventing  the 
falling  off  in  rate  of  a unimolecular  reaction  at  low  pres- 
sures, is  to  be  attributed  to  their  ability  to  transfer 
some  of  their  energy  to  vibrational  energy  of  the  reacting 
molecule,  then  some  parallelism  is  to  be  expected  between 
their  efficiency  in  this  respect  and  the  removal  of  excess 
vibrational  energy  in  association  reactions.  This  paral- 
lelism is  borne  out  by  the  experimental  evidence. 
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Rosen  discussed  the  case  of  collision  between  a 
molecule  and  an  atom  or  another  molecule  >rtiere  some  vi- 
brations are  excited  in  the  former  and  a binding  to  the 
latter  is  produced.  A stable  system  may  form  after  col- 
lision only  if  the  potential  energy  of  the  particles  as  a 
function  of  the  distance  between  them  has  a minimum  and  if 
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the  sum  of  this  potential  energy  and  the  energy  of  relative 
motion  is  less  than  the  top  of  either  side  of  the  potential 
valley.  Clearly  the  free  particler.  have  energy  above  the 
outer  wall  of  the  potential  valley,  and  Rosen  shows,  by 
means  of  an  approximate  solution  of  the  wave  equation, 
that  in  some  cases  a portion  of  this  energy  of  relative 
motion  may  flow  to  other  parts  of  the  system. 

For  the  purpose  of  determining  the  order  of  magnitude 
of  the  lifetime  of  a molecule  having  an  excess  of  energy 
Rosen  supposed  that  there  are  three  atoms  in  a line,  the 
two  outside  atoms  each  being  connected  to  the  inside  atom 
by  a potential  of  the  :'lorse  tyte.  Ke  then  considered  cne 
case  v^ere  one  vibrator  in  a relatively  low  state,  with 
vibrational  quantum  number  m,  is  supposed  to  drop  bacK  to 
the  lowest  state  (vibrational  quantum  maaber  zero),  the 
energy  being  taken  up  in  the  transition  into  the  continuum 
representing  dissociation  of  the  other  vibrator.  Rosen  con- 
cluded that  if  m > 4 there  is  very  little  chance  of  the 

molecule  decomposing  before  it  exchanges  energy  at  collision. 
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Rice  pointed  out,  however,  that  it  is  not  necessary 
that  the  exchange  of  energy  should  take  place  in  one  step, 
but  that  it  can  take  place  throu^^h  a series  of  exchanges  of 
smaller  amounts  of  energy  between  the  vibrators  without  dis- 
sociation taking  place,  until  one  of  thorn  is  in  a state  where 
its  probability  of  dissociation  is  great.  The  probability 
of  dissociation  ultimately  taking  place  will  then  not  fall 
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off  as  rapidly  with  increasing  m as  Rosen  has  predicted. 

2 7 

Amongst  the  most  recent  papers  are  two  by  Bauer 
in  which  he  discussed  the  recombination  of  two  atoms  in 
the  presence  of  a third  body  by  considering  the  probability 
coefficient  of  the  reverse  process,  the  dissociation  of  a 
molecule  by  collision.  Various  mechanisras  for  dissociation 
were  briefly  discussed,  and  an  attempt  wis  made  to  estimate 
the  general  order  of  magnitude  of  the  cross  section  for  the 
excitation  of  a molecule  into  the  vibrational  continuum  of 
its  electronic  ground  state.  This  cross  section,  which  was 
calculated  by  a method  depending  upon  correcting  the  Born 
approximation  by  means  of  the  exact  one  dimensional  prob- 
lem, was  computed  for  the  recombinatioxi  of  two  hyd’^s^ren 
atoms  in  the  presence  of  a third  during  a.  collision.  The 
three-body  coefficient,  obtained  from  the  cross  section  of 
the  reverse  dissociation  process  by  the  principle  of  detail- 
ed balance,  is  much  smaller  than  the  empirical  values.  An 
estimate  of  the  number  of  thermal  collisions  required  to 
de-excite  the  first  vibrational  state  of  hydrogen  by  col- 
lisions with  protons  is  More  nearly  in  line  with  empirical 
values. 

2 8 

Messiah  has  extended  a mass  tensor  approximation  to 
the  computation  of  the  scattering  of  slow  neutrons  by 
molecules  taking  Into  account  the  molecular  vibrations. 

The  calculations  were  made  for  the  case  of  and  CH^  at 
room  temperature  for  neutron  energies  below  the  threshold 
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of  excitation  of  higher  vibrational  states.  Ihe  agreement 
with  experiment  is  very  satisfactory  in  the  expected  range 

of  validity.  The  asymptotic  form  of  the  total  cross  section 

* 

for  large  values  of  the  incident  energy  was  derived,  and 
numerical  calculations  show  that  the  asymptotic  behavior 
is  quickly  reached. 
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B.  Measurements  on  Carbon  Dioxide 

,1 

a.  Sonic  Method  * 

I 

As  was  noted  in  the  preceding  section,  Pierce  was  the 
first  investigator  to  make  precision  measurements  of  the 
velocity  of  sound  in  various  gases,  and  the  increase  of  | 

velocity  with  increasing  frequency  which  he  observed  was 
correctly  interpreted  shortly  thereafter  by  Herzfeld  and 
Rice  as  resulting  from  the  finite  time  required  for  the 
equilibration  of  energy  between  vibrational  and  trans- 
lational degrees  of  freedom. 

Since  then  very  many  measurements  have  been  made  and  I 

} 

a diversity  of  methods  for  taking  sound  velocity  and  ab-  j 

i 

j 

sorption  data  invented.  We  will  malce  no  attempt  to  sum-  i 

marize  these  investigations,  but  rather  refer  the  reader  ’ 

I 

to  three  papers  >diich  together  admirably  cover  the  period 
from  the  investigations  of  Pierce  up  until  relatively  recent 
times.  - 

The  first  of  these  papers  is  a review  published  by 
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Richards  in  1939  and  entitled  ’’Supersonic  Phenomena”. 

This  paper  opens  with  a complete  derivation  from  first 
principles  of  the  equations  governing  the  propagation  of 
plane  sound  waves,  and  then  proceeds  to  review  experimental 
methods  and  observations.  The  coverage  of  this  paper  is 
quite  remarkable  as  may  be  suggested  by  the  fact  that  its 
bibliography  lists,  with  at  least  a short  comment  on  each, 
no  less  than  three-hundred  forty-eight  entries. 
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The  second  paper  is  by  Kittel  and  is  entitled  "Ultra- 
sonic Research  and  the  Properties  of  I'latter".  This  excel- 
lent review  paper  covers  propagation  in  gases,  liquids,  and 
solids,  experimental  methods,  applications  to  other  fields, 
and  unusual  phenomena  associated  with  liquid  h.elium.  The 
bibliography  concentrates  on  papers  since  194.1,  end  thus 
forms  with  that  of  Richards  a very  nearly  complete  listing 
of  the  pertinent  studies  of  ultrasonic  phenomena  published 
prior  to  194^* 

3 1 

The  third  review  is  by  Markham,  Beyer,  and  Lindsay 
and  was  published  in  1951.  This  paper  is  entitled  "Absorp- 
tion of  Sound  in  Fluids",  a broad  field  which  it  covers 
quite  well.  Emphasis  was  placed  on  a discussion  of  basic 
principles  and  concepts  rather  than  on  detailed  calculations. 
A listing  of  the  more  reliable  experimental  data  was  in- 
cluded. The  excellent  bibliography  covers  important  work 
in  the  field  published  prior  to  1951. 

3 2 

In  1950  Laiubert  and  Rowlinson  published  a very  in- 
teresting study  in  wiiich  they  reported  measurements  of 
ultrasonic  dispersion  in  a series  of  organic  vapors  over  a 
range  of  frequencies,  temperatures,  and  pressures.  They 
concluded  that,  for  polyatomic  molecules  in  general,  there 
is  little  hindrance  to  the  interconversion  of  translational 
and  vibrational  energy,  especially  at  high  temperatures. 

Slow  Interconversion,  leading  to  ultrasonic  dispersion,  was 
found  only  for  small  or  rigid  molecules,  where  internal 
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rotations  or  vibrations  of  low  frequency  are  absent. 

33 

Recently  Zartman  has  been  able  to  improve  the 
ultrasonic  interferometer  so  as  to  attain  greater  sensi- 
tivity and  reproducibility.  These  improvements,  coupled 
with  scrupulous  care  in  preparing  samples,  should  go  far 
towards  increasing  the  consistency  of  measured  values.  It 
seems,  however,  that  this  technique  has  very  nearly  reached 
the  upper  limit  of  its  precision. 

Out  of  this  wealth  of  work  the  following  basic  conclu- 
sions fundamental  to  this  thesis  have  been  reached  with 
respect  to  dispersion  and  absorption  of  sound  in  pure  gases: 
i)  Dispersion  and  absorption  of  soxind  are  to  be  expect- 
‘ ed  in  any  polyatomic  gas  in  ’idiich  there  is  an— 
appreciable  vibrational  excitation,  due  to  the 
inevitable  inefficiency  of  the  conversion  of  trans- 
lational energy  into  quanta  of  vibrational  energy 
by  collision. 

ii)  The  effect  of  pressure  variations  is  in  accord  with 
the  kinetic  theory  of  heat;  an  increase  of  pressure 
increases  the  number  of  collisions  and  displaces 
the  dispersive  region  to  higher  frequencies.  Double 
collisions  appear  to  be  the  usual  mechanism  by 
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which  vibrational  energy  is  excited  ; all  of  the 
measurements  with  carbon  dioxide  have  indicated 
that  the  relaxation  time  is  proportional  to  the 
pressure  of  the  gas  indicating  a bimolecular  re- 
laxation process. 


. ... 
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Hi)  Generally  the  effect  of  increasing  the  temperature 
is  to  increase  the  tremsition  probability  (note 
exception,  H^O  with  COj). 

iv)  Very  small  traces  of  foreign  gases  may.  have  an 

enormous  effect  on  the  efficiency  of  the  collision 
process,  and  consequently  on  the  location  of  the 
dispersive  region  (see  later  data  on  carbon  dioxide). 
The  effect  of  the  presence  of  these  traces  is  always 
to  shift  the  dispersive  region  to  higher  frequencies; 
any  other  result  would  necessitate  a revision  of  the 
collision  theory. 

v)  In  nearly  all  cases  sonic  measurements  show  several 
energy  states  adjusting  with  a common  frequency 
over  the  range  measured.  That  is,  experimentally 
one  obtains  only  one  dispersive  region.  It  has 
long  been  considered  unlikely  that  the  several  states 
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have  equal  relaxation  times  in  these  cases  , but 
rather  it  has  seemed  more  probable  that  they  fall 
out  together  because  of  rapid  redistribution  of 
energy  within  the  molecule.  This  point  has  now, 
however,  been  resolved  by  a new  method  discussed 
[ ^ later. 

Over  three  score  studies  on  the  velocity  and  absorption 
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of  sound  in  carbon  dioxide  are  reported  in  the  literature  ; 

37 

those  of  Eucken  and  co-workers  on  carefully  purified 
carbon  dioxide  cover  the  widest  ranges  of  temperatures. 


frequencies,  and  pressures  and  show  the  greatest  over  all 
prev,ision  and  consistency.  They  have  also  made  detailed 
studies  of  the  effects  produced  by  -tnown  anounts  of  foreign 
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gases.  Me  shall  use  their  data  for  comparison  purposes, 
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and  refer  the  reader  to  the  paper  by  Kantrowitz  for  a 
comparison  with  the  data  of  other  workers.  Following  the 
formula  of  Landau  and  Teller  we  here  plotted  in  figure  1 
the  log  of  the  relaxation  time  (in  molecular  collisions) 
versus 

At  293°K.  the  '’collision  number"  for  pure  carbon  di- 
oxide is  57,000,  while  between  water  and  carbon  dioxide  it 
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is  105.  Hence  collisions  between  carbon  dioxide  and  water 

molecules  are  over  500  times  as  effective  in  de-exciting 

the  carbon  dioxide  as  are  collisions  between  carbon  dioxide 

molecules.  For  details  on  the  effect  of  changing  the  tem- 

; o 

perature  see  work  by  V/idom. 
h.  OhTick  TnbA  Method 

40 

Recently  Smiley,  Winkler,  and  Slawsky  have  reported 
the  results  of  shock  tube  investigations  on  carbon  dioxide. 
They  recorded  by  means  of  interferograms  the  relaxation 
region  behind  plane  shock  waves  in  a shock  tube  containing 
carbon  dioxide.  In  these  interferograms  the  fringe  dis- 
placement is  a direct  measure  of  the  density  change.  The 
peak  temperature  in  the  carbon  dioxide  shock  is  about  900°K., 
at  \diich  temperature  the  lowest  vibrational  modes  of  carbon 
dioxide  make  a large  contribution  to  the  enthalpy.  They 


Figure  1 


Plot  of  the  logaritlim  of  the 
Relaxation  Tine  (in  noiecular 
co.Tlisions)  ^ T“’/^for  Carbon 
Dioxide 
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obtained  a value  of  approximately  0.1^isec.  at  normal  pres- 
sures, whereas  a value  of  6psec.  would  be  expected  for 
bone  dry  carbon  dioxide.  They  claim,  however,  that  the  ob- 
served value  is  in  rough  agreement  with  that  expected  for 
the  known  water  vapor  content  of  about  3f5.  The  brief  note 
published  indicates  that  this  work  will  be  extended;  it  is 
hoped  that  attendant  upon  this  extension  greater  accuracy 
in  specifying  the  experimental  conditions  will  be  forth- 
coming. Specifically,  more  definite  information  with  re- 
spect to  purity  of  reagents  and  precision  of  the  readings 
would  be  desirable,  as  well  as  studies  to  determine  whether 
the  relaxation  rate  is  really  proportional  to  the  displace- 
ment for  such  a large  disturbance  from  equilibrium. 

c.  Infrared  Method 

The  newest  method  to  be  applied  to  the  study  of  relaxa- 
tion phenomena  in  pure  and  mixed  gases  involves  a process 
which  may  be  looked  upon  as  the  reciprocal  of  the  sound 
dispersion  method.  In  this  new  method,  modulated  radiation 
of  such  a spectral  composition  as  to  be  absorbed  by  one  or 
several  degrees  of  freedom  is  passed  through  the  gas.  This 
energy  is  absorbed  through  vibrational  excitation,  and  if 
there  were  immediate  equilibration  between  all  degrees  of 
freedom  it  would  appear  instantly  as  a rise  in  translational 
temperature,  and  hence  in  pressure.  However,  if  there  Is  a 
lag  in  the  establishment  of  equilibrium  there  will  be  a lag 
between  irradiation  and  rise  in  trcnslational  temperature. 


The  pressure  variations  in  the  gas  are  followed  by  a micro- 
phone, and  by  studying  the  phase  relation  between  the  micro- 
phone output  and  the  impinging  modulated  beam  one  can  deter- 
mine the  energy  transfer  lag  between  those  nodes  which 
absorb  the  radiation  and  those  of  trar\3lation,  responsible 
for  the  pressure.  This  method  was  suggested  independently 
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by  Gorelik  and  Slobodskaya  in  Russia  and  Cottrell  in 
England. 

The  most  interesting  feature  of  this  method  is  that 

it  enables  one  to  determine  the  relaxation  time  of  individ- 

* 

ual  vibrational,  or  even  rotational,  inodes.  Slobodskaya 
applied  this  method  to  the  much  disputed  case  of  carbon 
dioxide,  and  has  measured  the  lifetimes  of  vibrational 
quanta  corresponding  to  the  three  absorption  bands  2.7p, 
4.3p,  and  14.Su. 


A(n) 

T (sec. ) 

Probability  of 
Inelastic  CollisicTS 

4.3 

7 x 10-^ 

1.6  X 10-’ 

2.7 

4 x 10-*^ 

2.S  X 10"’ 

14.8 

1.6  X 10“* 

6.9  X 10“’ 

In  the  case  of  sound  dispersion  measurements,  which 
give  only  a single  relaxation  time,  we  would  expect  to  get 
a value  equal  to  or  slightly  larger  than  that  of  the  lowest 
vibrational  frequency.  The  value  of  1.6  x 10~*  sec.  for 
the  14. Sp  band  obtained  In  this  work  agrees  well  with  that 
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See  comment  in  section  I.D.a. 


of  Kneser  who  obtained  1x10-“  sec.  from  sonic  measurements 
under  the  same  conditions. 

Another  great  advantage  of  this  method  is  the  very 
small  soUTiple  required  and  the  ease  of  isolating  it  from 
contamination.  While  Slobodskaya  reports  a precision  of 
10^  within  a series  of  measurements,  it  is  probable  that 
the  accuracy  can  be  increased  with  further  work.  To  cate 
not  enough  is  known  about  this  method  to  indicate  its 
ultimate  possibilities,  but  it  is  clear  that  they  must 
certainly  be  quite  great. 

d.  Impact  Tube  Method 

Since  this  method  is  the  one  utilized  in  this  thesis 
it  will  be  described  in  detail  elsewhere.  V/e  here  quote, 
however,  the  results  obtained  by  two  previous  investigators 
who  have  applied  this  method  to  the  measurement  of  the  re- 
laxation time  of  carbon  dioxide.  Kantrowitz^^^^ obtain- 
ed as  a final  average  collision  niumber  3<,600  at  105®F.  At 
a temperature  of  65°F.  and  one  atm.  Griffith  obtained  a 
value  of  2.02n  sec.  for  the  relaxation  time  of  carbon 
dioxide  containing  an  unknown  concentration  of  HjO. 
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C.  Measurements  on  Nitrogen  Tetroxide 
a.  Sonic  Method 

In  1910  Nernst  and  Keutel  first  suggested  that  it 
might  be  possible  to  determine  the  velocity  of  a chemical 
reaction  by  measuring  the  velocity  of  sound  through  an 
equilibrium  mixture.  Using  high  enough  frequencies,  one 
should  obtain  the  velocity  corresponding  to  the  passage  of 
the  sound  wave  through  the  mixture  without  reaction,  but 
at  sufficiently  low  frequencies  one  should  expect  the  equi- 
librium to  adjust  itself  by  reaction  to  the  periodic  com- 
pression and  rarefaction,  and  hence  decrease  the  velocity 
of  sound.  At  intermediate  frequencies  one  would  then  ex- 
pect a "critical  region"  with  intemediate  velocities. 

In  his  Berlin  Inaugural-Dissertation  Keutel  reported 
the  results  of  applying  the  Kundt-tube  method  of  measuring 
sound  velocities  to  the  NO2-N2O4  system,  which  he  studied 
at  various  pressures.  He  concluded  that  the  gas  does  not 
behave  as  a mixture  of  two  non-reacting  ases,  but  that  the 
reaction  between  them  was  not  so  fast  as  to  maintain  com- 
plete chemical  equilibrium  during  the  passage  of  the  sound 
wave.  He  also  concluded  that  the  specific  heat  cannot  be 
measured  by  this  method  because  of  the  displacement  of  the 
equilibrium  is  unknown,  but  surmised  that  with  improved 
technique  and  r.ore  accurate  thermodynamic  data  it  would 
prove  possible  to  determine  the  rate.  He  did  not,  however, 
derive  the  relationship  between  the  sound  velocity  and  the 
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rate  constant. 
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Einstein,  in  192C,  developed  the  theory  of  the  effect 
of  freouency  on  the  velocity  and  absorption  of  sound  in  a 
mixture  of  chetaically  reacting  ideal  gases  of  negligible 
absorption. 

The  dissociation  of  the  colorless  gas  nitrogen 
tetroxide  to  form  the  highly  colored  red  nitrogen  dioxide 
provided  one  of  the  earliest  examples  for  the  application 
of  thermodynamics  to  the  calculation  of  the  effect  of  toin- 
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perature  on  equilibrium.  It  was  also,  as  we  have  noted, 
the  first  chemical  reaction  to  which  the  sonic  method  of 
rate  determination  was  applied.  There  have  been  many  sub- 
sequent attempts  to  utilize  the  sonic  method  to  determine 
the  rate  of  this  reaction,  both  because  the  reectjon  is  it- 
self of  intrinsic  interest  due  to  the  wealth  of  information 
about  its  other  parameters,  and  because  the  low  hiat  capacity 
of  the  reactants  and  products  compared  to  the  high  heat  of 
reaction  make  the  chances  of  success  particularly  favor- 
able.  The  maximum  decrease  in  sound  velocity  has  been 
calculated  for  this  case  to  be  approximately  3.3ji  at  25®C. 

4 9 

and  1 atm.  This  theoretical  calculation  is  for  the  low 
frequency  limit  vhere  the  chemical  reaction  follows  com- 
pletely. We  shall  now  review  briefly  the  sonic  measurements 
which  have  been  reported  on  this  system. 

The  first  measurements  made  on  this  system  were  those 
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of  E.  and  L.  Natanson,  who  measured  the  velocity  of  sound 
for  only  a single  frequency,  and  hence  could  make  no  esti- 


mate  as  to  the  position  of  the  critical  range  of  frequencies 
where  the  velocity  would  be  affected  by  the  reaction.  Keutel, 

5 1 

and  later  Selle,  concluded  that  the  "critical  region"  lay 
within  the  audible  range;  a conclusion  not  borne  out  by 

5 2 

later  work.  The  early  experiments  of  Argo  , snd  later 
those  of  Grilneisen  and  Goens  demonstrated  that  the  criti- 
cal range  lies  above  the  audible  range  of  frequency. 
Einstein’s  equations  were  first  employed  by  Richards  and 
Reid,  who  applied  them  to  their  velocity  measurements  on 
the  NO^-NaG^  system,  and  calculated  that  the  rate  constant 
for  the  dissociation  reaction  N2O4  = 2NG2  is  x ].0^ 

A _ 1 

t 0.5  X 10  sec.  at  25°C.  and  260  mm.,  with  an  activation 

energy  for  dissociation  of  13.9  t 0.9  kg. cal ./mole. 

Using  frequencies  in  the  ultrasonic  range  measurements 
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have  also  been  made  by  Olson  and  Teeter,  Teeter, 
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Kistiakowsky  and  Richards,  , and  Knescr  and  Gauler. 

Olson  and  Teeter  concluded  that  the  critical  range  had 
actually^  been  passed  at  a frequency  of  51,570  cycles  per 
second  in  the  gas  at  565  mm.  pressure  and  25°C.  In  later 
investigations  on  this  system  Teeter  concluded  that,  due  to 
the  interference  of  vibrational  heat  capacity  relaxation, 
it  was  not  possible  to  determine  the  rate  constant  in  this 
way.  Kistiakowsky  and  Richards  have  carefully  studied  a 
wider  range  of  conditions  and  have  found  that  to  within 
0.5^  there  is  no  systematic  change  of  velocity  with  fre- 
quency between  10  and  SOkc.  For  the  highest  frequency  the 
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measurements  were  carried  down  to  a pressure  of  300  mm.  of 
Hg.  These  investigators  found  that  the  greatest  soiirce  of 
error  was  evidently  snail,  almost  undetectable  traces  of 
impurities.  The  problem  was  restudied  by  Kneser  and  Gauler 
who  agreed  with  Taeter  on  the  perturbing  effect  of  the  vi- 
brational relaxation,  and  also  noted  the  fact  that  the  ideal 
gas  law  could  not  be  applied  to  a gaseous  system  that  con- 
densed near  room  temperature. 

This  oroblem  of  ideality  had  already  been  considered 
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by  Luck,  who  extended  Einstein's  treatment  to  the  case 
of  mixtures  of  real  absorbing  gases.  He  concluded;  "The 
departure  of  a reacting  gas  mixture  from  ideal  properties 
cannot  change  the  type  of  sound  velocity  dependence  upon 
frequency  in  the  mixture,  but  can  affect  nunerical  values, 
and  so  can  probably  explain  the  discrepancies  between 
experiment  and  the  earlier  theory. 

"The  practical  effect  of  absorption  on  the  precise 
measurement  of  sound  velocity  in  the  dispersive  frequency 
region  is  to  render  such  measurement  impossible.  Precise 
velocity  measurements  in  the  dispersive  region  are  essential 
for  reaction  rate  determination  since  very  small  velocity 
differences  are  involved,  but  the  absorptiori  maxinura  should 
be  sufficient  in  actual  gases  to  reduce  sound  intensity  to 
less  than  2%  of  its  initial  value  in  14  wave-lengths,  making 
the  usual  interferometric  velocity  measurements  both  very 
difficult  and  very  inaccurate. 
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”0n  the  other  hand,  it  may  be  seen  from  the  graph  that 
absorption  measurements  offer  a promising  means  of  deter- 
mining the  reaction  rate.  The  curve  of  absorption  versus 
frequency  shows  a reasonably  sharp  maximum.  ** 

Unfortunately,  absorption  coefficients  as  a function 
of  frequency  are  difficult  to  determine  with  accuracy. 

The  whole  problem  has  been  considerably  clarified  by 

tl  5 9 

a long  paper  by  Damkohler,  who  considered  the  propagation 
of  a sound  wave  in  a gaseous  system  in  which  any  number  of 
reactions  are  taking  place.  In  this  paper  all  other  loss 
factors  such  as  thermal  conductivity,  thermal  relaxation, 
and  viscosity  were  taken  into  account.  He  concluded  that 
for  the  case  of  NOj-NjO^  there  may  well  be  interference  be- 
tween the  relaxation  effects  due  to  reaction  and  those  due 
to  the  vibrational  heat  capacity. 

The  conclusion  with  which  most  of  the  above  investi- 
gators concur  is  that  the  available  sonic  aata  can  at  best 
give  a minimum  value  for  the  velocity  constant  for  the 

A _ 1 

NaO^  dissociation  reaction,  viz. . 5x10  sec. 

b.  Dissociation  ■under  a Temperature  Gradient 
60 

Dirac  suggested  in  1924  that  if  a gas  such  as 
nitrogen  dioxide-tetroxide  was  subjected  to  a temperature 
gradient,  the  concentrations  necessary  for  equilibrium  would 
be  different  at  different  points  so  that  there  would  be 
diffusion  between  the  two  kinds  of  molecules,  a steady  state 
being  attained  when  the  rate  of  diffusioii  of  double  molecules 
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into  any  region  is  equal  to  their  excess  rt.te  of  decoin- 
position  over  their  rate  of  formation  within  that  region.  | 

He  develops  equations  vhich  would  enable  one  to  compute 
the  reaction  rate  In  terms  of  other  physical  constants  of 
the  gas  from  a detailed  knowledge  of  the  concentration  or 
tem])erat\ire  gradient  under  steady  state  for  a fixed  over-  i 

i 

all  temperature  gradient.  However,  for  any  practical  c.ise  | 

the  quantities  to  be  measured  fall  off  sn  rapidly  over  the  I 

region  of  interest  that  this  method  would  be  of  no  value  1 

1-1 

in  determining  rates  of  the  order  of  1x10  sec.  , or  f; ster. 
c.  Method  of  Brass  and  Tolman 

I 

llais  method  has  been  succinctly  described  by  Brass  and  I 
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Tolman  as  follows:  i 

”The  general  nature  of  the  method  was  to  cause  an  equi-  ^ 

librium  mixture  of  the  two  gases  to  flow  with  a sudden  drop  •' 

in  pressure  through  a perforated  diaphragm,  and  take  measure- 
ments of  the  temperature  of  the  gas  both  before  passage 
through  the  diaphragm  and  at  successive  positions  in  the  i 

path  of  flow  beyond.  The  sudden  drop  In  pressure  and  con- 
centration on  passage  through  the  diaphragm  should  be 
followed  by  an  increase  in  the  degree  of  dissociation  of 
the  tetroxide  and,  since  the  reaction  is  endothermnl , thin 
should  be  accompanied  by  a drop  in  the  temperature  of  the 
gas.  Hence,  knowing  the  velocity  of  flow  and  the  tempera- 
ture at  different  distances  beyond  the  diaphragm,  it  was 
hoped  that  conclusions  as  to  the  rate  of  dissociation 
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could  be  drawn.” 

This  method  has  some  advantages  over  the  acoustic 
method  since  the  equipment  required  is  much  less  compli- 
cated and  yields  information  as  to  the  rate  more  directly. 
Unfortunately  there  Is  a great  deal  of  uncertainty  due  to 
lack  of  knowledge  about  the  behavior  of  jets  issuing  from 
the  fine  holes  in  the  diaphragm  and  difficulty  in  making 
precision  teraneratiAre  measurements  on  the  flowing  gas  since 
the  thermocouples  not  only  drain  heat  from  the  system  but 
also  provide  local  regions  where  friction,  compression, 
and  possibly  catalytic  effects  may  disturb  the  postulated 
uniform  gradient  in  reaction. 

The  best  estimate  of  the  uiiimolecular  dissociation 
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rate  constant  based  on  this  experiment  is  o.2xl0  sec.  at 
1 atm.  and  25®C.,  but  this  value  the  authors  .state  may  be 
in  error  ”even  as  to  order  of  magnitude”.  There  is,  how- 
ever, a definite  positive  effect  which  the  authors  feel 
can  only  be  explained  by  a high  rate  of  dissocation  of 
nitrogen  tetroxide. 

d.  Shock  Tube  .-Method 

Using  the  conventional  shock  tube  method  Carrington 
6 2 

and  Davidson  have  passed  a shock  wave  through  an  equi- 
librium mixture  of  NCj  and  Passage  of  the  shock  wave 

provides  a rapid  adiabatic  non-isentropic  compre.ssion,  and 
the  subsequent  dissociation  of  N2O4  is  followed  photo- 
electrically  with  light  of  such  wavelengths  as  to  be  ab- 
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sorbed  preferentially  by  either  NO?  or  NjO*.  Using  this 
method  they  studied  the  dissociation  of  N2O4  in  the  presence 
of  a large  excess  of  nitrogen  over  the  temperature  range 
-20®  to  28°  and  over  the  nitrogen  pressure  range  0.5  to  4 
atm.  Around  1 atm.,  the  rate  law  is  = k(N204)(N2); 

k = 2.0x10  exp  (-11000/RT)  liter  mole-’  sec."’.  At  high- 
er nitrogen  pressures,  the  rates  are  less  than  the  values 
predicted  by  the  above  equation  and  it  is  estimated  that 
the  limiting  first-order  rate  constant  for  the  dissociation 
of  N2O4  is  k = 10’*  exp  (-12900/RT)  sec."’ 

This  method  is  analogous  to  one  limit  of  the  method 
applied  in  this  thesis,  and  will  be  discussed  extensively 
in  this  light  in  sections  II  and  V.  It  has  the  advantages 
of  reouiring  simpler  equipment  and  only  small  amo\mts  of 
gas.  The‘"e  may  be  some  disadvantage  involved  in  not  being 
able  to  vary  the  rate  of  energy  input  as  this  means  that 
displacements  from  equilibrium  will  be  large  for  measurable 
shocks. 


II.  Theoretical  Analysis  of  Impact 
Tube  Experiments 

A.  General  Wsture  of  this  Experiment 

Impact  tube  experiments  are  based  upon  the  sudden 
change  in  energy  distribution  experienced  by  a rapidly 
flowing  fluid  when  it  impinges  upon  a stationary  ob^iect 
of  such  small  dimensions  as  not  to  alter  greatly  the 
general  flow  pattern.  In  this  case  the  object  is  an  impact 
tube  at  the  face  of  vdiich  the  gas  flow  is  stopped  and  the 
energy  of  mass  motion  is  rapidly  converted  into  random 
molecular  motion,  which  is  then  redistributed  among  the 
various  other  degrees  of  freedom.  iCnowledge  of  the  steady 
state  pressure  at  the  face  of  the  impact  tube  can  be  made 
to  yield  information  about  the  nature  of  the  energy  ex- 
change processes  talring  place.  As  additional  information 
one  raust  have  the  pressure,  temperature,  and  velocity  of 
the  gas  before  impact  and  the  size  of  the  impact  tube. 

The  production  of  the  gas  streaiii  by  flow  through  a 
nozzle  and  its  subsequent  compression  are  shown  diagram- 
raaticelly  in  figure  2.  Figure  3 shows  the  velocity  profile 

with  respect  to  distance  and  tine  along  the  central  stream- 

♦ 

line  approaching  a source  shaped  impact  tube.  Initially 


The  streamlines  around  the  impact  tube  may  be  computed 
according  to  the  method  outlined  on  p.  212,  of  ref.  1. 
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the  gas  is  in  equilibrium  at  a pressure,  p^,  and  tempera- 
ture, Tq,  and  has  a negligible  flow  velocity,  u^.  Since 
the  pressure  in  the  expansion  chamber  is  maintained  at  less 
than  Pq,  the  gas  flows  through  the  nozzle,  and  under  certain 
conditions  (design  and  operating  details  are  in  following 
section)  one  obtains  over  a region  around  the  center  of  the 
orifice-a  non- turbulent  stream  of  gas  at  a pressure,  pi, 
temperature,  T,,  and  flow  velocity,  Ui . Upon  approaching 
the  impact  tube  the_  gas  along  these  streamlines  is  slowed 
and  brought  to  stagnation  at  the  face  of  the  impact  tube. 

If  one  imagines  that  the  energy  of  mass  motion  is  converted 
instantaneously  into  random  molecular  motion,  the  pressure 
and  temperature  change  as  a consequence  to  new  instantan- 
eous values  Pa,  T^,  which  relax  eventually  to  equilibrium 
at  P3,  T3 , u=0.  Under  steady  state  conditions  there  is  no 
net  gas  flow  across  the  face  of  the  impact  tube  and  hence 
cne  may  imagine  the  face  of  the  impact  tube  to  be  covered 
by  a thin  diaphragm  balanced  for  equal  pressures  inside  and 
outside.  Thus  the  pressure  measured  by  a manometric  device 
attached  to  the  impact  tube  is  the  stagnation  pressure  at 
the  face  of  the  tube. 

If  the  time  over  which  pressure  changes  occur  in  flow 
through  the  nozzle  is  large  compared  to  the  equilibration 
time  of  the  gas  then  the  adiabatic  expansion  through  the 
nozzle  may  be  considered  to  occur  by  a series  of  infinitesi- 
mal equilibrium  steps  and  the  entropy  change  will  be  zero. 
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The  enthalpy  of  the  gas  will  be  lessened  by  an  ainoiint  equal 
to  the  kinetic  energy  of  mass  motion  of  the  gas. 

On  encountering  the  face  of  the  impact  tube  the  kinetic 
energy  of  mass  motion  is  returned  to  the  gas  as  random 
molecular  motion.  For  an  experiment  of  this  type  the  com- 
pression time  depends  upon  the  ratio  of  the  size  of  the 
impact  tube  face  to  the  streaming  velocity  of  the  gas,  the 
gas  flow  not  being  effected  by  the  presence  of  the  impact 
tube  until  it  is  within  a distance  equal  to  a few  times  the 
breadth  of  the  impact  tube  face  (see  figure  3).  \<ith  gas 
flow  velocities  amounting  to  several  tenths  of  a Mach  number 
and  impact  tubes  having  a circular  cross  section  of  approxi- 
mately 0.03  cm.  diameter  compression  times  of  the  order  of 

6 

1x10"  sec.  can  readily  be  attained. 

If  this  compression  time  is  long  compared  to  the  equi- 
libration time  of  the  gas,  equilibriuin  will  be  .maintained 
during  compression  and  the  latter  will  be  isentropic;  the 
final  temperature  and  pressure  will  clearly  equal  those  of 
the  reservoir.  If,  on  the  other  hand,  this  compression  time 
is  very  short  compared  to  the  equilibration  tine  of  the  gas, 
the  energy  input  will  occur  instantaneously  compared  to  the 
equilibration  time  of  the  processes  involved,  and  the  maxi- 
mum entropy  increase  will  occur.  In  order  to  obtain  infor- 
mation about  the  rate  of  the  equilibration  processes  it  is 
necessary  to  have  the  compression  time  and  equilibration 
time  of  the  same  order  of  magnitude,  for  then  the  fraction 
of  the  Ttiaximiun  entropy  increase  \Aiich  is  produced  may  be  used 
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to  determine  the  equilibration  time  in  terms  of  the  known 
compression  time.  After  a very  rapid  compression  the  tem- 
perature vdiich  corresponds  to  the  translational  energy  of 
the  molecules  rises  to  a peak  value,  Ta,  and  the  pressure 
becomes  Pa.  In  general  this  pressure  will  not  equal  p^,  and 
Griffith  has  shown  that  for  a non-reacting  ideal  gas  Pa  < Pq. 
That  this  is  also  true  for  a reacting  mixture  of  ideal  gases 
will  be  shoTTn  in  part  C of  this  section.  During  the  subse- 
quent adjustment  period  the  temperature  corresponding  to  the 
translational  modes  decays  to  some  equilibrium  value,  T3, 
■vdiile  the  other  modes  (i.e,  vibration,  reaction)  gain  energy 
and  follow  a different  path  _with  energy  distributions 
corresponding  to  a special  ’’temperature"  “Under  the 
conditions  of  this  experiment  this  adjustment  occurs  at  a 
constant  pressure,  Pz,  and  the  temperature  change  is  compen- 
sated for  by  an  increase  in  the  density,  or  by  an  appropri- 
ate shift  in  the  degree  of  dissociation  for  a chemically 
reacting  mixture. 

♦ 

Figure  4 is  a schematic  drawing  showing  these  changes. 
In  this  figure  represents  a characteristic  reaction  time 
of  the  proper  order  of  magnitude  to  make  this  experiment 
yield  significant  results  and  c is  the  average  velocity  of 
random  molecular  motion. 


* 


This  figure  based  in  part  on  figure  1 of  ref.  44b, 
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B.  Theory  for  a Nonreactinc  Ideal  Gas 

In  this  subsection  the  theoretical  equations  (govern- 
ing the  impact  tube  experiment  as  applied  to  a nonreacting, 
ideal  gas  are  developed.  To  a large  extent  this  develop- 

1 44 

ment  is  based  on  the  work  of  Kantrowitz  and  Griffith. 

The  general  case  of  a real  non- dissociating  gas  has  been 
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discussed  in  some  detail  by  Bauer.  The  important  equa- 
tions (7),  (15),  and  (19)  are,  however,  developed  by  an 
independent  method,  which  the  author  feels  is  at  once  clear- 
er and  less  intuitive  than  earlier  approaches.  An  inaccu- 
racy in  the  earlier  work  is  also  pointed  out.  The  guiding 
principle  in  formulating  this  section  has  been  to  develop 
the  equations  in  such  a fashion  that  their  extension  to  the 
case  of  a reacting  mixture  of  gases  in  the  following  section 
could  be  easily  followed. 

In  order  to  specify  completely  the  state  of  a typical 
non-dissociating  gas  like  carbon  dioxide,  the  pressure  and 
a single  temperature  suffice  provided  the  system  is  in  equi- 
librium. The  nature  of  this  experiment  is  such  that  at  the 
face  of  the  impact  tube  the  translational  energy  is  increased 
so  rapidly  that  rhe  system  does  .not  remain  in  thermodynamic 
equilibrium.  The  rotational  modes  follow  this  change  in 
translational  energy  within  a comparatively  short  time,  and 
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on  the  time  scale  of  this  experiment  one  nay  assume  that 
they  are  in  equilibrium  with  the  translational  modes.  The 
vibrational  modes  are  not  directly  coupled  with  che  trans- 


lational  inodes  and,  although  in  general  there  is  coupling 
between  the  vibrational  and  rotational  modes,  it  is  usually 
very  weak  so  that  >Aiereas  the  latter  follow  rapid  changes 
in  translational  energy,  the  former  do  so  only  after  an 
appreeiable  (of  the  order  of  microseconds)  time  lag.  How- 
ever, the  vibrational  modes  of  a molecule  are  rather  strong- 
ly coupled  to  each  other  and  maintain  an  equilibrium  energy 
distribution  even  under  these  conditions.  Accordingly,  one 
may  define  two  temperatures:  T as  the  temperature  in  the 

Boltzmann  formula  corresponding  to  the  population  distri- 
bution among  translational  and  rotational  states,  and 
as  the  temperature  in  the  Boltzmann  formula  corresponding  to 
the  population  distribution  among  vibrational  states.  The 
specific  heat  may  then  be  regarded  as  arising  from  two 
sources,  vibrational  motion,  and  all  the  rest.  The  vibra- 
tional contribution  is  where  £yib 

energy  of  vibrational  motion.  Cp  refers  to  the  specific  heat 
at  constant  pressure  not  including  vibrational  degrees  of 
freedom.  For  = T one  has  the  usual  total  s))ecific 

I 

heat,  ^vib*  ideal  gas  equation  of  state 

will  be  assumed  throughout. 

The  entropy  change  involved  in  each  step  of  the  impact 
tube  method  for  the  idealized  case  of  instantaneous  com- 
pression will  now  be  disc\issed.  This  will  clearly  give  the 
maximum  possible  entropy  gain.  The  enthalpy  changes  will 
also  be  followed. 
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The  expansion  through  the  nozzle  is  sufficiently  slow 
compared  to  the  rate  of  equilibration  so  that  the  process 
proceeds  reversibly.  Ihe  expansion  is  also  adiabatic 
(C;  = 0)  . 

2 

ASqi  = 0 and  AHqi  = -hu,/2  . (.1) 

Here  M is  the  moieciilar  weight  of  the  gas  and  u,  is  the 
unit  mass  velocity  of  the  stream  in  state  1. 

For  a truly  instantaneous  compression  the  entropy 
change  is  equal  to  zero  since  the  non-lagging  modes  follow 
reversibly  and  the  lagging  modes  are  in  nowise  effected. 

The  energy  of  mass  motion  is  returned  to  the  gas*  non- 
lagging  inodes. 

2 

AS  ^ 2 “0  and  AH » 2 ~ Mu  ^ / 2 . ( 2 ) 

During  the  equilibration  period  following  such  a com- 
pression, the  experiment  is  so  arranged  (see  figure  2)  that 
the  pressure  remains  at  P2  vdiile  the  non-lagging  degrees 
of  freedom  go  from  T2  to  Tj  and  the  lagging  degrees  of 
freedom  go  from  T,  to  T3 . The  entropy  increases,  but  there 
is  no  change  in  the  enthalpy. 


A52  3 


and  AH  2 3 


C.  (3) 


Now  consider  going  directly  from  che  initial  state 
(0)  to  the  j’lnal  state  (3).  The  pressure  changes  from 
to  P2,  while  the  temperature  of  all  modes  goes  from  Tq  to 
T3;  hence 


^03 


03 


i 


and  AH 


= 0.  (4) 
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•Since  any  reversible  nnth  may  be  used  to  evaluate  the 


entropy  change, 


AS03  “ ^01  12  ^^2  3, 


r>nd  substituting 


^rom  eouations  (l),  (2),  (3),  and  (4)  one  gets 


Rln  = C, 


I 1 

In  + C In  -O 
T,  P Ta 


Note  thrt  T3  cancels  out  since  C = C„  + C ..  . Since  all 

p P vib 

of  these  steps  occur  adiabatic-a?  ly  the  tot'.  1 energy  of  the 
system  is  conserved  and  hence  considering  states  0 and  2 


Cp(Ta-T„)  . C^ib(T.-To)  = 0. 


vSo].ving  (6)  for  T2  and  substituting  into  (^)  gives 


£0  = ) ^vib 

P2  ” I C., 


ri.b  Tq 


c;/R 


n 

fp 


rib/^ 


From  A?Qt  = 0 one  nay  deduce  the  familiar  adiabatic 
relation  for  an  ideal  gas: 

Rln  = Cp  in  ^ . (S) 

Combining  equations  (6),  (7),  and  (3)  then  allows  one  to 
compute  the  maximum  'instantaneous  compression)  pressure 
defect  (pQ-Pa).,  . Equation  (6)  is  needed  to  give  T2 . 

For  any  nonreacting  gas  with  lagging  and  non-lagging 
degrees  of  freedom  the  tine  rate  of  entropy  change  may  be 


written  as 

IS  _ 


dt  V)p^T,T^.vdt 


-3  4-  V \rr  ' 


P^'^vjb  ^<^l’vib^  o.T 
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The  coefficients  to  be  used  here  are: 


I < 


(iS)  , 

dp  T,T^ib  P.T^ib 

(^)  = f c' 

^P.Tvib  ^ ’’ 


/ S \ _ 1 p 

'^vibp.T  ■ "fvlb 


The  hydrodynaniicpl  momentum  conservation  equation  yields 
dvi  ^ 

^ ~ 2 dt  ^ 'diich  may  be  rewritten  for  an  Ideal  gas  as 

^ - -P-  r M-l  i r n ' 

dt  ■ RT  2 dt  ^ ■ 

Substitution  of  (10)  and  (11)  into  (9)  yields 

d^  1 M dui  . £e  ^ + ^vib  ^^'vib  y 

dt  ■ f 2 dt  T dt  Tviij  dt  * ^ 

Treating  the  time  rate  of  enthalpy  change  in  a simi- 
lar manner : 

M = r— ) ^ dT 

dt  dt 

* (3^2  T . (IS 

<>Tvib  p,T  dt 


2 dt 


For  an  ideal  gas  (r^H/dp)fp  m = C,  and  substituting  the 

appropriate  values  for  the  other  coefficients  gives 

^ = p'  ^ + r 1^’vlb  _ M du^  /, 

dt  ^p  dt  ^vib  dt  " ■ 2 


Using  (I4)  to  eliminate  the  derivatives  of  u^  and  T from 
(12)  then  yields  the  important  relationship 
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r 


^ = r I 1 1 1 ^^vib 

dt  ^vib  T J dt 


(15) 


One  ..lay  now  find  cuite  simply  the  entropy  gain  lor 
instantaneous  compression.  After  such  a compression 
du^/dt  = 0 and  equation  (i4)  yields 


t 


dt 


+ C 


<^^vib  _ 


vib  dt 


= 0 


(16) 


This  may  be  integrated  using  as  the  boundary  condition 
" ^ = T-3  to  give 


T 


'vib  ni 


Vib 


(17) 


Combining  (17)  with  (15)  and  integrating  over  the  ecuili- 
bration  period  following  the  instantaneous  compression  one 
gets 


^vib  ’ 


vib  ,p 

^vib 


C_ 

T'3 


^"^vib 


(13) 


is  taken  as  the  average  of  (TT^^^^)  over  the  v/hole  pro- 
cess, and  from  figure  4 u good  value  is  seen  to  be 
T = ^(^o  specific  heats  and  equilibration  ti'aes 

will  be  taken  at  this  temperature  as  well.  Sijice  the  pro- 
cess is  adiabatic 


- Toj  * ^ = C.(T,  - 


To) 


and 


(IS)  becomes 


AS 


i.c. 


1 ^vlb 

2 


(M 

\o 


U1  / 


(19) 


This  formula  is  identical  to  that  found  by  Karitrowitz 
and  Griffith  following  a kinetic  formulation.  The  kinetic 
analysis  applies  to  the  case  of  non- instantaneous  com- 
pression as  well,  and  will  now  be  discussed. 

As  was  noted  in  Section  I.,  the  details  of  the  col- 
lision process  must  be  exaiained  in  order  to  find  an  expres- 
sion for  the  rate  at  'tdiich  energy  exchange  will  occur  be- 
tween the  various  modes.  An  assumption  is  introduced; 
namely,  that  after  a small  perturbation  from  equilibriuifl 
the  time  rate  of  adjustment  of  the  vibrational  energy  is 

59,65E 

proportional  to  its  displacement  from  equilibrium. 


(20) 


10^13,66,67 

Semi-classical  treatments  of  this  problem  have 

shown  "t*  to  have  the  following  functional  dependence: 


p'’(TiS) 


(f)-^'^%xp  {(Vs)'/'(|V^^j  . (21) 


Here  V is  the  vibrational  frequency,  s the  molecular 
radius,  and  m the  molecular  mas.^.  When  changes  in  the  tem- 
perature are  small  T may  be  considered  constant,  and  after 
a time  T the  temperature  difference  dirainishes  to  1/e,  or 
0.363,  of  its  original  value.  For  a first  order  process 
this  time  is  clearly  independent  of  tht'  magnitude  of  the 
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startin[;  temperature  di.f f eronce,  and  is  frecuentiy  referred 
to  as  the  "relaxation”  time. 

Kr.ntrowitz,  and  Griffith  as  well,  used  a somewhat 
different  assumption  instead  of  (20),  writing 


^ ^vib  '^vih 
dt 


Cvib 


(?■<') 


That  is,  they  postulated  that  the  time  rate  of  ad.justrierit 
of  vibrationaT  energy  is  proportional  to  the  difference  be- 
tween the  actual  vibrational  energy  and  the  vibratioiial 
energy  for  equilibrium  partition  at  the  translational  tem- 
perature T.  The  difference  between  these  two  assumptions, 
(20)  and  (22),  clearly  depends  on  the  ratio  of  (13-!'^^^) 
to  (T-Ty^^) . Since  the  relaxation  process  is  adiabatic 
Cp(T-T3)  + Cy^^j(Tyi|3-T3 ) = 0,  and  it  follows  that 

_J? — ^ . hence  relexation  times  calculated  according 
T - Tvib  Cp 

to  the  two  different  assumptions  will  be  related  as  follows: 


V/r^  = Cp/Cp  . (a) 

II  5 9 6 5a 

The  work  of  Damkohler,  Manes,  and  others  indicates 
clearly  that  the  assumption  (20)  is  to  be  preferred.  How- 
ever, since  the  various  formulas  and  earlier  treatiaents 
have  been  developed  according  to  assumption  (22),  we  shv^ll 
for  purposes  of  comparison  continue  with  the  latter  except 
\diere  otherwise  stated. 

Following  Kantrevritz,  the  variable  6 , vdiich  represents 
the  excess  energy  in  the  lagging  modes  over  the  energy  at 
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equilibrium  partition  at  the  temperature  T,  is  now  intro- 


duced as  a measure  of  the  disequilibrium. 


(24) 


In  experiments  vrtierein  changes  in  T are  small  compared  to 
T itself,  (15)  may  be  rewritten  using  (Z2)  and  (24)  as 

jr.  _ 6dt  ^"'^vib -V  £ ^ dt 

-•  - ^ • (25) 

^^vib  'A  ^vib  ^ 


is  here  define<^  as  in  equation  (13).  Then 

...  r_li_dt^  . 

" "k  Cvih 


(26) 


Since  ^(t)  deperic's  on  the  compression  profile  and  charac- 
teristics of  the  gas,  the  integral  ;iust  be  taken  over  that 
part  of  the  flow  for  which  is  desired.  By  combining 
equations  (14),  (22),  and  (24)  so  as  to  elliainate  the  tem- 
perature,  r.n  equation  for  6 is  obtained. 


dt  + £p  _L  _ Sib  M ^ 

dt  ■ cir  2 dt 

p K p 


(27) 


The  integral  is 

C(t,)=exp(-J  I exp(/ 

'k  Cp  % Cp 

C_yb 

T J A.  U O • 


(23) 


Knowledge  of  the  velocity  history  of  the  flow  and  the 
specific  heats  fcr  the  gas  now  enable  one  to  relate  pres- 
sure differences  to  relaxation  times. 
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This  approach  may  now  be  utilized  as  an  alternate 
method  for  deriving  (19).  Equation  (27)  may  be  rewritten 


as 


dt 


''P  ''P 

As  the  time  required  for  stopping  becomes  shorter  the 
second  term  will  become  small  comnared  to  the  other  terms, 
and  in  the  limit,  vrtien  becomes  negligible  in  compari- 

son  to  the  other  two  terms,  the  change  in  6 is 

P 

In  the  case  being  considered  6 was  originally  zero  and  the 
velocity  was  u, . A sudden  stopping  then  produces  a value  of 

4 = r?yib\  M 2 

^o  ^ ' 2 

V/ith  u = 0,  (29)  gives  for  the  approach  to  equilibrium 

The  entropy  change  for  instantaneous  compression  is,  accord- 
ingly. 


(^)  at  . (^)  = (?^)(LZ^)  at. 

tv  C dt 


(29) 


ASq3  ) 


- /'^  6^  dt  _ f 

J ^vib  J 


° ^ 
^p  ^vib 


2 Cp  Cvib  T 


( 


'vib 


Mu, 


2 C, 


(30) 


which  is,  of  course,  identical  with  equation  (i9)  since  for 
instantaneous  compression  cancels  out. 
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( ' 


oince  the  expansion  axid  compression  take  place  adia- 
batically  AH03  = 0,  which  for  a nonreacting  ideal  gas 
means  that  T3  = Tq.  The  overall  entropy  change  is  there- 
fore given  by 

As  = R In  P0/P3  . (j;l) 

’■/hen  the  difference  (PQ-P3)  = Ap  is  small  compared  to 


and  from  (19)  or  (30) 

^Pj . c . _ /^vib  \ 

r\  “ _ _ _ ' _ ' * 


^Cp  RT  2 C. 


(33) 


It  is  convenient  at  this  point  to  substitute  a set  of 
dimensionless  variables.  Again  using  Kantrowitz’  notation 


define  the  following: 
‘'K-Cp  d 


t _ u,  t 
d 


(34) 


£'  = 


(W<h|  u,P 


u = u/u. 


Using  (30),  (32),  and  (34)  in  equations  (26)  and  (23)  one 
obtains 

t«  t'  t' 

e’(t’)  = exp(-  j"  dt')  I exp(;^;  dt')f^|"dt« 


and. 


.00 


^i.c.  ' 


— lit'  . 


(35) 


(36) 


j 
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The  scuare  bracketed  portion  of  (35)  arises  from  the  fact 
that  the  energy  inout  on  compressj on  occurs  over  a finite 
tiiue  interval  rather  than  inotantaneouoly . Equation  (35) 
thus  represents  a folding  of  the  energy  input  curve  into 
the  relaxation  expression.  From  (36)  one  sees  that 
AF/APs  r,  depends  only  on  ‘t'*-  for  a given  flow  pattern.  This 
is  the  great  simplification  obtained  by  running  an  experi- 
ment in  which  only  small  temperature  changes  occur. 

In  order  to  actually  calculate  relaxation  times  from 

observed  pressure  defects  it  is  necessary  to  know  the  value 

of  the  integral  appearing  in  (36),  hence  of  t’(t').  The 

latter  requires  a detailed  knowledge  of  the  coapression 

streamlines  in  the  region  around  the  face  of  the  impact  tube 

£i.e.  of  (t;3  which,  in  turn,  depend  on  the  geometry  of 

the  impact  tube.  Tne  flow  around  a source- shaped  impact 

tube  can  be  described  analytically,  and  Kantrovrits  has 

evaluated  the  above  integral  for  this  special  case.  For  the 

square- tipped  imyiact  tribes  used  in  our  experiments  no 

analytic  solution  is  available,  but  Griffith  has  carried 

out  the  necessary  computations  numerically  using  Southwell's 

68 

Relaxation  Method  to  obtain  the  flow  pattern.  This  I'-iethod 
rests  upon  the  accuracy  with  which  a aifference  equation 
can  be  made  to  approximate  a differential  equation.  Assuming 
incompressible  flow  impinging  axially  upon  a square-tipped 
impact  tube  having  a hollow  center  of  inside  diameter  equal 
to  0.620  times  the  outside  diameter,  he  obtains  the  values 
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shown  in  figure  5 for  the  dependence  of  velocity  on  distance 
and  the  values  shown  in  figure  6 for  the  dependence  of  € ' 
on  the  ticie.  Upon  integrating  over  the  path  for  various 
't ' he  obtains  the  fcllowing  values: 


'^'k' 

1/10 

1/3  1 3 

10 

0.0328 

0.226  0.480  0.723 

0.903 

For  very  short  relaxation  times  the  approxi.nation  formula 
Ap/aPt  - = 0.867  V may  be  used.  This  is  valid  for 
T;  » ^ O.OE.  This  completes  the  theory  necessary  for  re- 
ducing observed  pressure  defects  to  relaxation  times. 


FIGURE  5 


VELOCITY  ALONG  AXIS  FOR  FLOW  AROUND  A PITOT  tube 


C.  Theory  for  a Reacting  Mixture  of  Ideal  Gases 


In  this  subsection  the  theoretical  equations  governing 
the  impact  tube  experiment  as  applied  to  a reacting  mixture 
of  ideal  gases  are  developed.  This  constitutes  an  exten- 
sion of  the  application  of  this  type  of  experiment,  and 
was  first  discussed  by  Bauer  . For  the  purposes  of  the 
present  analysis  consideration  will  be  limited  to  a gas 
wiiich  dissociates  according  to  the  equation  Az  = <.A,  but 
the  formulas  obtained  may  be  readily  extended  to  other  cases 
of  pressure  dependent  chemical  reactions.  It  will  be  as- 
sumed that  the  species  Aj  and  A individually  obey  the  ideal 
gas  law. 

In  order  to  specify  completely  the  state  of  a typical 
dissociating  gas  like  nitrogen  tetroxide,  the  pressure  and 
a single  temperature  suffice  provided  the  system  is  in 
equilibrium.  As  has  been  noted  in  part  A of  this  section, 
the  nature  of  the  impact  tube  experiment  is  such  that  at 
the  face  of  the  impact  tube  the  translational  energy  of 
the  gas  is  increased  so  rapidly  that  a condition  of  dis- 
equilibrium results.  One  then  needs  two  additional  vari- 
ables to  specify  completely  the  state  of  the  gas:  n,  the 
degree  of  dissociation  and  0,  the  temperature  in  the 
Boltzmann  formula  corresponding  to  the  pupulation  distri- 
bution among  vibrational  states.  This  latter  definition 
parallels  that  used  for  a non-dissociating  gas  and,  as 
before,  T is  talcen  as  the  temperature  in  the  Boltzmann 
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formula  corresponding  to  the  population  distribution  among 
translational  and  rotational  states.  (Note:  At  first 

though  one  might  be  inclined  to  postulate  different  0’s 
for  the  lagging  modes  of  NO2  and  Rigorously  one 

should  assLime  a different  0 for  each  individual  lagging 
degree  of  freedom.  Since  this  exact  assumption  would  lead 
to  almost  intractable  mathematical  expressions  (see  Bauer, 
ref.  63)  we  make  the  simplest  assumption  of  a single  9 for 
all  of  the  lagging  vibrational  modes.  A senarato  variable 
is  used  to  describe  the  reaction  which  may  be  expected  to 
adjust  at  an  appreciably  different  rate  from  vibration.) 

a.  Thermodynamic  Considerations. 

The  entropy  change  involved  in  each  step  of  the  im- 
pact tube  method  for  the  ideclii^ed  case  of  instantaneous 
compression  will  now  be  discussed.  This  will  clearly  give 
the  maximum  possible  entropy  gain.  The  enthc^py  changes 
will  also  be  followed.  Since  ^ and  AS  of  reaction  are 
given  for  changes  at  25°C,  and  pressures  for  each  component 
of  one  atmosphere  (AH°98  and  AS%s)>  ^^ne  must  include  suit- 
able correction  terras  when  the  changes  occur  under  con- 
ditions other  than  these.  All  quantities  will  be  written 
per  original  mole  of  A-  unless  otherwise  specified. 

The  expansion  through  the  nozzle  is  sufficiently  slow 
compered  to  the  rate  of  equilibration  that  the  process  pro- 
ceeds reversibly.  The  expansion  is  also  adiabatic. 

ASo,  = 0;  AHo<  = - I u,'  . 


(37) 
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Here  M is  the  molecular  weight  of  Aj,  and  Ui  is  the  unit 
mass  velocity  of  the  stresm  in  state  1. 

For  a truly  instantaneous  compression  the  entropy 
change  (/5s3iz)  is  equal  to  zero  since  the  non-lagging  modes 
follow  reversibly  >diile  the  lagging  modes  are  in  nowise  - 
affected.  The  energy  of  mass  motion  is  returned  as  molecu- 
lar enthalpy. 

^ * (3o) 

During  the  equilibration  period  afterjsuch  a compres- 
sion the  experiment  is  so  arranged  (see  figure  2)  that  the 
pressure  remains  at  p 2 while  the  non-lagging  specific  hest 
goes  from  T 2 to  theiraggih^  heat'bf-the'vi- 

-brational  modes  goes  from  T,  _to  Tj,  and  the  lagging  degree 
of  dissociation  g^es  from  ai_to  03,  The  entropy  increases, 
but  there  is' no  change  in- the  enthalpy. 

Vl)f 

^0(2)  ^ ^d(l)J' 

+ (l-ai)R  ln(^;  P2)  + 2ai  HlnC^i^  p^) 

+ ^°(a3-a,)-(l-o'3)Rln(±^^  Pa) 

1+G3 

.2a3Rl„(_|a  p,) 

+ {(1-c.O  + 2c,  j in  Tj/T,, 

^ = 0.  (39) 

X • V..  • 
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Here  C refers  to  the  total  specific  heat,  refers 

to  the  non-lagr;irtg  specific  heat  (translation  and  rotation), 
and  refers  to  the  la^itjing  specific  heat  (vibroT.ion)  . 

The  jubscriut  1=1  refers  to  species  A,  and  i = P refers 
to  species  A2.  The  sneclfic  heats  are  taken  per  for-nala 
weight  of  the  substance  indicated.  This  expression  fer 
^23);  h'^en  Goriruted  using  a path  such  that  the 

O 

change  in  n occurs  under  conditions  wliere^^Sjoe  I'^rcperiy 


applies. 

Since 

this  ar 

gu’aent  i 

3 used 

reper  tedly 

the  St. 

for  this 

case  are  presented  in 

detail 

below. 

State  T 

State 

> 

u: 

a, 

Gl 

f'3 

'J? 

*23 

T: 

Tj 

To 

To 

To 

‘ 0 • 

T3 

6: 

T, 

To 

m 

■^0 

To 

'p 

' 0 

3 

PNO2 : 

2a, 

l+a,P2 

1 

2^3  ^ 

1+Ui 

X 

X 

l+a,P2 

PNjO*: 

1-a, 

Lp2 

1+c:, 

1-a,  _ 

Up  2 

1+a, 

1 

1 

1-c,  . 

l + i'l3 

1-n, 

; — 
1+03 

It  should  be  noted  at  this  point  that  the  specific  hert 
terns  of  ecuaticjn  (39)  demarkod  by  curly  brackets  are 
strictly  spoaliing  tenperature  dependent  and  represent  aver- 
age value:-  ove'^  a teii:pers tiire  interval  with  1/T  as  the 
weighting  junction.  Curly  brackets  w:il  be  used  throughout 
this  section  to  signify  such  special  averages. 

Now  consider  going  directly  fror.  the  initial  state  to 
the  final  state.  The  changes  are: 
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P 

T = 0 
a 


Po 

T 

•^o 


P2 

T3 

a 


3 > 


and  hence: 


^0>) 


= (l-ao)Rln(^  Po)  . SooRlnCf^  Po> 

-^iS°(ao-a3)  - (l-a,)  Rln  pj 

-2a.3Rln(|^~  p,)  + {(l-a3)Cp(2) 

+ 2a3  J In  Tj/Tq 


(40) 


and  AH03  “ 2 ) ^ 'vT3~Tq)  + AH°(c3~c:o^  “ ^ 

Since  any  reversible  path  may  be  used  to  evaluate  AS, 


ASo3)j^^^^  = ^ ^^23) 

and  upon  substituting  from  equations  (37),  (33), 
(40)  one  gets  after  combining  terms: 
(l+a,)Rln(po/p2  )j^  = |(l-a,  )cj^.2)  + 


i.o.  ^ 

(39),  and 


j In  T^/T2 


+ I )C0^2)  ^®’^0(l)  J ^o/^’ 

+ AS°(ao-cx,)-(l-a,)Rln(~^  . \^) 

^ l+^o 

_ 2a,Rln(|^  . i^)-(a„-a,  )Rin  Kp(To>  . (4I) 

Note  that  T3  has  canceled  upon  combining  terns.  Kp(T)  is 
the  pressure  independent  equilibrium  constant  evaluated  at 


T. 
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KpC-r)  = p . (42) 

Here  Is  the  ecuilibrium  degree  of  dissociation  at  p,  T. 

Since  all  of  these  steps  occur  adiabaticaliy,  one  may 
write  ^02 )s  « - Oi  hence 

[ {(I'O,  )Cg(2)  + 2a,Cg(]^j|  T,  - |(l-ao)C9(2)  + 2’o''0(l)}^o] 

^ [{a-.)S(2)  ^ 

= 0 . (43) 

Writing  the  explicit  equation  for  ASqi  = u yields  the 
adiabatic  relation: 

♦ f 

(1+a,  )Rln(Po/pi ) = Cp^)  | In  T^/T, 

+ AS°(ao-a,)  - (1-0,)  B In  (“^  • 

- 2o,  R In  (i^  . iii!.)  . (oo-a,)  R Ln  (T,)  . (44) 

Combination  of  equations  (4I),  (4^),  (43),  and  (44) 
allows  one  to  co.mpute  the  maximum  (instantaneous  compres- 
sion) pressure  defect,  (Po-Pi)^  „ > for  a given  Po,  Tq,  p, . 
Equation  (43)  is  needed  to  give  Tj. 

One  sees  then,  that  the  thermodynamic  approach  can 
only  give  limiting  values:  when  the  compression  is  very 

slow  compared  to  the  equilibration  rate  AS03  = 0 and 
P2  = Po>  ^he  compression  is  very  rapid  compared  to  the 


In  ec.  (44)  and  thereafter  AS^a  will  be  abbreviated  to 
AS°  and  AH298  vill  be  abbreviated  to  AH®* 
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( 


\ '• 


equilibration  rate  r 0 and  (p^-Pz)  > ^-nd  is  ^^iven 

by  (41)*  In  order  to  find  the  entropy  increment  for  inter- 
mediate cases  one  must  consider  the  details  of  the  kinetic 
process  and  find  the  entropy  increment  as  a function  of 
the  time  rate  of  enthalpy  input  and  of  the  rates  of  equi- 
libration of  the  lagging  modes. 


b.  Time  Dependent  Entropy  Changes 

For  the  general  case  of  a mixture  of  two  reacting 
gases  related  by  Aj  = 2A  with  both  lagging  and  non-lagging 
degrees  of  freedom,  the  time  rate  of  entropy  change  may  be 
written  as 


dt 


aP  T,0,a  at  ^ 


(iS) 


p,T,a 


dO 

dt 


+ 


/'c)S\  da 

p,T,0  dt  • 


(45) 


The  coefiicients  to  be  used  here  are: 


/c)3\  , _ _ — - C l'*'a)R 

T,9,a  ^ P 


(individuailv 
ideal  gases; 


(46) 


In  Tq/T 
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All  of  these  coefficients  are  per  initial  mole  of  H2O4. 

The  last  coefficient  was  obtained  by  the  method  indicated 
in  the  discussion  of  equation  (39),  and  the  specific  heats 
are  the  same  type  of  weit;hted  average  values  noted  previ- 
ously. The  hydrodynamical  momentum  conservation  equation 
yields 


dx  _ £ du^ 

dt  ‘ ' 2 dt  ’ 

written  for  this  case  as 


dp  _ D 
dt  ■ (i+a)RT 


( 

2 dt'' 


which  may  be  re- 

U'7) 


T is  used  here  rather  than  9 because  only  the  translational 
energy  is  involved  in  the  pressure.  Substitution  of  (46) 
and  (47)  into  (45)  yields 


dt 


1 M ^ 

‘ T 2 dt 


" " 2aCg,^,(e)]a| 

^ (f?ST  P)  * {=;(2)-c;{l)j  1"  VT 


In  T^j/e  + dS° 


] 


da 

dt 


(43) 


Treating  the  time  rate  of  enthalpy  change  in  a similar 
manner; 


dH  _ ^ 

dt  ■ ^c)P  T,e,a  at 


p,0,a 


p,T,a 


dt 


+ (^)  ^ M du^ 

p,T,9  'it  ” 2 dt  * 


(49) 
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For  a mixture  of  ideal  gases  and  with  a constant,  (^) 

T,9,o 

= 0,  and  substituting  the  appropriate  values  for  the  other 
coefficients  gives 

" [=9(2)  - "=9(D j (V6)j  5? 

= - ^ if  • (50) 

Using  (50)  to  eliminate  the  derivatives  of  u^  and  T fro;d 
(43),  yields  the  important  relaaion 

dt  ” f ^■“*=9(2)^®^  * ■‘“=e(i)^®^|  ST 

- i (T„-T) 

* {=9(2) -2=9(1)]  (^d-e)]  i? 

" {=9(2)-2=9(d]  V9  - ^So]f2.  (5i) 

This  equation  corresponds  to  equation  (15)  for  a nonreact- 
ing gas.  Comparison  of  the  two  indicates  the  added  com- 
plexity \rtiich  arises  from  the  need  to  specify  an  additional 
variable,  in  this  case,  case  a,  as  a measure  of  the  chemi- 


cal disequilibrium.  The  parallel  with  this  earlier  work 
may  be  continued  if  one  sets  du^/dt  = 0 in  (50)  and  inte- 
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grates  the  resulting  differential  equation,  using  as  a 
boundary  condition  T = 0 = T3,  a = 03 . This  can  be  readily 
performed  in  this  case  only  if  the  specific  heats  are 
ass’jined  constant,  and  the  coefficient  of  ^ is  assumed  to 
have  a constant  value  of  to  a first  approximation;  then 
one  gets 

{(l-a)c;(j)(T)  + 2aC;(i)(T)}  T + {(l-a)Cg(2) O) 

+ ZaCg^)(e)^  9 + AH°a  = )Cj,(  j,  (Tj  ) 

2a3Cp^j(?3)^  Tj  + |^(l-a3)CQ^^j(T3) 

+ 203  (T,)|  T3  + AHo  03  . 

Now  by  analogy,  one  uses  this  equation  to  substitute  for  T 
in  (51)  giving 

f = <(l-a)Cg(2)(e)  ^ 2aCg^){9)]  ||  - §| 

- fl?7^  [*  ^ l‘'k2)-K(i)l  <VT(e,u)) 

" {Cg(2)-2C9(,)|  (V9)]  ff  . [-H  in  P) 

+ j ^ , where  T(0,a)  is  given  by 


T(9  a)  = - jg-^)C9(2)(e)  ^ 2°Cq(i)(9)^ 
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^ 

[(l-a)C|,(.)(T)  . 2=c;^(,)(T)} 


a 


^ f(l-ai)C^(;)(T,)^2a,c;(^)(I,)| 

{(l-a)C;(2)W  + 2ac'(^)  (T)  j 

|(l-c<3  )C0^2)  (*^3  ) ^^^3^9(1)  (^3 

{(l-a)c;(2)(T)  + 2aCp^)  (T)i 


Unfortunately  this  substitution  leaves  0,  a aivl  their  dif- 
ferent iris  inconveniently  scrambled  so  that  the  next  step 
in  the  earlier  procedure  of  integr?  ting  between  limits  is 
not  feasible  until  0 and  a have  been  expressed  as  functions 
of  the  time.  To  obtain  9(t)  and  a(t)  one  must  consider 
the  detailed  kinetics  of  the  process.  Then  substitution  of 
(51)  into  ^00 


ASq3  ~ AS 


.3=^ 


(gf)  'it  - 


(52) 


and  integration  over  the  compression  path  gives  the  entropy 
gain  in  terms  of  the  kinetic  parameters.  Or,  reversing  this 
procedure,  experimentally  determined  values  of  the  entropy 
gain  could  be  used  to  deduce  values  of  the  kinetic  parameters. 


c.  Consideration  of  tbe  Kinetic  Details 


In  order  to  carry  out  the  integration  indicated  in 
(52)  one  needs  the  tine  dependence  of  p,  T,  0,  and  a.  To 
determine  these  time  dependent  functions  one  needs  four 
independent  differential  equations  of  vhich  tvo  are  pro- 
vided by  the  hydrodynamic  condition  (47)  ana  the  enthalpy 
condition  (50).  An  equation  identical  to  the  vibrational 
relaxation  equation  (20)  of  the  previous  section  mi,^ht  be 
assumed  to  hold  here  also.  A fourth  eauation  may  be  derived 
from  a consideration  of  the  kinetics  of  the  reaction.  If 
one  writes  simply 

Aa  y 2A  , then 

^d 

(^)  = 'dn\ 

reaction  'dt'reaction 


1+a 


^ (at)r  = f [Va,  - 


P 

= k. 


eq 


K 4+a-'  RT  4+a  ^ 

L eq  J 


The  boundary  conditions  for  solution  are  clearly: 

at  tQ:  a = a,  T = T,  0 = T,  P = Pi 

t— ►oo:  0 = 03  T = T3  0 = 13  p = P2. 

Conceptually  this  represents  a comolete  solution  to  the 
problem  since  the  four  independent  differential  equations 

63 

plus  boundary  conditions  specify  a solution.  Bauer  has 
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outlined  an  iterative  procedure  by  which  in  principle  one 
may  compute  p,  T,  Q and  a as  functions  of  the  time  in  terms 
of  the  parameters  Uq,  T^,  p, , the  two  rate  constants,  and 
the  heat  of  activation  for  dissociation.  Evaluation  of 
(52)  for  ranges  of  values  in  these  para..:eters,  and  compari- 
son of  the  calculated  results  with  experimentally  deter- 
mined values  would  enable  one  to  determine  the  parameters. 

While  in  principle  such  an  iterative  procedure  may  be 
utilized,  in  actual  practice,  this  method  would  be  difficult 
to  carry  out  and  there  would  also  be  the  question  of  find- 
ing the  most  rapidly  converging  procedure. 


d.  Introduction  of  Relaxation  Times 

An  alternate  method  of  approach  follows  from  the  ex- 
tension of  the  concept  of  a ^relaxation  time"  to  chemical 
equilibration.  One  then  considers  only  the  over-all  rate 
and  glosses  over  details  of  the  mechanism.  It  is  possible 
to  apply  this  simplifying  concept  because  one  is  dealing 
with  a system  which  is  displaced  only  slightly  from  equi- 
libr lum. 

A number  of  investigators  have  discussed  the  general 
problem  of  the  return  of  a system  to  equilibrium.  Manes, 
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Kofer,  and  V/eller  have  shown  that  in  general  for  re- 
versible reactions  close  to  equilibriuoi  the  net  over-all 
reaction  rate  is  proportional  to  tJi,  where  X is  any  one  of 
the  thermodynamic  functions.  That  the  approach  to  equi- 
libriura  is  first  order,  irrespective  of  the  type  of  rate 

laws  which  govern  the  opposing  reactions,  has  previously 

II  39 

been  pointed  out  by  Daaikohler.  Gilkerson,  Jones,  and 
Gallop  have  shown  that  the  conclusions  of  Manes  £t  al. 
may  also  be  arrived  at  as  a consequence  of  the  theory  of 
absolute  reaction  rates. 

Bauer  pointed  out  that  the  criterion  of  "small  dis- 
placements from  equilibrium"  appears  to  be  a necessary  but 
not  a sufficient  condition  for  the  conclusions  reached  by 
Manes.  C] early  for  a perturbation  no  matter  how  small, 
which  involves  a metastable  interraediate,  tha  rate  of  re- 
turn to  equilibrium  would  not  be  proportional  to  the  dis- 
placement from  equilibrium,  as  measured  by  the  difference 
in  free  energy.  Certainly  a single  relaxation  time  would 
not  suffice  to  characterize  the  path  followed  on  equilibra- 
tion for  such  a system.  Bauer  suggested  that  a sufficient 
condition  would  be  either  that  (a)  one  consider  only  old 
systems  - i.e.,  systems  in  which  the  perturbation  occurred 
sufficiently  long  ago  so  that  the  equilibration  rate  has 
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settled  down  to  its  final  decay  rate,  or  (b)  the  pertur- 
bation must  be  induced  in  that  part  of  the  system  which  is 
associated  with  a process  having  the  longest  tine  constant. 

In  the  discussion  which  follows  we  shall  therefore 
only  postulate  the  applicability  of  a first-order  decay  law 
to  the  equilibration  process  for  adjacent  steps  in  a sim- 
plified mechanism.  Since  the  enthalpy  input  is  continuous 
rather  than  instantaneous  the  mathematical  analysis  is  so.ne- 
what  cumbersome  and  is  dealt  with  in  detail  in  Appendix  I. 
For  a perturbation,  h(t),  the  departure  from  ecuilibrium, 

( t) , is  given  by 


-t/r 

(t)  = e r 


4. 

f ^h(t* 

j at* 


I 


tVt  ♦ 
e dt  , 


(53) 


where  h(t)  =0  for  t < 0,  and  r is  the  ratio  of  "departure” 
to  "perturbation"  which  has  been  assumed  to  be  additive  and 
linear. 

In  the  impact  tube  experiment  the  perturbation  consists 
of  conversion  of  energy  of  mass  flow  into  enthalpy  of  the 
gas,  and  follov/ing  equation  (50) 

(54; 

Substitution  of  (54)  into  (53)  yields 


^h(  t ) _ M du/ 
<?t  '"2  dt 


■n(t)  = I 


1 


du‘ 

d7" 


t*/r 

dt 


(55) 


This  is  the  generalized  derivation  of  Krntrowitz's  equation 
(11)  or  Griffith's  equation  (28). 
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At  this  point  a comparison  between  the  impact  tube 
experiments  and  shock  tube  experiments  (see  Section  I, 
Parts  B.b.  and  C.c.)  is  easily  drawn.  The  shock  tube  ex- 
periments are  clearly  one  limit  of  the  above  description; 
that  is,  in  that  case  h(t)  is  a step  function,  and  corres- 

o 

ponds  to  the  instantaneous  compression  limit,  h , with 


The  effect  of  approaching  such  a limit  is  illustrated  by 
the  histograms  of  figure  7.  For  a given  T this  liinit  may 
be  approached  by  going  to  larger  u,  or  smaller  d and  thus 
decreasing  the  tine  over  which  the  enthalpy  injection 
occurs.  Or,  for  a fixed  value  of  d/u, , one  nay  consider 
systems  having  increasing  't  in  which  case  in  the  limit 
the  enthalpy  injection  occurs  over  times  infinitesimal 
compared  to  T.  If  we  define  the  characteristic  com- 

pression tine  for  the  experiment,  as 

Xo  ^ d/u,  , (57) 

then  the  impact  tube  exi.crimout  approaches  the  shock  tube 
experiment  as  Xq/X  — ^ 0. 

It  is  also  clear  thrt  (55)  is  simplified  considerably 
in  this  limit,  but  at  the  same  time  this  analysis  points 
to  the  obvious  restriction  of  shock  tube  experiments  to 
Instantaneous  perturbations,  which  must  be  at  least  as 
strong  as  the  wealvest  shock  which  can  be  accurately  measured. 
In  turn,  this  means  that  the  departures  from  equilibrium 
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for  the  shock  tube  experiment  will  be  large  compared  to 
those  encountered  in  the  ordinary  range  of  the  impact  tube 
experiment  where  equilibration  proceeds  apace  with  per- 
turbation. For  the  impact  tube  experiment,  the  adjustment 
of  u also  permits  considerable  variation  in  the  magnitude 
of  the  total  entropy  injection. 

One  must  now  formulate  expressions  for  the  departures 
from  equilibrium  in  terms  of  the  tine  and  kinetic  parameters 
describing  the  rate  of  vibrational  energy  and  reaction  equi- 
libration. 

Define  the  following  measures  of  departure  from  equi- 
librium: 

= excess  energy  in  non-lagging  modes 
over  final  ecuilibrium  value 
''^0  = excess  energy  in  lagging  (vibrational) 
modes  over  final  equilibrium  value 
= excess  energy  in  reaction  node  over 
final  ecuilibrium  value. 

It  is  clear  that  since  at  final  equilibrium  T = T3, 

0 = Tj,  and  a = 03  one  may  write: 

= Cp  (T-T,) 

■Vlg  = Cg  (e-T,)  (53) 

% =AH°  {n-n,)  . 

The  specific  heat  terms  used  here  are  defined  as 
Cp  ^ {(l-a)Cp(2)  + 

Cq  r ^ -•^‘^9(1)]  * 
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These  are  in  general  functions  of  both  a and  T,  •v^ile  AP® 
js  a function  only  of  T. 

Let  us  now  consider  the  equilibration  process  as  taking 
place  in  two  steps.  That  is,  the  energy  of  mass  icinetic 
motion  which  is  injected  into  the  gas  at  the  face  of  the 
impact  tube  will  be  pictured  as  going  first  into  trans- 
lation, from  which  it  flows  into  the  vibrational  modes  and 
then  into  the  chemical  degree  of  freedom.  Diagramraatically 


/ Translational  + \ ^ 

'Rotational  Energy‘S 


/Vibrational  \ / Chemical  \ 

' Energy  ' 5=~  ' Energy  ' * 

^2 


This  picture  is  further  clarified  by  reference  to  figure 
8.  Here  T;j_  is  an  intermediate  t'^mperature  corresponding 
to  the  maximum  temperature  which  would  be  reached  by  the 
vibrational  modes  if  the  chemical  reaction  lagged  indefi- 
nitely, From  the  above  it  follows  that 

dT|(TI  y 

*dT  ~ ^T  ”^0 


^ = -kl  + k2 


(60) 


= -kj  + k2  - k\  + ki  . 

These  three  simultaneous  differential  equations  may  be  ar- 
ranged in  the  follov/ing  pattern: 

(^  + k,)'n^  + (-kI.)'>^Q  +0  =0 

(-k,)  ^2  + k,.>y^  + (-kz)'^^  = 0 

0 + (-^  + kp  = 0 • 


(61) 


Temperature 


Figure  8 

Relaxatiori  Curves  following 
on  Instantaneous  Compression 

Case  of 


I 
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Nov  if  one  introduces  tric.l  solutions  of  the  form 


'ii  (t) 


T,e,a 


> 


(6;) 


and  substitues  these  into  (61),  then  the  condition  tliat 
such  solutions  exist  is 


ki-p  -k|  0 

-k,  kj+kl-p  -kj 

0 -kj  kj-p 


= 0 , (63) 


the  roots  of  which  are  the  Pj's  ?n  (62),  This  means  that 
for  a system  of  the  type  bein^  discussed  there  are  three 
effective  relaxation  times:  l/pi , I/P2,  1/P3« 

At  this  point  suppose  we  simolify  the  kinetic  assump- 
tion (60)  by  wi’iting 


dt 


- 


(64) 


This  clparly  implies  that  for  each  1 only  one  of  the  Aj^j’s 
of  equation  (62)  has  an  appreciable  magnitude  depending 
upon  the  form  of  the  initial  perturbation.  Generally,  this 
is  not  valid;  in  our  case,  due  to  the  different  natures  of 
the  two  equilibration  steps,  it  is  not  a poor  approximation. 
’''\^(t)  is  given  as  a result  of  the  folding  operation  by 
equation  (55),  and  for  an  instantaneous  compression,  h°  is 
given  by  (56);  or,  its  equivalent 

(65) 


—I 


h°  = Cp(Tj-I,)  . 
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C 


The  instantaneous  dlsequilibrations  produced  by  such  an 
instantaneous  compression  are 

"^T  = CpCTi-T,)  "ng  = Cg(T,-Tj)  3 /ai'>(a,-a3).  (66) 

Hence  on  dividing,  one  gets  for  the  ratio  r. 


X I2-T, 


^ T^Cg 
9 Tj-T,  c' 


= (67) 


o 


If  one  now  defines  the  following  dimensionless  vari- 
abl  es: 

t'  ; t4-  , , u'  = u/u,,  ^^(t)  f , 

(63) 


then  (55)  may  be  rewritten  as 

r,  ^ -v/r[ 


^i(t) 


M ii  u 

2 


:’/t;  ' 

L 


du<^ 

dt’ 


dt» 


(69) 


Note  that  the  integral  depends  only  on  and  on  the 
shape  (not  the  size)  of  the  impact  tube. 

These  formulas  may  now  be  applied  to  the  calculation 
of  AS, 3.  Equation  (5I)  may  be  approximated  to  first  order 
as 


dt 


4 (a-T)  ^ ^ (T-T^)  |S 

dt  dt 

|T«- 


. To-'n  , 

since  ^ <<1* 


(70) 


ind 


I ^^|«1  . That  is. 


T I T 

' -^o  ' 'o’  '*0 

(70)  follows  from  (51)  by  neglecting  squares  and  higher 

powers  of  the  above  terms  in  the  expansion.  T represents 

an  average  temperature  and  may  be  taken  as  Tq  to  this  same 


accuracy. 
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Conbining  (5S)  with  (6a)  yields 


(T-Tq)  = \/Q 


—I 


(T-e)  = ^ ^ 

c 


Cp 


e 


dt 


^0  d 


da  1 
dt  = aH°  dt 


1_  126 
0 9 


1 

AH®  r. 


(71) 


Substitution  of  these  equations  into  (70)  and  then  into 


(52),  yields 
00 


’^T  'le  '^9  ■’^T  ^a  , 

-r-z + — ; V dt  . 

p \ '^a 


(72) 


With  the  aid  of  this  (63)  reduces  to  dinierisionless  varicbles: 

d 1 (Ta-T3)C9(T,-T3)  / , 


AS 


03 


u. 


rv 


''8 


dt' 


d 1 Cq(T,-T3)  f 


dt’ 


d 1 (T2-T3)AH® 


Ui  T‘ 


X 


a 


(ch  -<’W  / ^ 


'Vt  7a 


dt' 


(7 


O 


This  soives  the  problem  of  computing  ASqs  as  a function 
of  the  kinetic  parameters,  since  (69)  may  be  used  to  evalu- 
ate “^^(t)  for  a square-ended  impact  tube  (figure  6),  and 
these  ''?^(t)  may  then  be  inserted  into  (73)  to  compute 
ASq3  • 
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The  complexity  of  the  calculations  is  reduced  con- 
siderably if  one  assumes  that  = "^0  = • This 

implies  that  “ ^0  “ that,  by  virtue  of  (64), 

the  ratio  of  the  rates  of  approach  to  zero  of  two  >Cj^(t)’s 
is  proportional  to  the  ratio  of  their  displacements;  i.e., 

anAt)/dt 

d.ij(t)/dt  = • (74) 

This  is  to  all  appearances  a reasonable  assumption  to  use 
as  a first  approximation.  (Note  that  this  assumption  aJ so 
removes  the  restriction  imposed  by  (64),  which  implied  a 
restriction  on  the  relative  magnitudes  of  the  of  (62), 
since  the  terms  add  for  this  less  gcnercal  case.)  Refer- 
ence to  (69)  and  the  values  of  the  constants  therein  as 
given  by  (66)  and  (67)  shows  that  this  assumption  also  re- 
quires that 

"^x(t)  =n’(t)  . (75) 

This  clearly  follows  from  (74)  on  substitution  of  the  di- 
mensionless forms.  Hence,  the  consequence  of  assuming  equal 
't'»s  is  that  (73)  may  be  rewritten  as: 

/iSo3  = ^ iy  [Cg(T3-T,)(Ta-T,)  - AJi'’(':a-T3)(a,-03)]  X 

'V7»(t»)  dt»  . (76) 

o 

An  interesting  result  may  be  obtained  by  going  back  and 
computing  ASo3)<  - under  the  sajne  degree  of  aporoximation. 
For  the  instantaneous  compression  case  the  initial  magni- 


o 
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tudes  of  disequilibrium  are  given  by  (66),  and  for 

gets  for  their  time  dependence 
-t/T^ 


^l(t)  = Cp(T2-T3)  e 
Ue(t)  = Cq(Ti-T3)  e 


(77) 


O * t / 

"»Z^(t)  = /!SH°(ai-a3)  e 


Substitution  of  these  equations  into  (72)  yields 

eo 


’'’=>1.0.  '¥  I 


— t 


C^(T2-T3)Ce(T,-T3) 

M ^ 


- 2/ 
^0 


2 


(Ti~T3)  ^ C^(T2-T3)AH°(a,-cx3) 

+ _ j « 


-2t/r 

e dt, 


or 


^os)i  c " [cg(T3-T,)(Tj-T,)-aH»(Tj-T3)(a,-a3)]  | . 

(73) 

Comparison  of  (73)  with  (76)  now  permits  one  to  write 


AS' 


00 

= / 


Tl'(t')  dt'  , 


(79) 


o. 


\diere 


AS 


* - AS03/A303).. 


(30) 


Equation  (79)  is  identical  in  turn  with  Griffith's 
equation  (29)  or  Kantrowitz*  equation  (I6).  One  sees  that 
AS’  depends  only  on  for  a given  flow  pattern.  This  is 


the  simplification  which  follows  from  the  fact  that  in  the 
impact  tube  experiment  only  small  departures  from  equi- 
librium are  introduced. 

The  method  to  be  used  for  analyzing  the  data  is  now 
quite  clear.  Experimentally  one  determines  (Pq-Pz)  for  a 
given  set  of  flow  conditions.  This  corresponds  to  an 
entropy  increase  which  is  to  be  compared  with  the  instan- 
taneous compression  entropy  increase  computed  for  the  same 
set  of  flow  conditions.  The  resulting  value  of  A.S'  is 
converted  to  a value  of  ' using  equation  (79),  and  from 
this  vclue  and  the  known  ratio  d/u,  the  relaxation  time, 

'V , may  be  computed. 
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III.  Equipment  and  Procedure 
A.  Equipment 

The  experimental  procedures  required  for  the  impact 
tube  exceriment  are  of  three  types.  They  involve  the 
preparation  in  adequate  amounts  of  the  pure  gases,  the 
setting  up  of  a gas  circulation  system  \vhich  maintains 
the  desired  non- turbulent  flow  under  conditions  such  that 
the  temperature,  pressure,  and  velocity  of  the  streaming 
gas  are  at  known  and  controllable  values,  and  the  con- 
struction of  an  impact  tube  with  attendant  pressure 
sensing  device. 

Both  the  carbon  dioxide  and  the  nitrogen  tetroxide 
used  in  these  experiments  were  purchased  in  tanks  as  the 
liquified  gases  from  The  Matheson  Company. 

Adequate  drying  of  the  carbon  dioxide  presented  a 
serious  problem.  Eventually  a high  pressure  tank  was 
constructed  in  which  the  gas  to  be  dried  was  allowed  to 
come  to  equilibrium  with  phosphoric  anhydride  at  a pres- 
sure of  about  60  atmospheres.  The  carbon  dioxide  used  in 
the  reported  results  was  Matheson  "Bone  Dry*'  grade  (purity 
specified  at  99.956|4;  remainder  water).  Twelve  pounds  of 
this  gas  were  held  over  three  pounds  of  reagent  grade 
phosphoric  anhydride  for  one  month  at  60  atmospheres  be- 
fore being  used.  Due  to  the  pronounced  dependence  of  the 
relaxation  time  of  carbon  dioxide  on  moisture  content,  the 
consistency  of  our  data  is  good  Indication  of  the  purity 
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of  the  gas  used.  Of  course,  there  was  no  way  to  prevent 
the  gas  from  picking  up  minute  amounts  of  impurities  during 
flow,  from  the  tubing,  valves,  etc, 

..liile  the  carbon  dioxide  flowed  readily  from  the  tank 
at  room  temperature,  such  was  not  the  case  for  the  nitro- 
gen tetroxide  vrtiich  boils  at  21.3°C.  and  has  a heat  of 
vaporization  of  99.0  cal./g.  at  il.0°C.  In  order  to 
assure  a sufficiently  rapid  flow  of  this  gas  it  was  neces- 
sary to  heat  the  tank.  The  heater  was  constructed  from 
three  600w.  glow  coil  elements  using  for  a core  a brass 
tube  5 1/4"  inside  diameter  and  1?^"  long,  with  asbestos 
and  sand  as  the  insulating  materials.  The  heater  could 
be  operated  at  600,  900,  1200,  and  ISOO  w.  input,  and  the 
air  temperature  between  the  cylinder  and  heated  jacket 
was  roughly  250°C.  during  the  course  of  most  runs.  It  was 
observed  that  the  gas  flow  increased  rapidly  when  the 
nitrogen  tetroxide  reached  the  vicinity  of  its  critical 
temperature  of  15-'^.2°C.  The  gas  was  taken  from  No.  3 
cylinders  (each  having  a capacity  of  10  lbs.)  and  was 
specified  as  having  a minimum  purity  of  99,^  with  only  0.01^ 
of  water  and  traces  of  lower  nitrogen  oxides.  A record 
was  kept  as  to  which  cylinder  was  used  for  each  run,  for 
although  each  cylinder  gave  seif  consistent  results  there 
appeared  to  be  a consistent  difference  between  runs  made 
using  different  cylinders.  Even  with  the  above  arrange- 
ment the  highest  pressures  reported  were  reached  only  with 
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difficulty  and  the  data  obtained  suffered  sonewliat  from 
this  limitation. 

After  leaving  the  tank  both  gases  were  passed  through 
an  "after-heater”.  In  the  case  of  the  carbon  dioxide  this 
heater  was  used  to  compensate  for  the  cooling  iue  to  the 
Joule-Thorason  effect.  For  the  nitrogen  tetroxide  it  was 
needed  to  insure  complete  vaporization.  From  here  the 
gas  passed  through  tygon  tubing  to  the  system  past  a 
thermocouple  and  through  a bubbler,  the  latter  serving 
both  to  keep  the  input  pressure  from  greatly  exceeding 
atmospheric  and  as  a simple  indicator  for  adjusting  the 
gas  flow.  (See  figure  9,  vriiich  also  shows  the  equipment 
involved  in  the  following  discussion.) 

To  sweep  the  apparatus  free  of  moisture,  for  purposes 
of  ciligament,  and  the  taking  of  null  readings  a tank  of 
nitrogen  was  also  attached  to  this  input.  The  nitrogen 
used  was  Cerford  "dry  nitrogen"  from  Airco,  specified  as 
99.99%  pure  with  Hj  as  the  sole  impurity. 

The  gas  then  passed  through  approximately  10  feet  of 
glass  tubing  before  entering  the  first  of  two  12L.  reser- 
voirs, which  were  connected  through  a pressure  regulator. 
When  the  "after-heater"  was  adjusted  so  that  the  input 
temperature  of  the  gas  was  very  nearly  room  temperature 
the  time  of  holdup  in  the  connecting  tubing  and  reservoirs 
was  sufficient  to  insure  that  the  gas  had  come  to  thermal 
equilibrium.  When  working  with  nitrogen  tetroxide  care 
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was  taken  to  have  the  room  temperature  stabilized  at  25°C., 
or  slightly  higher,  thus  insuring  no  condensation. 

The  pressure  regulator  was  of  the  diaphragm  valve 
type,  being  a Reynolds  gas  regulator  utilized  commercially 
for  maintaining  constant  flame  height  despite  variation  in 
pressure  in  the  gas  mains.  An  arbitrarily  variable  pres- 
sure of  a standard  reference  bulb  was  substituted  for  the 
’’atmospheric”  side,  and  hence  it  was  possible  to  set  the 
regulator  so  that  the  pressure  in  the  second  12L.  bulb  was 
held  constant  (to  ±0.1  mm.  Hg)  at  any  value  less  than  our 
supply  pressure  of  one  atmosphere.  A serious  difficulty 
arose  because  of  the  corrosion  of  the  diaphragm  and  valve 
by  the  nitrogen  tetroxide.  Eventually_ it  was  found  that  a 
diaphragra  cut  from  an  0.010”  thick  sheet  of  Visqueen 
polymer  had  the  necessary  flexibility  and  was  only  slowly 
attacked  by  oxidizing  agents.  A portion  of  the  valve  was 
constructed  of  stainless  steel  vdiich  appreciably  reduced 
leakage  due  to  corrosion. 

From  the  second  reservoir  bulb  the  gas  passed  through 
a stopcock  and  then  into  the  Impact  apparatus  j^roper.  This 
equipment  is  shown  in  cross-section  in  figure  10,  with  the 
gas  enter ir.g  at  H.  In  this  apparatus  the  gas  flows  up  the 
glass  tube  and  over  a thermocouple  which  measures  the  tem- 
perature, Tq,  before  expansion.  In  all  runs  this  tempera- 
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ture  was  very  nearly  25®C.  and  differed  by  at  most  a couple 
of  degrees  from  room  temperature.  Comparison  of  this 
thermocouple  with  that  in  the  bubbler  gave  a measure  of 
the  extent  of  temperature  change  undergone  by  the  gas  in 
the  equipment. 

The  gas  next  passed  through  a 200  mesh  monel  metal 
screen  sandwiched  between  two  pieces  of  coarser  ■'nonel 
screening,  (initially  copper  screening  was  used,  but  cor- 
rosion by  the  nitrogen  tetroxide  was  quite  rapid.)  This 
acted  as  a field  of  point  sources  and  led  to  a smooth 
laminar  flow  up  the. polished  stainless  steel  tube  to  the 
nozzle.  Part  way  up  this  tube  small  holes  ooniiected  it  to 
a jacket  so  that  the  input  pressure  could  be  applied  to  a 
manometer  which  measured  the  initial  pressure,  p^;  this 
pressure  could  also  be  cippiied.  to  the  outside  of  a metal 
bellows  as  indicated  in  the  drawing.  IVhen  steady  state 
conditions  were  established  no  gas  flo-vred  through  these 
holes  ^nd  therefore  their  presence  clearly  could  not  effect 
the  flow. 

In  designing  the  nozzle  through  uhich  the  gas  flows 
into  the  expansion  chamber  it  is  important  to  make  the 
time  rate  of  energy  change  as  small  as  possible  so  that 
complete  equilibrium  is  ii-aintainea  at  all  times.  From  the 
equations  of  section  IIB  it  ?s  clear  chat  d-S/dt  is  mini- 
mized when  dT/dt  is  smallest;  that  js,  vh'en  the  accelera- 
tion occurs  with  a small  rate  of  temperature  drop.  For  the 


work  with  carbon  dioxide  a 1/2’'  long  nozzle  shaped  like 

those  of  previous  investigators  was  used.  For  the  work 

with  nitrogen  tetroxide  a new  noz.-’J.e  was  constructed.  This 

7 o 

was  designed  following  the  equations 

ur  = cons’ z. 

6 

r X = cons’ t.  , 

where  r is  the  radius  of  the  hole  and  u the  flow  velocity 
of  the  gas  at  a distance  x from  the  base,  of  the  nozzle. 

The  design  was  based  on  a length  of  1.6"  and  an  orifice 
diameter  of  0.0700"  before  introducing  the  boundary  layer 
correction.  The  outlet  velocity  was  taken  to  be  17.2  ni./sec. 
and  the  von  Doenhoff  bounaary-layer  displace'.;ent  thickness 
was  computed  using  a value  of  120x10“^  poises  for  the  vis- 
cosity and  an  average  density  of  1.279x10“^  gm./cm.^. 

These  values  correspond  to  optimum  conditions  for  p^  = 

0.500  atra.,  Pq/Pi  = 1*44,  Tq  = 25°C.  The  find  nozzle 
dimensions  were; 
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0.01 
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0.60 
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0.03 
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0.30 

.0397 

0.07 

.0590 

1.00 

.0385 

0.10 

.0556 

1.20 

.0376 

0.20 

.0495 

1.40 

.0370 

0 . 4^ 

. 0A4<^ 

1.60 

.0366 

Corrosion  of  the  brass  nozzle  used  with  carbon  dioxide  was 
found  to  be  quite  rapid  when  nitrogen  tetroxide  was  used, 
and  the  nozzle  could  not  be  gold  plated  witnout  running  the 
risk  of  introducing  rough  spots  in  the  carefully  shaped 
throat.  An  attempt  to  machine  the  longer  nozzle  of  dimen- 
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sions  noted  above  out  of  stainless  steel  proved  frustrat- 
ingly  unsuccessful  due  to  breakage  of  the  reai.;ing  tool. 

It  was  observed,  however,  that  duraluininuia  (44ST)  with- 
stood the  action  of  th:-  gas  for  a considerable  period. 

Hence  the  nozzle  was  iiade  out  of  a stainless  steel  jacket 
surrounding  a duraluminum  plug,  in  which  the  throat  had 
been  machined.  V/ith  this  nozzle  no  corrosion  problems 
were  encoiintered. 

The  expansion  chamber  was  pumped  on  at  the  point 

marked  ”gas  outflow”  in  figure  10  by  a two  stage  hi-vac 

* 

pump  of  very  large  capactiy.  It  is  necessary  that  the 
pump  used  should  be  able  to  maintain  a pressure  drop  of  a 
fevctor  of  two  across  the  nozzle  when  p^  = 1.0  atm.  For 
the  nozzle  already  described  this  requires  a pumping  rate 
of  approximately  20  L.  (S.T.P.)  per  minute.  A mercury 
manometer  was  also  attached  to  the  expansion  chamber  at 
this  point  to  permit  me? surement  of  p, . It  was  found  that 
for  a fixed  Pq  any  expansion  ratio  could  easily  be  selected 
by  adjusting  a stopcock  interposed  between  the  chamber  and 
the  pump. 

A thermocouple  was  also  mounted  in  the  expansion 
chamber  alongside  the  region  of  smooth  flow.  In  later 
designs  this  was  removed  as  it  served  no  useful  purpose, 
and  always  gave  essentially  the  same  value  as  did  the  in- 
put thermocouple. 

_ 

Manufactured  by:  International  Machine  Works,  Worth 

Bergen,  N.  J. 
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The  requireaents  for  the  impact  tube  are  cuite  si.nple, 
however,  considerable  practice  is  required  to  meet  them. 

They  are: 

i)  The  tube  imist  be  round  and  have  the  center  hole 
concentric  with  the  outei  edge, 

ii)  The  ratio  of  the  inside  diameter  to  outside  dlemeter 
should  be  close  to  0.620  since  this  is  the  value 

4 4 

for  which  theoretical  calculations  are  avaiirble. 

iii)  The  tip  surface  must  be  at  right  angles  to  the 
axis  of  the  tube. 

iv)  The  tube  should  smoothly  widen  to  its  maximum 

diameter,  but  this  must  occur  over  a short  distance 
so  that  the  gas  flow  through  the  tip  to  reach  equi- 
librium will  not  be  greatly  impeded. 

The  tubes  may  bo  constructed  from  ordinary  4 iom.  pyrex 
tubing  by  drawing  out  in  tvro  steps.  In  the  first  step  the 
tube  is  drawn  down  to  about  1/4  of  its  initial  diameter  and 
the  walls  thinned.  In  the  second  step  the  tube  is  then 
drawn  down  to  the  required  outside  diameter  of  from  0.02 
to  0.08  cm.  Two  steps  are  required  to  get  the  proper  wall 
thickness  (otherwise  only  very  low  i.d./o.d.  can  ordinarily 
be  obtained).  It  is  quite  important  that  the  tube  be  heat- 
ed evenly  if  a circular  cross-section  is  to  result.  If  the 
glass  is  twisted  while  being  drawn  the  symmetry  is  consider- 
ably improved  and  the  resultant  tip  is  more  often  coaxial 
with  the  base  of  the  tube.  The  tubes  may  then  be  broken 


crystal  of  carborun- 


off  by  scratching  slightly  with  a s.Tiall 
dum  and  applying  pressure,  llith  considerable  luck  as  well 
as  practice  this  breaking  can  be  done  so  as  to  yield  a 
perpendicular  non-chipped  surface,  liore  often  it  is  neces- 
sary to  grind  the  tip  flat  and  perpendicular  to  the  tube 
axis  on  a very  fine  stone.  V»‘ashing  the  tubes  with  a strong 
solution  of  detergent  followed  by  large  quantities  of  dis- 
tilled water  usually  proved  to  be  sufficient  cleaning.  They 
were  then  checked  and  measured  under  a microscope  with  a 
graduated  eyepiece  reading  to  ±0.0002  cm.  One  could  easily 
tell  wiiether  the  surface  was  really  smooth  and  the  inside 
and  oiitsido  surfaces  coaxial.  Measurements  of  both  di- 
ameters were  taken  at  two  positions  separated  by  90  degrees 
to  check  for  eccentricity.  The  impact  tube  was  then  mount- 
ed with  black  wax  in  the  brass  tube  projecting  from  the 
movable  stage  A so  that  it  was  connected  to  the  insidf’  of 
the  bellows  previously  mentioned.  A small  "sqaare"  was 
used  to  set  the  impact  tube  perpendicular  to  the  ground 
surface  at  B.  The  extension  of  the  tube  below  the  ground 
surface  was  adjusted  so  that  the  tip  was  between  one  and 
two  times  the  nozzle  diameter  from  the  nozzle  face,  as 
this  distance  was  found  by  Griffith  to  give  the  most  con- 
sistent results. 

The  only  portion  of  the  eqv'.ipment  remaining  to  be 

described  is  the  differential  pressure  sensing  device 

that  portion  of  the  equipment  which  measures  the  difference 
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between  the  initial  presrure,  Pq,  and  the  final  pressure, 

Pz,  -vdiich  exists  at  the  face  of  the  inpact  tube  where  the 
gas  has  been  brought  to  rest.  These  two  pressures  are 
applied  respectively  to  the  outside  and  inside  of  a stain- 
less steel  bellows.  The  bellows  used  had  a wall  thickness 
of  0.0085",  was  1 1/2"  long  by  3/4"  in  diameter,  and  had 
23  convolutions.  Its  flexibility  was  such  that  a nressure 
differential  of  one  atmosphere  produced  a change  in  length 
of  about  0.05".  Most  of  the  inside  of  this  bellows  was 
filled  with  a brass  plug  so  that  very  little  gas  need  flow 
through  the  impact  tube  to  adjust  to  a given  pressure  de- 
fect. Thus  its  time  constant  was  minim.ized.  The  top  of 
the  bellows  was  closed  with  a.  stainless  steel  plug  which 
carried  a stainless  steel  rod  1 1/2"  long  with  the  core  of 
a differential  transformer  mounted  at  its  tip.  This  entire 
unit  was  enclosed  in  a brass  shell  which  was  connected  to 
the  jacket,  and  thus  maintained  at  the  pressure  p^.  On  the 

outside  of  this  shell  and  around  the  soft  iron  core  the  dlf- 

♦ 

ferential  transformer  was  mounted.  As  nay  be  seen  in 
figure  10,  the  central  winding  of  the  transformer  was  power- 
ed by  a Sola  Voltage  Regulator  (1.04  amp.,  115  v.  output) 
which  was  fed  by  a Stabline  Voltage  Regulator.  A combi- 
nation of  resistances  was  chosen  so  as  to  load  the  Sola  at 
near  its  rated  capacity,  and  to  provide  7.0  v.  as  a safe 

* Schaevitz  Engineering,  Camden  11,  N.  J. 

Series  S,  Type  18L 

Sensitivity  0.0008v/.001"/volt  in.  at  60  CPS 


continuous  operating  input  for  the  differential  transformer. 
Sixty  cycle  current  was  utilized  in  preference  to  a higher 
frequency  (which  gave  a higher  output)  because  no  oscil- 
lator which  had  sufficiently  good  amplitude  stability  was 
available.  The  outputs  of  the  two  secondary  coils  of  the 
differential  transformer  were  connected  in  opposition  to  a 
Western  Electric  107A  battery  operated  amplifier  (gain 
circa  800)  which  was  shock  mounted  to  reduce  troublesome 
microphonics.  The  amplifier  output,  taken  from  the  primary 
of  the  output  transformer  (no  load  on  secondary) , was  fed 
into  a Simpson  meter  (model  260) . This  meter  responds  when 
set  in  the  ”a.c.  volts"  position  by  rectifying  the  applied 
voltage  and  applying  the  resultant  direct  current  to  a 
suitably  calibrated  microamiaeter.  3y  tapping  across  this 
microarameter  one  obtains  the  input  current  suitably  recti- 
fied for  use  on  a recorder,  and  at  the  same  tij.ie  the  meter 
serves  both  as  a rough  indicator  for  adjusting  the  dif- 
ferential transformer  to  minimum  position  and  as  a sensi- 
tivity control.  As  actually  used,  the  level  of  the  recti- 
fied voltage  was  controlled  at  the  Simpson  meter  and  its 
zero  was  set  by  means  of  a bucking  circuit;  the  net  voltage 
was  finally  recorded  by  a Brown  high  speed  instrument  (Model 
No.  153X12V-W6-30,  full  deflection  10  mv.). 

In  addition  to  being  very  poisonous,  nitrogen  tetroxide 
is  a highly  corrosive  gas.  Insofar  as  was  possible  all 
equipment  subject  to  its  attack  v^as  constructed  of  glass  or 


stainless  steel.  Another  metal  which  resists  its  eetJon 


qnite  well  is  duraluminum  (2.‘.nT),  which  has  an  ndventafie 
over  stainless  steel  in  being  much  easier  to  machine.  In 
several  instances  gold  plated  brass  was  used  in  our  equip- 
ment, and  this  proved  resistant  to  everything  except 
direct  exposure  to  a constant  stream  of  the  gas.  Stou- 
cocks  were  lubricated  ;n,th  a halocrrbori  grease,  which  re- 
sists attack  for  some  time,  but  this  grease  has  a high 
viscosity  at  room  temperature.  In  the  case  of  those  ground 
surfaces  where  the  gas  had  very  little  opportunity  for 
direct  attack  Apeizon  grease  (L)  war.  found  to  be  satis- 
factory over  long  periods.  The  metal  and  glass  parts  were 
sealed  together  using  black  Apeizon  wax.  In  designing  the 
equipment  one  must  be  careful  that  such  joints  are  well 
protected  from  direct  blasts  nitrogen  tetroxide  ^.rtiich 
oxidizes  exposed  regions  to  produce  a porous  crumb] y mass. 
V.'ith  suitable  care  in  design  such  joints  last  quite  well, 
however.  Tygon  tubing  was  the  only  clastic  tubing  used  in 
these  experiments,  and  while  it  became  more  rigid  with 
continued  exposure  to  the  gas  no  layer  of  corrosion  rroducts 
was  noted. 


B.  Procedure 


a.  Aligmaent 

There  are  two  alignmeiit  steps  to  be  tciken  before  the 
equipraent  becomes  ready  for  Cciibration.  The  first  is 
adjustnent  of  the  differential  transformer  and  its  core, 
and  the  second  is  positioning  of  the  impact  tube  with 
respect  to  the  flow  stream. 

As  can  be  seen  from  figure  10,  the  position  of  the 
brass  tube  surrounding  the  transformer  core  can  be  varied 
by  using  the  top  ring  of  three ■ thumbscrews  which  move  the 
top  portion  (at  A)  with  respect  to  the  impact  tube  carry- 
ing section.  If  the  core  is  initially  perpendicular  to 
the  base  plate,  and  the  bellows  free  from  strains  so  theat 
it  moves  only  axially  on  comnression,  this  lateral  motion 
permits  centering  the  core  so  that  nowheres  does  it  touch 
the  surrounding  brass  tube.  The  core  may  be  adjudged  to 
be  moving  freely  when  the  calibration  curve  (see  below)  is 
reproducible  and  shows  no  hysteresis.  Mounted  atop  the 
differential  transformer  carriage  is  a differentjal  mi- 
crometer (not  shown  in  detail)  constructed  by  attaching  a 
left  hand  threaded  section  of  42  threads  per  inch  in  op- 
position to  a comiv.ercial  right  hand  micrometer  having  40 
threads  uer  inch.  Uith  this  arrangement  the  probe  is 
advanced  0.0012  inches  for  one  full  turn,  and  one  can 
easily  select  one  three-hundredth  of  this.  Tne  probe 
presses  down  on  the  differential  trc.nsformer  in  opposition 
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to  a spring  mounted  below;  it  is  used  to  adjust  the  core 
to  approximately  J?  mv.  off  null  position  where  the  response 
curve  is  very  nearly  linear  and  the  sensitivity  highest. 

The  initial  adjustment  on  the  impact  tube  which 
determines  its  distance  from  the  nozzle  has  already  been 
described.  T>ie  lower  ring  of  thumbscrews  at  the  top  of 
the  impact  apparatus  (figure  10)  permits  one  to  move  the 
tube  in  the  horizontal  plane,  at  B,  without  disturbing  the 
pressure  sensing  imit.  Lateral  transver ces  showed  that  a 
considerable  fraction  of  the  nozzle  area  was  free  from 
turbulence  since  the  defect  measured  with  nitrogen  was 
quite  small  and  coristant.  Beyond  this  area  the  pressure  in 
the  impact  tube  changed  quite  pronouncedly.  The  angle  of 
attack  was  ad Jus  cable  via  the  lowest  ring  of  thumbscrews 
shown  in  that  diagram;  these  work  oij  t.he  ball-  joint  which 
has  as  its  center  the  tip  of  the  impact  tube.  Vie  found, 
as  had  Griffith,  that  the  observed  defects  were  constant 
a.s  long  as  the  angle  of  attack  was  within  5°  of  zero.  Also, 
it  was  found  possible  to  construct  impact  tubes  so  that 
their  tips  were  so  nearly  perpendicular  to  the  axis  of  the 
tube  that  no  angular  adjustment  was  needed.  Hence,  after 
this  was  ascertained,  the  ball- joint  section  was  replaced 
by  2 straight  glass  column  with  ends  ground  perpendicular 
to  its  axis,  and  with  only  a single  side  tube  for  outflow. 
This  simple^  design  was  also  substantially  more  resistant 
to  breakage  than  its  more  complicated  forerumier  and  re- 
moved a potential  point  of  leakage. 
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b.  Calibration  and  Sensitivity 

A pressure  differential  could  be  applied  across  the 
v-'i'-ssurc  sensing  device  by  opening  the  stopcock  at  the  top 
of  the  apparatus  (figure  lO)  to  the  air,  closing  the  stop- 
cock between  the  bellows’  chamber  and  expansion  chauber, 
and  evacuating  the  rest  of  the  inner  chamber  slightly.  For 
calibration  purposes  an  n-butyl  phthalate  open  end  manoneter 
attached  as  shown  in  figure  9 v;as  used  to  read  the  pressure 
difference.  The  ratio  of  the  density  of  n-butyl  phthalate 
to  that  of  mercury  was  found  lo  be  1:12.30;  hence,  a 
pressure  difference  could  easily  be  read  to  ± 0.02x10"^ 
atm.  (±  2x0.01  c.n.  of  n-butyl  phthalate). 

A typical  calibration  curve-  is  :ho'.vTi  in  figure  li. 

The  points  have  been  nui.ibered  in  the  sequence  in  which 
they  were  taken  (nos.  1,  6,  and  12  were  at  zero  pressure 
differential).  The  stability,  reproducility , and  absence 
of  hysteresis  arc  indicative  of  proper  core  centering  and 
high  stability  in  the  electrical  circuit.  v/il:h  the  Brown 
recorder  set  for  full  scale  deflection  at  10  mv.  it  was 
possible  to  read  to  ± 0,02  mv.,  and  hence  a pressure  dif- 
ference could  be  determined  to  witliin  the  equivalent  of 
± O.O4  mv.  Reference  to  the  calibration  plot  indicates 
that  the  response  was  about  3.3  mv./cm.  of  n-butyl  phthal- 
ate, hence  pressure  differentials  could  be  read  to  a pre- 
cl  Sion  of  '^1.2x10"  atm.  In  some  later  runs  the  sensi- 


tivity was  increased  by  about  35^  by  finer  adjustment. 


Figure  11 

Typical  Calibration  Plot 
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The  greatest  difficult^''  in  this  setup  was  a gradual  wander- 
ing of  the  zerb  point,  which  therefore  had  to  be  redeter- 
mined after  every  few  readings,  so  that  the  true  zero 
could  be  interpolated.  Higher  quality  electrical  com- 
ponents would  doubtless  reduce  this  difficulty  to  the 
point  where  it  would  cease  to  be  significant.  The  re- 
sponse time  of  the  whole  pressure  sensing  device  was  limit- 
ed by  that  of  the  recorder  (2  seconds  for  full  scale  de- 
flection) . 

In  view  of  the  observed  stability  and  reproducibility, 
the  sensitivity  of  this  unit  could  quite  easily  be  increased 
by  using  an  amplifier  with  greater  gain.  Other  possi- 
bilities are:  increasing  the  cross-section  of  the  bellows, 

using  a longer  bellows  with  more  convolutions,  using  a 
bellows  of  more  flexible  material,  using  a recorder  of 
greater  sensitivity,  or  using  a differential  transformer 
of  higher  voltage  output  (step  up  ratio). 

c.  Nitrogen  Background 

In  order  to  check  flow  conditions  one  must  make  a run 
with  a gas  whose  rotational  and  vibrational  relaxation 
times  are  very  much  smaller  than  the  compression  time  of 
the  equipment.  Nitrogen  is  such  a gas  since  the  vibration- 
al contribution  to  its  specific  heat  is  negligible  at  room 
temperature  and  its  rotational  specific  heat  has  a relaxa- 
tion time  shorter  than  10"^  sec.  at  room  temperature  and 

lb 

atmospheric  pressure. 


The  results  of  a typical  run  with 
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nitrogen  are  shown  in  figure  12.  The  sharp  break  at  an 
expansion  ratio  of  about  2.2  is  due  to  the  setting  in  of 
t’lrbulence  and  will  be  discussed  more  fully  in  section  IV. 
The  smooth  curve  and  low  value  of  the  defect  compared  to 
that  for  carbon  dioxide  indicates  a regular  background 
effect  which  can  be  corrected  for.  The  fact  that  the 
nitrogen  defect  curve  is  small  and  linear  below  Pq/pi  = 
2.0  and  passes  through  the  origin  is  in  line  with  the 
fact  that  ordinary  hydrodynamic  effects  such  as  misalign- 
ment of  the  impact  tube,  teciperature  or  pressure  inhomo- 
geneities, turbulence,  etc.  would  be  expected  to  produce 
a total-head  defect  which  varies  directly  as  (Pq/pi  - 1) . 
This  -^is,  of  course,  not  the  case  for  a defect  depending 
on  a heat  capacity  lag. 

d.  Procedure  for  Carbon  Dioxide  Runs 

The  following  stepwise  procedure  was  followed  in 
making  all  observations  on  carbon  dioxide  recorded  in 
section  IV. 

i)  For  greater  stability  the  electrical  equipment 
supplying  the  input  to  the  differential  transformer 
was  not  turned  off  between  runs.  The  system  was 
checked  for  leaks  before  each  set  of  runs  iising  an 
auxiliary  mercury  diffusion  puiup  circuit  with  acLeod 
and  Pyrani  gauges. 

ii)  The  battery  powered  amplifier,  rectifier,  buck- 
ing cjrcuit,  and  recorder  were  turned  on  about  one 
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hour  before  the  start  of  the  actual  run.  The  end  of 
thc-ir  required  war:a-up  period  was  ascertained  by  a 
trece  of  constant  deflection  on  the  recorder.  During 
this  period  the  theruiocouple  cold  junction  was  filled 
and  the  thermocouples  checked,  the  pump  serviced  and 
started,  the  continuously  monitored  input  carefully 
adjusted  to  7.00  v.,  the  position  of  the  transfor..ier 
with  respect  to  the  core  ad  pasted,  the  barometric 
reading  recorded,  p^^  control  set  to  approximately  the 
value  desired,  and  the  room  brought  to  an  even  ttm- 
Dc-rature. 

iii)  A complete  calibration  curve  of  recorder  dis- 
placements versus  n- butyl  phthalrte  manometer  read- 
ings was  taken  for  later  use.  This  curve  also  in- 
dicated by  its  lack  of  hysteresis  that  the  core  was 
well  centered  and  served  as  an  extra  check  on  the 
leak-tightness  of  the  central  portion  of  the  equip- 
ment . 

iv)  A background  run  was  made  with  dry  nitrogen 

over  the  pressure  ratio  range,  1 . 2 (p^/p, ) 4 2 . 0, 

giving  the  magnitude  of  the  hydrodynamic  correction. 
This  also  served  to  check  for  proper  aligi-uaent  of 
the  impact  tube,  permitted  final  adjustment  of  pQ, 
and  swept  the  apparatus  free  of  foreign  gases. 

v)  The  system  was  flushed  with  carbon  dioxide  at  an 
intermediate  pressure  ratio  until  a constant  pressure 
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defect  vas  observed.  The  after  heater  vas  adjusted 
for  proper  temperature  of  the  input  gas. 

Vi)  Measurements  were  recorded  for  the  pressure 
defect  of  carbon  dioxide  over  the  pressure  ratio 
range,  1.2  ^ (Po/pi ) ^ 2.0,  starting  vrith  the  smallest 
ratios.  Theriaocouple  readings  were  taken  at  the 
bubbler  and  input  stations  while  the  gas  was  flowing. 


At  all  times  the  input  voltage  was  monitored  to  check  on 
its  consistency.  The  recorder  readings  at  zero  pressure 
defect  were  noted  before  and  after  steps  lii,  iv,  and  vi. 

If  there  were  appreciable  fluctuations  in  this  null  posi- 
tion, or  if  for  any  reason  a_ suspicion  arose  that  the 
apparatus  was  functioning  not  quite  satisfactorily,  opera- 
tions were  suspended,  and  the  run  startc-d^anew.  Ail  of 
these  steps  were  repeated  at  each  Pq,  the  runs  being  spaced 
by  sufficient  time  so  that  the  batteries  would  not  be  ad- 
versely affected. 


e.  Procedure  for  Nitrogen  Tetroxide  Runs 

The  procedure  followed  here  closely  paralleled  that 
described  above  for  carbon  dioxide.  Adoitional  time  was 
required  for  step  v since  the  tank  of  compressed  gas  had 
first  to  be  brought  to  a high  enough  temperature  to  main- 
tain the  desired  Pq  o.ver  the  range  of  expansion  ratios. 

The  highest  p^  which  could  be  maintained  over  the  whole 
range  of  expansion  ratios  was  0.7  atm.,  and  this  only  with 
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a compar? tively  full  tank.  Afte.’'  the  run  it  was  necessary 
t's  re.-iove  the  heated  tank  from  its  heater  and  cool  it  with 
cloths  soaked  in  ice  water  imtil  evolution  slowed  before 
closing  the  cylinder.  The  corrosive  gas  was  then  sweat 
from  the  equlp..;ent  with  dry  nitrogen. 

Due  to  the  heating  requirements  plus  corrosion  prob- 
lems, runs  with  nitrogen  tetroxide  proceeded  considerably 
less  sr;!!^othly  than  those  with  carbon  dioxide,  and  there- 
fore, as  will  be  seen,  the  results  showed  less  consis- 
tency: Otherwise  the  runs  followed  very  closely  the  came 

pattern,  and  the  added  experience  obtained  from  earlier 
work  helped  to  mitigate  the  increased  difficulties. 
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IV.  Carbon  Dioxide 

In  this  section  heat  capacity  lag  uieasurenent s taken 
on  carbon  dioxide  are  recorded;  their  reduction  to  relax- 
ation times,  and  a critical  interpretation  of  the  results 
are  then  pre;:ented. 


/• . The  Data 

The  purity  of  the  {jas  used  and  the  procedure  followed 
has  already  been  described  in  section  III.  k typicr 1 cali- 
bration plot  is  shown  in  figure  11,  while  figure  IZ  shows 
a typical  plot  of  nitrogen  and  carton  dioxide  defrcts 
expansion  ratio.  The  sharp  breaks  at  (Pq/p,)  = 2.2  corres- 
pond to  the  setting  up  of  shock  waves  at  the  nozzle  and  in 
front  of  the  impact  tube,  when  the  Mach  number  reaches 
unity.  The  decline  in  slope  of  the  pressure  defect  curve 
of  carbon  dioxide  above  (Pq/pi)  = 1.5  appears  to  be  cormaon 
to  all  of  the  runs,  independent  of  p^.  This  has  not  been 
satisfactorily  accounted  for;  it  may  be  due  to  slight  im- 
perfections in  the  throat  surface.  I’jlnts  above  expansion 
ratios  of  1.5  were  not  used  in  the  computation  of  relaxa- 
tion tJmes. 

From  a graph  of  the.  type  shown  in  figure  12  values  of 
the  carbon  dioxide  defect  minus  that  for  nitrogen  at 
selected  exnansion  ratios  were  obtained.  Attempts  to  secure 
useful  data  at  p^  = 0.2  atm.  and  below  were  frustrated  by 
the  fact  that  the  nitrogen  background  correction  then  be- 
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comes  an  appreciable  fraction  50%)  of  the  values  ob- 
served for  carbon  dioxide,  so  that  the  precision  of  the 
data  is  quite  poor.  The  results  of  measurements  at  p^  = 
0.967,  0.752,  0.534»  end  0.341  ?tms.  are  shovm  in  figure 
13,  where  the  relative  pressure  defects  are  plotted  vs 
the  expansion  ratios.  The  impact  tube  used  in  all  of  these 
runs  had  an  outside  diameter  of  0.0770  cm.,  with  (l.D./O.D.) 
=0.56  and  maximum  eccentricity  < 3%, 

It  is  important  tc  note  that  within  the  margins  of 
the  experimental  accuracy  these  curves  are  superposable. 

As  will  be  discussed  later,  this  is  due  to  the  bimolecular 
nature  of  the  collision  process.  It  may  be  noted  here, 
however,  that  such  consistent  data  were  not  obtained  with 
less  pure  or  slightly  moist  carbon  dioxide.  Indeed,  one 
may  consider  this  sort  of  a check  as  a sensitive  criterion 
for  the  purity  of  the  gas.  Previous  investigators  using 
the  impact  tube  method  were  restricted  to  the  fixed  value 
p,  = 1 atm.,  but  could  vary  p^,  and  hence  obtained  only 
one  such  curve.  They  had  no  such  internal  check  on  either 
the  purity  of  the  gas  or  the  bimolecular  nature  of  the 
reaction. 
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B.  Treatment  of  Data 


The  impact  tube  experiment  yielded  values  of  the  pres- 
sure defect  for  various  coiabinaticns  of  Tq,  and  p,  . 
These  data  have  been  used  to  determine  the  relaxation  time 
following  the  theory  developed  in  section  II.  The  experi- 
mentally determined  entropy  gain  was  compared  with  the 
entropy  gain  computed  for  instantaneous  compression  and, 
from  this  ratio,  the  ratio  of  the  unknown  relaxation  time 
of  the  gas  to  the  known  characteristic  compression  time 
of  the  equipment  was  determined.  A stepwise  description 
of  the  method  of  computation  follows;  table  I contair s 
tabulated  values  of  the  more  significant  quantities  com- 
puted. 

1.  The  relative  pressure  defects  v^  expansion  ratios 
were  tabulated  for  each  p^,  T^^  (columns  la,  lb,  Ic). 

2.  Using  the  adiabatic  expansion  equation  (8,',  T,  was 
calculated  for  each  p^,  T^,  and  p,  (column  ?:?>.),  As 

7 1 

a first  approximation  C was  taken  at  T^,  nid  hence- 
forth  at  T = + T,).  Two  successive  .''^rations 

were  sufficient  to  give  T,  to  the  same  accuracy  as  T^ 
(columns  2b,  2c). 

M 2 

3.  The  free  flow  energy,  , in  calories  per  mole  was 

calculated  using 

|u,  ’ = Cp(f)  |vT,  j.  , («1) 

and  the  free  flow  unit  mass  velocity,  u, , in  cm. /sec. 
computed  therefrom  (column  3). 
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4.  Fran  equation  (33)  the  pressure  defect  for  instan- 
taneous compression  war.  computed  (column  4)  • For 

I rj 

carbon  dioxide  C = and  C was  taken  as  the 

difference  between  Cp  and  Cp. 

5.  The  ratio  of  the  experimental  to  the  instantaneous 
entropy  gain  was  calculated  to  a first  approximation 
from 


6 . 


as;  , = ^ 5), 

total  Aoi.,.  RinES),  ^Pi.o/Po 

P3  C82) 


A large  scale  graph  of  the  values  of  ys  computed 
by  Griffith  for  a square  ended  impact  tube  (see  table 
of  values  in  section  IIB)  was  then  utilized  to  obtain 

t 

a first  approximation  to 

7.  A correction  was  then  applied  for  entropy  gain  in  the 
nozzle.  Although  a square  ended  impact  tube  was  used, 
the  correction  term  derived  by  Kantrowitz  and  based  on 
a source  shaped  impact  tube  was  applied.  Since  the 
correction  is  for  the  nozzle  contribution,  the  shape 
of  the  impact  tube  is  not  material.  Only  an  anproxi- 

t 

mate  value  of  ^ is  needed  as  the  term  is  a small 

n 

correction.  The  formula  utilized  was 


As’  = 4 
n 3 Kjj 


16 

K 


9X‘ 


n 


[l-exp(-  f K„)] 

[l-exp(-3Kn)] 


(S3) 


with 
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K - D „ li  p‘ 

" 1 ^ Cp 


(84) 


1 

Here  1 is  the  nozzle  length  for  the  experii-ients  on 
carbon  dioxide),  and  D = 3xradius  of  the  impact  tube. 
Now 

I . t 

, (85) 


Ab 


^total  “ ^ 


+ AS,; 


8. 


where  is  the  relative  entropy  gain  for  the 

t 

overall  process,  i.e.  the  result  of  step  5,  AS^  is 
the  relative  entropy  gain  in  the  nozzle,  and  AS*  is 
the  relative  entropy  gain  on  compression.  Subtracting 
AS^  [eqns.  (83,84)]  from  j eqn.  (82)]  gave  a 

t I 

second  approximation  for  AS  . From  this  a new 
was  obtained.  This  value  was  averaged  vath  the  previ- 
ous approximation  to  and  a new  AS^  computed  there- 
from (column  7a).  TTiis  was  then  subtracted  from 
total  step  5 and  a still  better  approximation  to 
As*  obtained.  This  procedure  was  repeated  until  AS* 
was  determined  to  within  the  accuracy  of  the  experiment- 
al results  (column  7b).  This  required  three  iterations, 
and  in  all  cases  only  the  first  term  in  (83)  was  ap- 
preciable. Only  final  values  of  the  relative  entropy 
gains  are  recorded  in  table  I. 

The  value  of  corresponding  to  the  final  corrected 
value  of  As*  (column  7c)  was  used  to  compute  accord- 
ing to  the  equation 

^ ^ (column  8)  , (86) 

p ’ 


Table  I 


(la) 

Fq/Pi 

(lb) 

Po 

(Ic) 

Pi  ^ 
(atm) 

(2a) 

Ti 

(°K.) 

(2b)  (2c) 

T C (T) 

(°K.)  (cal./aeg. 

D 

0 

= 0.967 

£i  tim  • j 

0 

II 

293°K. 

1.100 

0.11 

0.379 

291.7 

29/f.8 

8 . 34 

1.150 

0.30 

.341 

233.7 

293.4 

3.32 

1.200 

0.54 

.306 

2S6.0 

292.0 

8 . 80 

1.250 

0.33 

.17  A 

233.3 

290.6 

3.79 

1.300 

1.22 

• 744 

230.3 

239.4 

8. 73 

1.350 

1.64 

.716 

278.4 

233.2 

3.76 

1.400 

2.10 

.691 

276.1 

257.0 

3.75 

1.450 

2.60 

.667 

273.8 

285.9 

8.74 

Po 

= 0.752 

Si  • y 

0 

II 

299. 8°K. 

1.100 

0.14 

0.634 

293.5 

296.6 

3.86 

1.150 

0.36 

.653 

290.5 

295.2 

3 . 34 

1.200 

0.63 

. 626 

237.8 

293.8 

3.82 

1.250 

0.97 

.601 

235.1 

292.4 

3.  SI 

1.300 

1 .40 

.578 

2S2.6 

291 . 2 

8.80 

1.350 

1.84 

.557 

230.2 

290.0 

S.73 

1.400 

2.33 

.537 

277 . 8 

283.8 

3.77 

1.450 

2.80 

.519 

275.6 

287.7 

8.76 

Po 

= 0.534 

3 • y 

II 

0 

300. 2°K. 

] .100 

0.17 

0.485 

293.9 

297.0 

3.86 

1.150 

0.37 

.464 

290.9 

295.6 

3.84 

1.200 

0.67 

.445 

23H.1 

94 . 2 

3.33 

1.250 

1.04 

.427 

235.5 

292.3 

3.31 

1.300 

1.46 

.410 

283.0 

291.6 

3.30 

1.350 

1.34 

.395 

230.5 

290.4 

3.79 

1.400 

2.34 

.381 

273.2 

289.2 

8.77 

1.450 

2.32 

.363 

276.0 

283.1 

3.76 

Po 

= 0.341 

A t^d  • y 

To  = 

299. 2°K. 

1.100 

0.16 

0.310 

292.9 

296.0 

3.35 

1.150 

0.34 

.297 

239.9 

294.6 

3.33 

1.200 

0.62 

.284 

237 . 2 

293.2 

3.82 

1.250 

0.92 

.273 

284.5 

291.3 

3.80 

1.300 

1.32 

.262 

232.0 

290.6 

3.79 

1.350 

1.58 

.252 

279.6 

239.4 

8.78 

1.400 

2.12 

.244 

277.3 

283.2 

S.76 

1.45C 

2.43 

.235 

275.0 

237.1 

3.75 

mole) 


Table  I (Cont'd 

.) 

1 

1 

1 

(la) 

(3) 

(4) 

(5) 

(7a) 

(7b) 

Pq/P’ 

Ui  X 10”^ 

(cm. /sec . ) 

Po 

(cfiiorie 

^total 

s/de^. 

mole) 

as’ 

c 

Po 

= 0.967 

atm. 

, Tq  = 298®K. 

1.100 

10.3 

0.27 

0.41 

U.O4 

0.37 

1.150 

12.5 

0.60 

0.50 

0.07 

0.43 

1.200 

14.2 

1.00 

0.54 

O.OS 

0.46 

1.250 

15.7 

1.50 

0.55 

0.09 

0.46 

1 . 300 

17.0 

2.06 

0.59 

0.11 

0.43 

1.350 

13.1 

2 . 66 

0.62 

0.12 

0.50 

1.400 

19.2 

3.33 

0.63 

0.13 

0.50 

1.450 

20.1 

4.04 

0.64 

0.14 

0.50 

Po 

= 0.752 

atm. 

, To  = 299. 

8®K 

• 

1.100 

10.3 

0.27 

0.52 

0.06 

0 . 46 

1.150 

12.5 

0.60 

0.60 

0.08 

0.52 

1.200 

14.2 

1.00 

0.63 

0.08 

0.55 

1.250 

15.7 

1.50 

0.65 

0.09 

0.56 

1.300 

17.0 

2.06 

0.63 

0.09 

0.59 

1 

1.350 

18.1 

2.66 

0.69 

0.10 

0.59 

1.400 

19.2 

3.33 

0.70 

0.10 

0.60 

1.450 

ro.i 

4 . 04 

0.69 

0.10 

0.59 

Po 

= 0.534 

atm. 

, Tq  = 300. 

2°K 

• 

1.100 

10.3 

0.27 

0.63 

0.03 

0.55 

1.150 

12.5 

0.60 

0.62 

0.09 

0.53 

1.200 

14.2 

1.00 

0.67 

0.09 

0.58 

1.250 

15.7 

1.50 

0.69 

0.10 

0.59 

1.300 

17.0 

2.06 

0.71 

0.10 

0.61 

1.350 

18.1 

2.66 

0.69 

0.10 

0.59 

1.400 

19.2 

3.33 

0.70 

0.10 

0.60 

1.450 

20.1 

4.04 

0.70 

0.10 

0.60 

Po 

= 0.341 

atm. 

, Tq  = 299. 

2°K 

• 

1 

; 0 

1.100 

10.3 

0.27 

0.59 

0.07 

0.52 

1.150 

12.5 

0.60 

0.57 

0.07 

0.50 

1.200 

14.2 

1.00 

0.62 

0.03 

0.54 

1 

1.250 

15.7 

1.50 

0.61 

0.03 

0.53 

1 

1.300 

17.0 

2.06 

0.64 

O.OS 

0.56 

1 

1.350 

18.1 

2 . 66 

0.59 

0.03 

0.51 

1 

1.400 

19.2 

3.33 

0.64 

0.09 

0,55 

1.450 

20.1 

4.04 

0.61 

C.08 

0.53 

(la) 

Pq/P- 

Table  I (Cant’d.) 

(7o)  (3)  (9) 

'^K  ^ 10^  P 

^ (sec.)  (atm.) 

• 'X  Tj^xlO 
( atm . sec 

Po 

= 0. 

967  atm., 

To 

= 298°K. 

l.JOO 

0.64 

3.8 

0 . 923 

3.5 

1.150 

0.82 

4.0 

.Q04 

3.6 

1 . 200 

0.92 

3.9 

.886 

3.5 

1.250 

0.92 

3.6 

.870 

3.1 

1.300 

1.00 

3.6 

.856 

3.1 

1.350 

1.08 

3.6 

.342 

3.0 

1.400 

1.08 

3.4 

.829 

2.3 

1.450 

1.08 

3.3 

.317 

2.7 

Po 

= 0. 

752  atm. , 

To 

= 299. 3°K. 

1.100 

0.92 

5.4 

0.718 

3.9 

1.150 

1.17 

5.7 

.702 

4.0 

1.200 

1.32 

5.6 

.689 

3.9 

1.250 

1.38 

5.4 

.676 

3.7 

1.300 

1.57 

5.6 

. 665 

3.7 

1.350 

1.57 

5.3 

.654 

3.5 

1.400 

1.65 

5.2 

.644 

j . 4 

1.450 

1.57 

4.7 

. 656 

3.0 

Po 

= 0. 

534  atm., 

To 

= 300. 2°K. 

1.100  • 

1.32 

7.8 

0.510 

4.0 

1.150 

1.21 

5.9 

.499 

2.9 

1.200 

1.50 

6.4 

.490 

3.1 

1.250 

1.57 

6.1 

.430 

2.9 

1.300 

1.72 

6.1 

.472 

2.9 

1.350 

1.57 

5.3 

. 464 

2.5 

I.4OO 

1.65 

5.2 

.458 

2.4 

1.450 

1.65 

5.0 

.451 

2.3 

Po 

= 0. 

34I  atm. , 

To 

= 299. 2®K. 

1.100 

1.17 

6.9 

0.326 

2.3 

1.150 

1.08 

5 • ^ 

.319 

1.7 

1.200 

1.27 

5.4 

.312 

1.7 

1.250 

1.21 

4.7 

.307 

1.4 

1.300 

1.38 

4.9 

.302 

1.5 

1.350 

1.12 

3.8 

.296 

1.1 

1.400 

1.32 

4.2 

.292 

1.2 

1.450 

1.21 

3.6 

.288 

1.0 

100. 


where  d is  the  diameter  of  the  impact  tube. 

9.  The  mean  pressure  during  relaxation,  p,  was  taken  as 

since  <<i  (coluum  9). 

2 ° (Po-Pi) 

When  the  values  for  each  p^  are  averaged,  one  gets: 


, (atm.) 

^ave • ^ ) 

"^K(ave.)^^^  sec.) 

Ave.  dev 
in 

0.967 

0.866 

3.6 

6^ 

0.752 

0.673 

5.4 

0.534 

0.473 

6.0 

10% 

0.341 

0.305 

4.3 

11% 

The  precision  decreases  with  p^  essentially  due  to  the 
fixed  precision  in  the  differential  pressure  measurement 
by  the  recording  manometer,  and  to  the  proportionately 
larger  nitrogen  background  correction.  The  above  results 
at  higher  p^’ s show  overage  deviations  which  are  smeller 
than  that  for  the  series  of  points  quoted  by  Kantrowitz, 
for  Pq  > 1 atm. 
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C.  Critical  Interpretation 


To  check  on  the  origin  of  the  fluctuations  in  the 
values  of  the  relaxation  times  we  have  thus  obtained,  a 
variety  of  plots  were  made  in  an  attempt  to  pin  down 
obvious  trends  in  the  data.  In  figure  14  values  of 
are  plotted  Vs  Pq/Pi • The  three  lowest  pressures  check 
each  other  fairly  well,  \rtierein  there  is  a strong  sugges- 
ticai  of  a decreasing  trend  with  increasirig  expansion  ratio. 
Relaxation  times  computed  from  the  run  at  p^  = 0.967  appear 
to  be  independent  of  P^/Pi ; these  check  the  data  of 
Kantrowitz  more  closely  (the  latter  were  not  corrected  for 
entropy  loss  in  the  nozzle). 

In  figure  15  the  relaxation  times  are  plotted  ys.  p. 

One  may  conclude  firstly,  for  any  one  run,  there  is  an 
obvious  trend  with  p,  but  the  slope  of  this  trend  decreases 
markedly  with  increasing  pressure.  This  suggests  that 
there  are  deviations  from  laminar  flow  in  the  nozzle, 

>diich  become  aggravated  at  the  lower  pressures  and  higher 
flow  velocities.  It  is  quite  likely  that  the  nozzle  used 
for  these  runs  was  too  short,  and  its  throat  was  not  suf- 
ficiently smoothj  a longer  nozzle,  such  as  that  used  for 
the  work  on  nitrogen  tetroxide,  would  certainly  reduce  the 
correction  term.  Secondly,  the  relaxation  times  decrease 
with  increasing  average  pressure  (as  1/p),  except  for 
values  derived  from  the  lowest  pressure  run.  Assuming  that 
the  transfer  of  energy  is  bimolecular,  one  would  antici;:aie 
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that  the  product,  px  , which  is  proportional  to  the 
"collision  n’jmber",  would  remain  constant  as  long  as  the 
temperature  and  composition  are  not  varied.  This  is  demon- 
strated in  figure  16.  Such  a cheek  on  the  biiuolecular 
nature  of  the  energy  exchange  process  has  previously  been 

7 2 

made  by  means  of  sonic  measurements.  In  this  figure  one 
also  notes  a general  decrease  in  Pk  with  decreasing  p. 

This  trend  suggests  that  slight  but  constant  ^oLiounts  of 
impurity  were  picked  up  by  the  gas  from  the  flow  system. 

The  average  of  the  results  for  the  runs  at  the  three 
highest  for  carbon  dioxide  at  T = 291.5±0.9°K.  is: 

* P^avo  = 3.2x10-''  sec.  atm. 

(average  deviation  = iO/j) 

Upon  correcting  for  the  rate  expression  used  by  Kantrowitz 

as  noted  in  section  II  B, 

X p)  = 2.5  X 10"^  sec.  atm.  , 
ave 

73  4 

or  collision  number  = 2.2  x 10  . 

The  computation  of  a "collision  number"  from  ("t  x p)ave> 
using  the  mean  of  the  Initial  and  final  pressure  values,  is 
clearly  a first  approximation.  Refinement  of  the  kinetic 
theory  of  this  dependence  is  needed  to  deduce  an  appropri- 
ate weighting  function. 

The  relation  of  the  above  result  to  other  available 
data  can  be  seen  most  clearly  from  figure  1. 

It  should  be  pointed  out  that  it  is  not  simple  to 
compare  the  scatter  in  the  "collision  nuraber"  for  sonic 


Figure  16 
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experiments  with  that  for  our  method.  In  sonic  experiments 
one  does  not  usually  locate  the  point  of  maximujn  inflection 
of  a complete  velocity  versus  freupiency  curve  for  a series 
of  pressures.  Rather  one  substantiates  the  fact  that  the 

I 

reaction  is  of  second  order  by  showing  that  all  data  re- 
duce to  the  same  "s*‘  shaped  curve  when  the  velocity  squared 
is  plotted  versus  the  log  of  the  ratio  of  the  frequency 
to  the  pressure.  Strictly  spealting,  one  should  differenti- 
ate the  "order"  of  the  reaction,  which  indeed  has  been 
demonstrated  to  be  of  the  second  degree,  from  the  "molecu- 
larity"  which  is  inferred  from  the  inverse  pressure  depen- 
dence of 

This  work  on  carbon  dioxide  was  not  meant  to  be 
definitive,  but  rather  to  serve  as  a test  of  our  equip- 
ment and  to  develop  a definite  procedure  for  its  use. 

In  the  next  section  results  from  the  application  of  the 
impact  tube  experiment  to  a ralxture  of  reacting  gases  are 
presented. 


O 
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V.  Nitrogen  Tetroxide 

In  this  section  heat  capacity  lag  raeasurements  on 
nitrogen  tetroxide  are  recorded  and  analyzed.  LiLiitations 
of  the  impact  tube  method  are  then  discussed. 

A.  The  Data 

The  purity  of  the  gas  used  and  the  procedure  followed 
have  already  been  described  in  section  III.  Tl'ie  calibration 
plots  obtained  for  these  runs  were  similar  to  that  shown 
in  figure  11,  except  that  further  adjustment  to  the  null 
position  of  the  differential  transformer  with  respect  to 
its  core  increased  the  sensitivity  approximately  35/«.  A 
typical  plot  of  nitrogen  and  nitrogen  tetroxide  defects 
vs  expansion  ratio  is  shown  in  figure  17.  It  s)iould  be 
noted  that  for  a given  expansion  ratio  the  pressure  defect 
(uncorrected  for  background)  found  for  nitrogen  tetroxide 
is  approximate  1/3  that  for  carbon  dioxide.  Corrected  for 
background  the  ratio  becomes  nearly  1:4.  This  decrease  in 
magnitude  of  effect  is  only  partially  compensated  for  by 
the  increase  in  sensitivity  noted  above.  The  net  effect 
is  an  increase  in  the  per  cent  average  deviation  of  the 
nitrogen  tetroxide  results  •vdien  compared  to  those  for 
carbon  dioxide. 

In  contrast  with  the  carbon  dioxide  data  no  decline 
in  the  slope  of  the  pressure  defect  curve  above  intermediate 
expansion  ratios  was  noticed,  nor  were  sharp  breaks  found 
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in  the  curve  when  the  orifice  exit  velocity  of  the  gas 
approached  a Mach  number  of  unity.  The  elimination  of  the 
former  phenomenon  is  very  probably  due  to  the  use  of  the 
longer  end  more  carefully  machined  nozzle  which  is  des- 
cribed in  section  III.  The  absence  of  a sharp  break  in  the 
nitrogen  tetroxide  defect  curves  for  expansion  ratios 
around  two  is  possibly  a result  of  very  small  temperature 
inhomogeneities  due  to  operation  close  to  the  condensation 
temperature  of  the  gas.  The  measured  defect  did  increase 
quite  strongly  in  this  region  indicating  a marked  increase 
in  turbulence. 

Since  the  exact  point  at  which  this  turbulence  begins 
to  contribute  significantly  to  the  magnitude  of  the  ob- 
served defect  is  difficult  to  ascertain,  no  data  above 
(Pq/Pi ) “ were  used  to  determine  the  relaxation  time. 

In  general,  for  (Pq/Pi)  < 1*3  the  nitrogen  background  be- 
comes such  a large  fraction  of  the  defect  observed  for 

jj  — 

nitrogen  tetroxide  that  variations  in  this  background  are 
large  compared  to  the  difference  between  the  nitrogen 
tetroxide  defect  and  nitrogen  defect.  Hence,  values  below 
(Pq/Pi)  “ ^*3  were  also  omitted  from  consideration. 

From  a graph  of  the  type  shown  in  figure  17  values  of 
the  nitrogen  tetroxide  defect  minus  that  for  nitrogen  at 
selected  expansion  ratios  were  obtained.  Attempts  to  secure 
useful  data  at  significantly  less  than  0.3  atm.  were 
frustrated  by  the  smallness  of  the  difference  between  the 


defects  observed  for  nitrogen  tetroxide  and  those  observed 
for  nitrogen  compared  to  fluctuations  in  the  nitrogen  back- 
ground and  zero  point  reading.  The  results  of  measurements 
at  Pq  = 0.7,  0.5,  and  0.3  atm.  taken  from  smooth  curves 
through  the  corrected  pressure  defects  are  displayed  in 
table  II,  which  includes  data  for  two  different  tanks  and 
two  different  impact  tubes.  When  a run  was  repeated  under 
the  same  conditions  with  the  same  tank  of  gas  a given  set 
of  results  was  customarily  reproducible  to  within  ±10%. 

The  impact  tubes  used  had  the  following  character- 
istics: 

Tube  A: 

O.D.  = 0.0292  cm.,  I.D./O.D.  = 0.55, 
eccentricity  < l^j  . 

Tube  B: 

O.D.  = 0.0440  cm.,  I.D./O.D,  = 0.51, 
eccentricity  < 1%. 

For  the  purposes  of  easy  reference  each  of  these  eight 
runs  has  been  assigned  a roman  numeral  identification 
niimber . 

Runs  I,  II,  and  III  are  believed  to  be  representative, 
and  it  is  on  these  that  attention  will  be  primarily  focused 
In  the  following  analyses.  Runs  IV,  V,  and  VI  are  included 
to  illustrate  variations  due  to  sample  impurity  (tanks  N), 
while  runs  VII  and  VIII  are  included  to  ^illustrate  the  ef- 
fect of  changing  the  diameter  of  the  impact  tube. 


Table  II 


Pressure  Defect  Data  on  Nitrogen  Tetroxide 


Net*^  defect  in  cm.  n-butyl  phthalate 


Run 

Experiment 

No.  Tank 

Impact 

tube 

, £o 
\.atm)  p, 

= 1.30 

1.40 

1.50 

1.60 

1.70 

I 

8-29-52 

M 

A 

0.492 

0.54 

0.67 

0.30 

0.96 

1.13 

II 

9-1-52 

M 

A 

0.699 

O.SO 

0.86 

0.96 

1.11 

1.33 

III 

9-1-52(11) 

M 

A 

0.311 

0.31 

0 . 44 

0.5s 

0.73 

0.90 

IV 

9-4-52 

N 

A 

0.704 

1.14 

1.23 

1.33 

1.44 

1.56 

V 

9-4-52(11) 

N 

A 

0.502 

0.63 

0.75 

0.39 

1.06 

1.27 

VI 

9-4-52(111) 

N 

A 

0.307 

0.29 

0.37 

0.47 

n CQ 

w • > 

0.71 

VII 

9-16-52 

N 

B 

0.485 

0.24 

0.32 

0.43 

0.56 

0.73 

VIII 

9-18-52 

N 

B 

0.289 

0.23 

0.30 

CO 

• 

0 

0.49 

0.62 

4 

After  correction  for  nitrogen  baclc/^round 
(cm.  n-butyl  phthalate)/972. 8 = atn. 
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B.  Treatment  of  Data 

The  impact  tube  experiment  yields  values  of  the  pres- 
su} e defect  for  various  combinations  of  Pq,  Tq,  and  p, . 

These  data  have  been  used  to  determine  the  relaxation  time 
following  the  theory  developed  in  section  II.  The  experi- 
mentally determined  entropy  gain  is  compared  with  the  en- 
tropy gain  computed  for  instantaneous  compression  and, 
from  this  ratio,  the  ratio  of  the  unknown  relaxation  time 
of  the  gas  to  the  known  characteristic  compression  time  of 
the  equipment  is  determined.  Two  cases  will  be  considered; 
(i)  both  the  vibrational  degrees  of  freedom  and  the  chemi- 
cal reaction  lag,  and  (ii)  only  the  chemical  reaction  lags. 

A stepwise  description  of  the  method  of  computation  follows; 
tables  of  the  more  significant  quantities  computed  have 
also  been  included. 

1.  The  thermodynamic  quantities  needed  for  this  compu- 
tation are  the  heat  capacities  and  the  equilibrium 
constant  as  functions  --f  the  tempers ture. 
i)  For  NO2  the  fundamental  frequencies  used  to  calcu- 

I 

late  the  heat  capacity  are  - 648,  1320,  and 
1621  cm.'\  The  values  obtained  are: 

CpCNOj)  in  cal.  mole-*  degree"'  T°K. 

9.07  293.2 

9.03  293.2 


S.99 


288.2 


li)  For  U2O4  the  fundament?!  frequencies  used  to 

75 

calculate  the  heat  capacity  are  : = 1265, 

vij  = 1360,  7^3  = 752,  = 613,  = 1744, 

“ ? > *^8  ~ 283,  *^9  =*^11  ~ 330,  ~ 2 “ 

500  cm.-'.  2-^7  represents  the  torsional  mode 
and  is  assiuned  to  be  sufficiently  low  so  that 
the  equipartition  value,  H,  of  the  heat  capacity 
Is  reached  at  room  temperature.  The  values  ob- 
tained are: 


Cp(N204)  in  cal.  raole"’ degree"’ 
13.88 


X iV  • 


293.2 


13.73 
IS. 53 


00"3  O 


CO  • ^ 


The  heat  capacities  at  other  temperatures  have 
been  obtained  by  linear  interpolation. 


ill)  The  values  of  the  measured  ideal  equilibrium  con- 
• stants  (extrapolated  to  zero  pressure)  have  been 

76 

taken  from  the  paper  by  Verhock  and  Daniels. 


T°C. 


KgQ  (atm) 

0.1426 

0.3183 

0.6706 


The  logarithm  of  Kgq^  was  plotted  vs  the  recipro- 
cal of  the  absolute  temperature,  and  according  to 

the  equation 

in  K = - + ASl 

P HT  R 


109. 


a straight  line  was  drawn  through  the  points; 

AH°  and  were  calculated  from  its  slope  and 
intercept  respectively.  The  values  obtained  are: 
AH%8  = 14649  cal.  mole-’ 

AS%8  = 45.26  cal.  degree“’mole“’  . 

These  check  well  with  the  values  given  by  the 
authors;  they  were  recomputed  to  insure  that  5 
consistent  set  of  constants  were  being  used. 

2.  The  computation  of  the  flow  conditions  starts  with  a 
specified  set  of  values  for  p^,  Tq,  and  p, , to  deter- 
mine T,  and  u, , The  procedure  actually  foliowea  was 
for  a selected  (p^jT^)  to  compute  p,  and  u,  for 
sevor^sl  selected  T,  . A plot  of  T,  and  u,  v^  (Po/Pi) 
was  then  prepared,  and  values  for  any  given  expansion 
ratio  read  therefrom.  Values  of  p,  were  computed  from 
eqn.  (44).  written  in  the  form: 

(l+a,)  B in  Ea  + (l-.i,)  R in 

Pi  1+f-O 

+ in,  R in  . i^)  + (oo-o, ) H In  KpCT^) 

+ AS°(a,-aQ)  + |(l-a,)Cp^2)  ■*"  ^ 

= residue  . (87) 

The  value  of  Pi  consistent  with  reducing  the  residue 
to  zero  was  obtained  by  computing  the  residue  for 
several  p’s  very  close  to  pi  and  interpolating.  The 
quantity  (^^  u,^)  was  computed  from  (37),  rewritten  as 


no. 


&Ho,  = - I u," 

= AH^Coi-Oo)  + ^(l-a,)Cp^2j  + 2a,Cp^^j  (Ti-Tq)  . 

These  values  are  recorded  in  table  III. 

3.  One  now  must  compute  AS’.  For  small  changes  in  the 
entropy, 


^ . AP  + (#)  . (€) 


T,G,a 


p,9,a  <^P  9, a 


I)  (#)  AP  + (iS)  (^)  Ap  , 

9 p,T,a  T,c  p,T,9  ‘^P  T,0 


and  since  all  of  the  partial  coefficients  change  but 
slightly  for  small  changes  around  (pq,  T^,  a^) , to  a 
good  approximation 

ASOC  AP  . (3S) 
Then  one  may  write,  using  (80), 


AS'  = — '.  (39) 

(Po'Pzn.c. 

Values  for  (Po-Pa)  ^re  observed.  The  instantaneous 
compression  values  are  computed  from  equation  (41)  sJid 
tabulated  in  table  IV.  The  "corr.”  term  in  this  table 


is  precisely 

oorr.  5 + (1-c)  R In 

+ 2a,  R in  ■ -^)  . 

However,  to  a second  order  in  <<1, 

Oq  I 

may  be  approximated  by 


(90) 
this  term 


Table  III 

Summary  of  Co.;.ruted  Values: 

0 — 1 Transition,  T^  = 25°C. 


/ \ 
(atm. ) 

T,°C. 

/ N 

(atm. ) 

Pq 

Pi 

/ C C3.X  • \ 

S.iole' 

0 . 3000 

25 

0 . 3000 

1.0000 

0.3259 

0.0 

0.00 

22 

.2449 

1.2250 

.3189 

158,6 

1.20 

19 

.1993 

1.5053 

.3119 

317.1 

1.70 

16 

.1614 

l.S5o7 

. 3048 

475.7 

2.08 

0 . 5000 

25 

0.5000 

1.0000 

0.2530 

0.0 

0.00 

22 

.yj97 

1.2204 

.2518 

147.1 

1.16 

19 

. v34'^ 

1-4934 

. 2455 

293.6 

1.63 

16 

.2725 

1.8349 

.2391 

44'" . 5 

2.01 

0.7000 

25 

0.7000 

1.0000 

0.2201 

0.0 

0.00 

22 

. 5756 

1.2161 

.2144 

140.5 

1.13 

19 

.4720 

1.4330 

or\  o c 

• A.  W ^ 

230.4 

1.60 

16 

.3856 

1.8154 

.2027 

/v21.9 

1.96 

111. 


corr.  = -RCao-ci,)'  ||  • ^7^  " ^ 

T2  was  computed  from  equation  (33),  rewritten  as 
follows; 

AH, 2 = I u/  = {(l-“i)Cp(2)  + (^2-1,).  (92) 

In  table  IV  (Po~P2)i.c.  been  tabulated  for  two 
cases. 

i)  The  vibrational  degrees  of  freedom  and  the  chemi- 
cal reaction -both  lag;  then  CpCNOj)  = 4 

Cp(N204)  = 5 B,  and  is  obtained  by  subtract- 

r 

ing  Cp^'^v  from  . Note  that  the  torsional 

mode  (V7)  has  been  included  in  CpC.NzO^)  as  part 
of  the  non-lagging  heat  capacity.  This  follows 
from  the  comment  in  lii)  above. 

ii)  Only. the  chemical  reaction  lags,  in  which  case 

^p(i)  ~ *^p(i)  ^0(3)  ~ simplifies  the 

computation  somewhat  in  that  one  of  the  terms  in 
(41)  drops  out;  however,  one  must  also  note  that 

1 

Cp(i)  is  now  temperature  dependent. 

As'  was  computed  for  the  two  cases  according  to  (89) 
by  taking  * the  ratio  of  the  experimental  and  "instan- 
taneous compression"  pressure  defects.  Final  values 
have  been  recorded  in  table  V.  Figure  IS  shows 
(Po-P2)i  « vs  £2.  for  = 25°C.  and  various  p 's  for 

■*'*^*Pi  ^ 

the  case  "vdiere  both  the  vibrational  modes  and  chemical 
reaction  lag,  and  for  the  case  vhere  only  the  chemical 
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Table  V 

Evaluation  of  T for  both  vibrational  modes  and  chemical  reaction  lagging 
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Table  V (Con't.) 


Evaluation  of  "Xi  for  both  vibrational 
modes  and  chemical  reaction  lagging 


Impact  Tube  B 


Tank  N 


i 


( 


£o 

Pi 

( 

sec . 

(Po-P^^l.c. 

(gtm.xlO"*'^) 

v^sec.xlO  ) 

r’ 

r 

(sec.xlO"^^) 

Po  = 

0.30  atm. 

Run  VI 

II 

1.30 

1.37 

1.5 

3.21 

0.158 

0.210 

0.67 

1.40 

1.55 

2.3 

2.84 

0.134 

0.172 

0.49 

1.50 

1.69 

3.4 

2.60 

0.115 

0 . 144 

0.37 

1.60 

1.83 

4.5 

2.40 

0.112 

0.139 

0.33 

1.70 

1.94 

5.8 

2.27 

0.110 

0.3.37 

0.31 

Po  = 

0.50  atm. 

Run  VII 

1.30 

1.32 

2.4 

3.33 

0.103 

0.127 

0.42 

1.40 

1.49 

3.9 

2.95 

0.084 

0.102 

0.30 

1.50 

1.64 

5.6 

2.68 

0.079 

0.094 

0.25 

1.60 

1.77 

7.5 

2.4B 

0.077 

0.092 

0.23 

1.70 

1.S3 

9.5 

2.34 

0.079 

0.094 

0.22 

Evaluation  of  Z for  only 
reaction  lagging 

Impact  Tube  B 

chemical 

Tajil:  N 

Po 

( 

(Po-P2)i.c. 

(atm.xlO'*’^) 

( 5ec.xl0+**) 

AS’ 

lr' 

X 

(sec.xlO^^) 

Pq  ^ 

0.30  atm. 

Run  Vi I I 

1?30 

1.37 

0.7 

3.21 

0.338 

0.56 

1.30 

1.40 

1.55 

1.0 

2.84 

0.308 

0.49 

1.39 

1.50 

1.69 

1.5 

2.60 

0.260 

0.40 

1.04 

1.60 

1.33 

2.0 

2.40 

0.252 

0.^8 

0.91 

1.70 

1.94 

2.6 

2.27 

0.245 

0.37 

0.84 

Po  “ 

0.50  atm. 

Run  VII 

1?30 

1.32 

l.C 

3.33 

0.247 

0.37 

1.23 

1.40 

1.49 

1.7 

2.95 

0.194 

0.272 

0.80 

1.50 

1.64 

2.4 

2.68 

0.184 

0.254 

0.68 

1.60 

1.77 

3.2 

2.48 

0 . 180 

0.247 

0.61 

1.70 

1.83 

4.1 

2.34 

0.153 

0.253 

0.59 
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reaction  lags.  J'rora  these  graphs  values  of  (Po-Pz)^  r. 
may  be  read  at  selected  Pq/Pi  for  cOiUparison  with  the 
experinental  quantities.  A si;nilar  graph  was  used  to 
deduce  interpolated  values  of  u, . Figure  19  shows  the 
dependence  of  u^ , T^,  p, , a,,  and  T2  (vibrational 
modes  and  chemical  reaction  both  lag.^ing)  on  Pq/pi 
for  Pq  = 0.500  atm.  and  = 25°C. 

4.  To  compute  T*  from  AS'  one  uses  equation  (79).  For 
the  case  of  a square  ended  iijipact  tube  the  integral 
has  been  evaluated  by  Griffith  and  the  values  pre- 
sented at  the  end  of  section  II  B.  From  a graph  of 
these  values  plotted  as  log  'X'  vs  AS',  values  of  b' 
for  the  determined  values  of  AS'  have  been  obtained. 
These  are  presented  in  table  V.  Tq  was  found  from 
(57)  by  taking  the  ratio  of  the  measured  diameter  of 
the  impact  tube  to  the  computed  streaming  velocity  in 
state  1;  then  9:  was  computed  according  to  (68). 

Values  of  and  "V  arc  presented  in  table  V. 

This  completes  the  reduction  of  the  measured  pressure 
defects  to  relaxation  times  according  to  the  theory  out- 
lined in  section  li  C. 

0 


I 


i 
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C.  Discussion 

a.  Interpretation  of  the  Data. 

If  the  transfer  of  energy  is  bimolecuiar  and  only  one 
relaxation  time  is  operative,  the  product,  p.  x 1:  is  the 

average  pressure,  and  is  taken  to  be  ^(Pq+Pi)]  , ‘/hioh  is 
proportional  to  the  •'collision  i.uiiiber"  should  remain  con- 
stant as  long  as  the  temperature  and  composition  are  not 
varied.  In  table  VI  values  of  p and  (p  x Y)  ere  displayed 
for  runs  I,  II,  and  III,  where  T has  been  computed  for  both 
vibrational  modes  and  chemical  reaction  lagr.ing  and  for 
only  the  chemicnl  reaction  lagging.  Figure  10  Is  a j-lot 
of  p X "I:  vs  p for  these  two  cases.  If  one  of  the  procedures 
used  to  reduce  the  measured  pressure  defects  to  relaxation 
times  is  completely  correct  then  for  that  case  all  of  the 
points  in  figure  20  should  fall  on  the  same  horizontal 
line.  Clearly  this  is  not  the  case.  One  notes,  however, 
that  when. reduced  on  the  basis  of  both  vibrational  modes 
and  chemical  reaction  lagging  the  points  cover  a co)isider- 
ably  narrower  range  than  when  reduced  on  the  basis  of  only 
the  chemical  reaction  lagging.  This  argues  in  favor  of 
the  inclusion  of  both  types  of  lag  in  the  reduction  of 
these  data.  One  further  notes  from  table  VI  (or  figure  20) 
that  there  is  a strong  trend  in  the  values  of  T cora,jUted 
for  any  one  run.  This  dependence  cannot  be  due  to  the 
pressure,  since  p varies  by  at  raost  15%  within  a run.  ilor 
can  it  be  due  to  temperature  variations  since  T .=  -^(Ti+Tj) 


Evalucition  of  p 


III,  Po  = 0.30 
atm. 


1.30 

1.40 

1.50 
1 .60 
1.70 


I,  = 0.50 
atia. 


1.30 

1.40 

1.50 

1.60 

1.70 


II,  r = 0.70 
o . 
atm. 


1.30 

1.40 

1.50 

1.60 

1.70 


Table  VI 


p X Y for  Huns  I,  II,  and  III 


p 

P 

X for  both 

P X for  only 

6 

and  n lagging 

a lagging 

0.P66 

0.17 

0.52 

0.;*57 

0.13 

0.43 

0.250 

0.10 

0.31 

0 . 244 

0.09 

0.2b 

0.233 

0.03 

0.22  ! 

1 

0.442 

0 . 33 

1.35 

0.429 

0 . 22 

0.60 

0.417 

0.14 

0.42 

0 . /^0b 

0.11 

0.33 

0.397 

0.10 

0.27  1 

0.619 

0.50 

|i 

2.20 

0.600 

0.26 

0.92  j 

0.583 

0.16 

0.57  1 

0.569 

0.13 

0.44 

0.556 

0.11 

0.36 

II,  III) 


Figure  20 

^ 3C  p xa  p for  Tank  M (Runs  I, 

O Calculated  on  basis  of  both  vibrational 
modes  and  chemical  reaction  lagging. 

• Calculated  on  basis  of  only  chemical 
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varies  by  at  most  2%  (note  that  lower  T, ' s are  balanced  by 
accompanying  higher  Tz’s)  within  a run.  This  variation  is 
reduced  when  one  considers  both  vibrational  inodes  and 
chemical  reaction  lagging,  which  suggests  that  further 
modifications  of  the  simplified  procedure  used  to  reduce 
the  measured  pressure  defects  to  relaxation  times  are  in 
order.  The  next  step  would  be  to  consider  three  different 
relaxation  times,  as  initially  set  up  in  section  III  C. 

In  this  connection  it  is  valuable  to  note  that  a similar 
trend  was  noted  in  our  measurements  on  carbon  dioxide.  It 
is  prooable  that  the  use  of  three  different  relaxation  times 
associated  with  the  different  vibrational  frequencies 
would  also  effect  a ir.ore  complete  reduction  of  the  data 
in  that  case.  Further  conclusions  on  this  question  may 
be  derived  from  table  V.  It  is  clear  that  the  only  effect 
of  changing  the  diameter  of  the  impact  tube  is  to  change 
Having  allowed  for  this  factor  the  relaxation  times 
computed  for  all  impact  tubes  should  be  identical  Com- 
parison of  the  results  of  runs  at  identical  values  of  p^ 
and  Pq/Pi , and  with  the  same  tank  of  gas,  leads  to  the 
following  ratios; 
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Po 


Eo 

Pi 


X deteriiilned  vlth  impact  tube  A 
determined  with  impact  tube  B 


0.30  atm.  1.30  0.90 

1.40  0 . S6 

1.50  0.86 

1.60  0.32 

1.70  0.77 

0.50  atm.  1.30  2.16 

1.40  1.S8 

1.50  1.56 

1.60  1.35 

1.70  1.23 


Clearly  the  ratio  of 't’ s for  the  two  impact  tubevS  depends 
on  p.  This  again  indicates  that  two  or  more  relaxation 
times  are  necessary  to  properly  reduce  the  data. 

From  figure  20  an  estimate  of  the  avorai^e  relaxation 
time  may  be  made  by  drawing  the  best  horizontal  line  through 
those  points  calculated  on  the  assumption  that  both  the 
vibrational  modes  and  chemical  reaction  lag.  At  2'5^C,  and 
1 atm.  total  pressure,  this  gives 

'Vs  0.14  M-  sec , I 


b.  Sources  of  Error 

The  above  result  is  sub^lect  to  four  primary  sources 
of  error. 

i)  Since  reduction  of  the  data  is  not  complete  the 
above  value  represents  an  average.  The  average  deviation 
from  this  value  including  all  points  is  0.07.  The  elimi- 
nation of  one  point  whose  deviation  is  five  times  this 
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value  gives  an  average  deviation  of  0.05.  This  is  the 
single  largest  error  involved  in  this  experiment.  How- 
ever, reduction  of  this  error  depends  only  upon  the  use 
of  a more  complex  theory  for  converting  the  measured 
pressure  defects  to  relaxation  times.  The  very  regularity 
of  the  trend  for  each  run  in  figure  20  augurs  well  for 
the  success  of  a more  complete  theoretical  treatment. 

ii)  The  effect  of  impurities  is  indicated  in  table 
V by  data  for  comparable  runs  on  two  different  tanks.  The 
average  difference  between  X values  for  the  two  tanks  is 
approximately  however,  the  ratio  of  the  value  of 

obtained  with  tank  N compared  to  tank  M varies  from  less 
to  greater  than  unity,  depending  on  p^.  In  part  this  may 
be  involved  with  the  theoretical  problem  referred  to  above, 
or  more  likely  to  a greater  dilution  of  a fixed  amount  of 
impurities  desorbed  from  the  flow  system  at  higher  Pq's* 

The  purity  specifications  on  the  gas  as  supplied  by  the 
manufacturer  (see  section  III)  indicate  that  the  gas  con- 
tains only  very  small  amounts  of  impurities  which  might  be 
expected  to  be  troublesome.  It  should  also  be  noted  that 
before  use  the  system  was  always  swept  clean  with  "dry” 
nitrogen  and  before  data  were  tcdcen  the  nitrogen  tetroxide 
was  run  through  the  system  until  at  a fixed  Pq  and  Pq/pi 
a constant  pressure  defect  was  observed.  It  seems  improb- 
able that  impurities  would  be  present  in  the  flow  system  in 
such  amount  and  desorb  in  such  a manner  as  to  produce  such 


117. 


a constant  Pleasured  defect  for  any  appreciable  ti.ne.  It 
was  also  found  possible  at  the  end  of  a run  to  return  to 
any  stipulated  set  of  conditions  and  obtain  very  nearly 
(see  point  iv  below)  the  sane  value  for  the  pressure  de- 
fect. All  of  these  factors  argue  against  the  prcbcbility 
of  gross  impurities  being  introduced  by  the  flow  system, 
and  indicate  that  the  impurities  in  the  gas  as  supplied  by 
the  manufacturer  vary  from  tank  to  tanlc.  The  absolute 
magnitude  of  the  error  in  'f  caused  by  these  impurities  is 
impossible  to  ascertain,  but  it  is  probably  less  than  a 
factor  cf  two.  The  error  from  this  source  is  not  inherent 
in  the  experiment.  Its  removal  depends  only  on  obtaining 
adequately  purified  nitrogen  tetroxide  in  sufficiently 
large  amounts. 

iii)  An  error  is  introduced  into  the  computations  due 
to  the  lack  of  precise  values  for  the  necessary  thermo- 
dynamic parameters.  The  magnitude  of  the  error  introduced 
in  this  manner  may  be  investigated  by  co;r.puting  the  vaiue 

of  (£o)  using  equation  (41)  for  various  assumed  values 
Pa  i- c . 

of  the  thermodynamic  constants.  One  finds  for  the  r.inge 

of  values  used  in  this  thesis  that  a 0.1%  change  in 

changer,  the  value  of  (-•*?-)  by  approximately  10%.  A 

Pz  i.c. 

change  of  0.1%  in  AS°  also  causes  a change  in  (— ) of 

P2  i . c . 

this  same  order  of  magnitude.  These  pararoeters  are  known 
to  approximately  0.1>^  and  their  method  of  computation  (see 
section  V B)  assures  the  naximum  of  self-consistency.  Hence, 
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one  nay  conclude  that  lack  of  sufficiently  accurate  know- 
ledge of  the  thernodynanic  parameters  introduces  an  error 

into  the  computed  value  of  (— ) of  roughly  + 10%.  The 

P2  i.c. 

error  subsequently  Introduced  in  the  magnitude  of  T is  of 
this  order  of  magnitude.  It  should  be  noted  that  this  is 
not  an  error  in  the  data  obtained  from  this  experirsent, 
but  reflects  the  lack  of  adequate  knowledge  of  the  parameters 
necessary  to  interpret  the  experimental  data.  Ihat  the 
interpretation  of  the  experiment  depends  so  strongly  on 
the  v^alues  of  these  thermodynamic  parameters  is  not  sur- 
prising in  view  of  the  fact  that  ^ (po’P^'i.c.  ^ 

In  the  experiment  (Pq-Pj)  is  determined  directly  by  anply- 
ing  the  two  pressures  to  different  sides  of  the  pressure 
sensing  bellows.  The  calculated  values  are  obtair'cd  by 
subtracting  two  large  numbers,  and  an  error  of  0.01,il  in 
the  computed  value  of  pj  will  change  the  computed  value  of 
(Pq-P2>  by  roughly  10:v. 

iv)  The  errors  inherent  in  the  present  equipment  are 
those  due  to  lack  of  precision  in  reading  the  pressure, 
turbulence  background,  and  zero  point  wandering.  All  of 
these  have  been  considered  in  detail  in  section  111.  As 
indicated  there,  pressure  difft.rentials  can  be  road  to  a 
precision  of  1.0x10-'  atm.  Reference  to  figure  17  shows 
that  the  nitrogen  defect  over  the  region  used  amounted  to 
approximately  2.0x10-“^  atm.  Assuming  that  this  defect 
was  the  same  for  the  nitrogen  tetroxide  part  of  the  run 
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probably  involves  an  error  of  roughly  10^^  of  this  quan- 
tity, or  2xi0~*.  Errors  introduced  by  wandering  of  the 
zero  point  were  considerably  reduced  by  frequent  readings 
of  this  quantity.  The  maximum  value  of  this  uncertainty 
as  actually  observed  over  many  experiments  is  ±0.2  mv. 
at  the  recorder  or  ± 5x10*^  atm.  Hence  the  maximum  error 
in  determining  the  corrected  pressure  defect  of  N2O4  is 
± 8x10-^  atm.  Oince  a typical  corrected  pressure  defect 
(see  tabic-  II)  amounted  to  0.5x10-^  atm.  this  is  equiva- 
lent to  a.n  error  of  ± As  discussed  in  section  III, 

the  largest  part  of  the  error  which  is  due  to  zero  point 
wandering  could  be  greatly  reduced  by  the  use  of  more 
stable  electrical  components.  A reduction  of  this  source 
of  fluctuations  to  20%  of  its  present  value  would  increase 
the  precision  of  the  measurements  by  a factor  of  2,  i.e., 
to  ± S%. 

One  further  point  should  be  mentioned  with  respect  to 
errors,  and  that  is  the  possibility  of  entropy  gain  in  the 
nozzle.  The  nozzle  utilized  for  the  experiments  on  nitro- 
gen tetroxide  has  a length  of  1.6  inches,  and  for  a typi- 
cal exit  velocity  of  1.6x10^  cm. /sec.  the  time  taken  by 
the  gas  to  expand  through  the  nozzle  is  greater  than  250 
li  sec.  Comparing  this  with  the  value  of  X quoted  above 
the  time  taken  to  expand  through  the  nozzle  is  approximate- 
ly 1700  times  the  relaxation  time.  Clearly  the  entropy 
gain  in  the  nozzle  due  to  this  relatively  slow  expansion 
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must  be  neg] igible  compared  to  the  entropy  gain  on  com- 
pression at  the  face  of  the  impacc  tube. 

c.  Comparison  with  other  Methods 

As  noted  in  section  I,  a number  of  methods  have  been 
suggested  for  determining  the  rate  of  dissociation  of 
nitrogen  tetroxide,  but  only  one,  involving  the  use  of  a 
shock  tube,  had  any  pronounced  success.  Ihe  bimolecular 
and  unimolecular  rate  constants  determined  by  Carrington 
and  Davidson  using  this  method  have  been  quoted  at  the 

65 

end  of  section  I.  Manes  discussed  the  relationship 

between  the  relaxation  time  and  the  postulated  reaction 

mechanism  for  a general  case  w^iich  includes  a large  number 

6ib 

of  simple  rate  laws.  He  obtained  the  result  that  for  a 
reaction  of  the  type 

aA  + bB  + • • ♦ qQ  + rR  + • • • 
near  equilibrium,  for  most  cases  'V  is  given  by 

> >*ere 

C^  is  the  concentration  of  species  A and  r^.  is  the  rate 
of  the  forward  reaction  as  written  above.  Utilizing  the 

kf 

relation  K = r—  and  applying  the  above  formula  to  the 
unimolecular  case 

2NO2 
^r 

^ = kf  + kj.  40^02  • 


gives 
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Substituting  Carrin^-ton  and  Davidson's  unimolecular  rate 
constant  Tor  kj^  gives  'b  - 0.022  usee.  For  the  hirnolecular 
case 

M + 2NO2  + M 

kj, 

^ kj,  4 Cj^  CfjOj  + kf  (Cj^  + 2 and 

using  Carrington  and  Davidson's  binoleculr.r  rate  constant 
for  kf  gives  '^  = 0.99  usee.  One  can  easily  show  that  the 
condition  for  deactivation  by  collision  or  chemical  decoir- 
position  to  be  equally  probable  is  that  the  tota]  concen- 
tration of  molecules  satisfies  the  equation 

C»  = 1 . 

^unimoi. 

Substituting  in  the  values  of  the  rate  coristants  used  above 
.gives  Crp  = 0.5  This  checks  as  well  as  can  be  exT'ected 
with  Carrington  and  Davidson' s conclusion  based  on  a plot 
of  exp.erimentai  values  of  log  , vs  log  C^,  that  the 

reaction  is  nearly  second  order  at  0.05  i.e.  when  the 
total  concentration  is  decreased  to  one- tenth  of  C^. 

One  atmosphere  corresponds  to  a total  concentration 
of  0.04  and  hence  all  of  the  data  reported  in  this 
thesis  were  taken  in  the  range  vdiere  the  bimolccalar  rate 
law  applies.  It  should  be  noted  then  that  the  average 
X deduced  from  the  Impact  tube  experiment  is  approximately 
one- seventh  that  calculated  from  the  bimolecular  rate 
constant  based  on  the  ;hock  tube  •'xperiments . The  differ- 
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ence  may  be  interpreted  as  being  due  to  the  greater  ef- 
ficiency of  collisions  vith  nitrogen  dioxide  and  tetroxide 

I 

in  exciting  the  nitrogen  tetroxide  molecule  over  collisions 

with  nitrogen  which  was  used  as  the  carrier  gas  in  the 

shock  tube  experiments.  The  conclusion  is  that  the  two  * 

experiments  agree  in  the  area  of  comparison. 

The  most  important  consequence  of  these  two  experi-  i 

ments  is  that  both  definitely  indicate  the  presence  of  a 
lag  in  equilibration.  This  should  be  contrasted  with  the 
acoustic  method  which  has  led  to  many  conflicting  reports 
even  as  to  whether  there  was  or  was  not  a velocity  change 
(see  section  I) . 

I 

I 

It  is  interesting  to  compare  the  shock  tube  and  the  i 

impact  tube  methods  for  measuring  the  rates  of  rauid  i 

chemical  reactions  in  the  gas  phase  since  such  a compari-  j 

! 

I 

son  most  clearly  indicates  the  advantages  and  disadvantages  ! 

of  each. 

The  success  of  the  impact  tube?  method  is  due  to  its 
being  a differential  experiment.  The  shock  tube  method 
relies  on  the  rapidity  with  which  a shock  produces  a large  j 

change  in  energy.  This  points  up  the  primary  difference 
between  these  two  methods.  The  resulting  difference  in 
theoretical  treatment  has  been  partially  indicated  in 
section  II  C.  A further  difference  results  from  the  fact 
that  since  the  displacements  from  equilibrium  are  small  in 
the  impact  tube  experiment  approximations  can  bo  made  vdiich 
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are  not  applicable  to  the  large  displaceraents  of  the  shock 
tube  experiment.  On  the  other  hand,  the  amount  one  can 
learn  about  the  mechanism  from  such  small  displacements 
from  equilibrium  us  ultimately  less  than  chat  which  can 
be  deduced  from  a larger  range  of  perturbations. 

In  general,  the  theoretical  treatments  of  both  ex- 
periments are  hampered  by  a lack  of  adequate  hydrodynaj:ic 
theories.  In  the  impact  tube  experiment  this  manifests 
itself  in  the  inaccuracies  involved  in  converting  the 
measured  pressure  defects  into  relaxation  times.  These 
inaccuracies  have  already  been  discussed.  In  the  shock 
tube  experiment  this  difficulty  manifests  itself  in  the 
lack  of  ability  to  compute  precisely  the  conditions  behind 
the  shock  wave.  Carrington  and  Davidson  have  estiniated 
that  an  error  of  0.6)2  in  the  shock  velocity  means  an  error 
of  1.5°  in  the  calculated  temperature  and  an  error  of  13% 
in  the  rate  constant.  Even  if  the  shock  velocity  were 
knovn  more  accurately  other  workers  have  pointed  out  that 
simple  theories  do  not  predict  the  density  fluctuations 

77 

with  much  accuracy. 

In  the  shock  tube  experiment  the  limit  of  resolution 
is  set  primarily  by  the  time  constant  of  the  detector- 
amplifier.  One  is  thus  limited  to  relatively  strong  ab- 
sorptions in  the  visible  and  ultra-violet.  The  resolving 
time  js  then  limited  to  ~4  P sec.,  as  determined  by  the 


time  I'equired  for  the  shock  to  cross  the  observing  beam. 


This  could,  of  course,  be  reduced  by  using  narrower  slits 
if  light  sources  of  sufficient  intensity  were  available. 
This  limitation  prevented  Carrington  and  Davidson  from 
making  measurements  at  temperatures  greater  than  30'’C. 
since  the  reaction  became  too  fast.  The  limit  of  resolu- 
tion of  the  impact  tube  type  of  experiment  is  set  by  the 
minimum  diameter  of  the  impact  tubes  which  can  be  made 
and  the  maximum  non- turbulent  flow  velocity.  Impc.ct  tubes 
can  be  constructed  as  small  as  0.01  cm.  while  the  maximum 
useful  velocity  is  approximately  2x10^  cm. sec.*'.  This 
gives  a time  constant  of  ^^0.5  .'i  sec.  The  relaxation  time 
must  be  within  a factor  of  iO  of  this  for  an  appreciable 
defect  to  be  observed,  and  hence  the  limit  of  time  resolu- 
tion of  this  type  of  experiment  is  ■^0.05  n sec.  Relaxa- 
tion times  as  long  as  10“^  sec.  may  be  determined  by  study- 
ing the  entropy  gain  in  the  nozale. 

One  serious  limitation  on  the  impact  tube  method  is 
that  it  requires  large  quantities  of  gas.  The  shock  tube 
method  does  not  reauire  such  large  volUiUes.  Because 
trace  impurities  may  greatly  effect  the  relaxation  time 
the  gas  used  must  be  quite  pure,  and  with  the  large  volumes 
involved  in  the  impact  tube  experiment  this  restriction 
presents  a difficult  though  not  Insoluble  proble.p..  There 
is  also  the  possibility  of  working  out  a suitable  recycling 
technique. 


c 


The  shock  tube  'nethod  can  also  be  used  over  a some- 
what wider  range  of  pressures,  temperatures,  r-.nd  concen- 
trations than  the  impact  equipment  which  has  been  built 
to  date.  Further,  the  observed  uressure  defect  in  the 
impact  tube  experiments  depends  on  an  integral  of  the 
form 


CD 


T 


(HI) 


> 


vrtiereas  the  shocK  tube  method  hrs  u decldtu  advantage  in 
that  the  concentration  oi'  the  reoctants  is  measured  direct- 
ly as  a function  of  time.  To  its  disadvantage  is  the 
relative  complexity  of  the  equipment  required  compared  to 
that  utilized  in  the  impact  tube  experiment. 
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Further  Analysis  of  CO 2 and  N2O4 
Heat  Capacity  Lag  Data 


From  the  theoretical  analysis  presented  in  section 
VC.  two  conclusions  follow:  (a)  The  available  d?.ta  may 

be  interpreted  in  terms  of  a single  relaxation  time  only 
to  a rough  approximation;  three  and  possibly  more  re- 
laxation times  are  in  evidence,  (b)  A general  treatment 
for  the  case  of  plural  relaxation  times  has  been  outlined. 
Hence,  before  continuing  with  exj.eriraental  work  in  which 
this  method  is  applied  to  other  homogeneous  fast  gas 
reactions,  the  data  recorded  in  this  report  will  be  re- 
viewed with  the  object  of  deducing  values  for  the  several 
relaxation  times.  The  degree  of  success  which  we  shall 
meet  will  determine  whether  extensive  further  development 
of  the  (isentropic  flow  + impact  tube)  method  will  be 
undertaken. 


The  distribution  list  may  be  found  at  the  end  of 


Section  I. 
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APPENDIX  I 


Consider  a system  in  equilibrium  and  having  one  re- 
laxation tine.  Suppose  to  this  system  one  applies  a step 

, . 6 hg  for  t > 0 

perturbation  hg(t;  = s q fQi-  t ^ 0 * define  a 

parameter  '^^gCt)  which  is  a measure  of  instantaneous  de- 
parture from  equilibrium,  ihe  step  perturbation  imposes 

o 

an  initial  departure  from  equilibrium  of 


That  is  the  ratio  of  departure  from  equilibrium  to  per- 
turbation, r,  is 


o 


The  decay  curve  will  then  be  described  by  the  differential 
equation 


_ ns 

“dt“  * - 


where 


is  the  relaxation  tine  of  the  system.  The  integrated 
form  of  this  equation  will  be 

= V3  . 

If  one  now  assumes  that  the  departure  is  additive  in 
terms  of  the  perturbation,  i. e. . independent  of  previous 


1 r* 

X ««>  / • 


perturbation  history  a given  perturbation  always  produces 
the  so.'.ne  departure,  the  general  case  may  easily  be  treated. 
Suppose  h(t’)  has  a general  form  as  shown  below.  It  is 


clear  that  assuming  the  additive  property  this  may  be  iooK- 
ed  upon  as  a series  of  step  perturbations  spaced  in  time 
t’ . Over  a time  interval  dt*  aroimd  t’  the  additional 
perturbation  will  be 


This  additional  pertur- 


bation will  contribute  to  the  departure  from  equilibrium  at 
all  times  subsequent  to  t’  and  the  contribution  will  amount 
to 


L c^t» 


]<^f 


Hence  in  general  one  may  write 


-(t-f)/r 

e 


t’=t 


where  one  has  summed  over  all  the  perturbations.  This  nay 


be  rewritten  as 
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'^(t) 


If  one  now  further  assumes  that  r is  independent  of  time 
the  above  nay  be  written  as 


’i(t) 


o 

/'ns 


I -ii:) 

“s 


<^(  t’ ) 
" c^t» 


dt»  . 


For  the  case  where  the  perturbation,  although  extended  in 
tire,  starts  with  a step  function,  the  last  equation  assumes 
the  form 


o 


f Mt<) 


'^(t)  = e 
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